Fermi gases in one dimension: From Bethe Ansatz to experiments* 
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This article reviews theoretical and experimental developments for one-dimensional Fermi gases. 
Specifically, the experimentally realized two-component delta-function interacting Fermi gas - 
the Gaudin-Yang model - and its generalisations to multi-component Fermi systems with larger 
spin symmetries. The exact results obtained for Bethe ansatz integrable models of this kind 
enable the study of the nature and microscopic origin of a wide range of quantum many-body 
phenomena driven by spin population imbalance, dynamical interactions and magnetic fields. 
This physics includes Bardeen-Coopcr-Schricffcr-likc pairing, Tomonaga-Luttingcr liquids, spin- 
charge separation, Fulde-Fcrrcl-Larkin-Ovchinnikov-like pair correlations, quantum criticality and 
scaling, polarons and the few-body physics of the trimer state (trions). The fascinating interplay 
between exactly solved models and experimental developments in one dimension promises to yield 
further insight into the exciting and fundamental physics of interacting Fermi systems. 
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I. INTRODUCTION 

Fundamental quantum many-body systems involve the 
interaction of bosonic and/or fcrmionic particles. The 
spin of a particle makes it behave very differently at 
ultracold temperatures below the degenerate tempera- 
ture. There are thus fundamental differences between 
the properties of bosons and fermions. However, as 
bosons arc not subject to the Pauli exclusion princi- 
ple, they can collapse under suitable conditions into the 
same quantum groundstate - the Bose-Einstein conden- 
sate (BEC). Remarkably, even a small attraction between 
two fermions with opposite spin states and momentum 
can lead to the formation of a Bardeen-Cooper-Schricffcr 
(BCS) pair that has a bosonic nature. Such BCS pairs 
can undergo the phenomenon of BEC as temperature 
tends to absolute zero. Over the past few decades, ex- 
perimental achievements in trapping and cooling atomic 
gases have revealed the beautiful and subtle physics of 
the quantum world of ultracold atoms, see recent re 
view articles dBloch et all 120081 120121: IChin et alJ.l201C: 



One-dimensional 
gas tube 



Dalfovo et allll999tlGiorggini et aZj,[2008l:lLeggetfl.l2001 
Lewenstein et adl2007tlRegal and Jirj , l2006tlZhall2009l) . 

In particular, recent experiments on ultracold bosonic 
and fermionic atoms confined to one dimension (ID) have 
provided a better understanding of the quantum sta- 
tistical a nd dynamical e f fects i n quantum many-bod y 
systems (jCazalilla et all l201ll ; lYurovskv efall 12008th 
These atomic waveguide particles are tightly confined in 
two transverse directions and weakly confined in the axial 
direction. The transverse excitations are fully suppressed 
by the tight confinement. As a result the trapped atoms 
can be effectively characterised by a quasi-lD system, 
see Fig. [TJ The effective ID inter-particle potential can 
be controlled in the whole interaction regime. In such a 
way, the ID many-body systems ultimately relate to pre- 
viously considered exactly solved models of interacting 




Fully paired Partially polarized 
wing core 



FIG. 1 Experimental confinement of two-comp onent ultra- 
cold 6 Li atoms trapped in an array of ID tubes ([Liao et all 
l2010h . The system has spin population imbalance caused by 
a difference in the number of spin-up and spin-down atoms. 
From Bloch (2005). 



bosons and fermions. This has led to a fascinating in- 
terplay between exactly solved models and experimental 
developments in ID. Inspired by these developments, the 
study of integrable models has undergone a renaissance 
over the past decade. Their study has become crucial 
to exploring and understanding the physics of quantum 
many-body systems. 



A. Exactly solved models 

1. The virtuoso triumphs of the Bethe ansatz 

The study of Bethe ansatz solvable models began when 
Bethe (1931) introduced a particular form of wavefunc- 
tion - the Bethe ansatz (BA) - to obtain the energy eigen- 
spectrum of the ID Heisenberg spin chain. After laying 
in obscurity for decades, the BA emerged to underpin 
a diverse range of physical problems, from superconduc- 
tors to string theory, see, e.g., Batchelor (2007). For 
such exactly solved models, the energy eigenspectrum of 
the model Hamiltonian is obtained exactly in terms of 
the BA equations, from which physical properties can 
be derived via mathematical analysis. From 1931 to the 
early 1960s there were only a handful of papers on the 
BA, treating the passage to the thermodynamic limit and 
the ex tension to the anisotropic XXZ Heisenberg spin 
chain (Ides Clqizeaux and Pearsonl 19621 Griffiths! . Tl964t 



iHulthenl . Il938t lOrbachl Il959t IWalkerl . Il959h . Yang and 
Yang (1966a) coined the term Bcthc's hypothesis and 
proved that Bethe's sol ution was indeed the gro undstate 
of the XXZ spin chain ( Yang and Yand . Il966al lbllcl) . 
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The next development was the exact solution of the 
ID Bose gas with delta-function interaction by Lieb and 
Liniger (1963), which continues to hav e a tremendous 
impa ct in quantum statistical mechanics ( Cazalilla et all 
1201 lh . They diagonalised the hamiltonian and derived 
the groundstate energy of the model. This study was 
further exte nded to the excitations above the groundstate 
(iLiebl . I1963D . McGuire (1964) considered the model in 
the context of quantum many-body scattering in which 
the condition of non-diffractive scattering appeared. 

Developments for the exact solutio n of the I D Ferm i 
gas with d elta-function interaction ( Gaudinl . Il967allbl ; 
lYand . Il967l) are discussed in the next subsection. A 
key point is Yang's observation ( Yand . fl967|) that a gen- 
eralised Bethe's hypothesis works for the fermion prob- 
lem, subject to a set of cubic equations being satisfied. 
This equation has since been referred to as the Yang- 
Baxter equation (YBE) after the name was coined by 
Takhtadzhan and Faddeev (1979). Baxter's contribution 
was to independently show that such relations also ap- 
pear as conditions for commuting transfer matrices in 
two-dimensional latt ice models in statistical mechanics 
(jBaxterl . [l972al . Il982h . Moreover, the YBE was seen as a 
relation which can be solved to obtain new exactly solved 
models. The YBE thu s became celebrated as t he mas- 
ter key to integrability ( Au-Yang and Perkl . 1 19891 ). 

The study of Yang-Baxter integrable models flour- 
ished in the 70's, 80's and 90's in the Canberra, St 
Petersburg, Stony Brook and Kyoto schools, with far 
reaching implications in both physics and mathemat- 
ics. During this period the YBE emerged as the un- 
derlying structure behind the solvability of a number 
of quantum mechanical models. In addition to the 
XXZ sp i n chain , examples include the XYZ spin chain 
( Baxter!. 1972bl). the t— J model at supersymmetric cou- 
pling ([Essler and Korepinl . Il992fc iFberster and Karowski , 



1993libr and the Hubbard model (Essler et 
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tribal 



2005 



1968: 
I1972D . 



Frahm and Korepinl. Il990l Il99l| iLieb and 
Ogata and Shibal . Il990i IShastrvl . Il986allbl : IShib; 
Three collections of key papers have been publish ed 
(|jimbol . 119901: iKorepin and Esslerl . Il994t iMattisl . fl993l) . 

strong l y correlated elec- 
ISchollwockl 120041: 
spi n exchange inter- 



Further exa mples are 
tron syst e ms ( Giamarchi 



ITakahashil. Il999t iTsvelikl 



2005; 



200 
19951) 



Montorsil . Il992l ; ISutherlandl 



actio n ([Essler et al. 

l2004h . Kondo physics of quantum impurit ies coupled to 
conduction electrons in equ ili brium ( Andrei et all 

19831) and o ut of 

Mehta and Andrei 120061; 



Il983t iTsvelik and Wiegmann 
equilibrium (Dovonl |200 



iNishino and Hatano 




Nishino et oil. l2009h. the 



BCS model dCambiaggio et all 119971; iDukelskv et < 
2004 iDunning and Linksl l2004t iLinks et all 12005 : 



Ric hardson! 
196 



1963a. b, 1965; Ri chardson and Sherman 



, .von Delft and Ralph!. 1200 1 | ). m o dels wit h lonp 
range interact i ons ( CalogeroT " 19691; Gaudinl. 11976 



iHaldand . Il988t IShastrvl. 119881 ISutherlandl. Il97ll). tw o 
Josephson coupled BECs (IZhou et all l2002l l2003a[) . 
BCS-to-BEC crossover (jOrtiz and Dukelskvl . 120051) . 



atomic-mol e cular B ECs ( Foerster and Ragoucvl l2007t 
IZhou et all . l2003b|) and quantum degenerate gases of 
ultracold atoms ( Cazalilla et al. , 2011 [_ Korepin et all 



Il993t IPethick and Smithl . 120081: lYurovskv et all l2008h . 

A significant development in the the ory of quantum 
integrable systems is the algebraic BA dFaddeevl . Il984l 
iKulish and Sklvaninl Il982t ISklvanin et aUll979l) . essen- 
tial to the so called Quantum Inverse Scattering method 
(QISM), a quantized version of the classical inverse scat- 
tering method. The QISM gives a unified description 
of the exact solution of quantum integrable models. It 
provides a framework to systematic ally constru c t and 



solve quantum many-body systems (lEssler et al. 20051 
IKorepin et all . 119931 ITakahashil . Il999l; iThackerl ll98lh V 
Other related threads are the quantum transfer 



matr ix (QTM ) (iDestri and de Vegal Il992t iKliimper 



19921; iSuzukil . Il985l) and T-systems (jKuniba et al. 
1994aM l 201 ll ) from which one can derive temperature- 
dependent properties in an exact non-perturbative fash- 
ion. Applic ations of this approach include the Heisen- 
berg m odel (|Shiroishi and Takahashil . l2002h , higher-spin 
chains (iTsuboil l2003l 12004 and" integra bl e qua ntum 
spin ladders (|Batchelor et all 120031 l2004allbl . |2007| ). T- 
systems and integrability in general also play a funda- 
mental role in the gauge/string theories of high energy 
physics (jBeisert et all 120121 iKuniba et al] , l201ll ) . 

Yang-Baxter integrability has also played a crucial role 
in initiating and inspiring progress in mathematic s. Par- 
ticularly to the theory o f knots, links a nd braids dJones , 
19851 iKauffmanl. Il987t IWadati et all Il989b IWul. 11992: 



Yang and 2006[) and the devel opment of quantum 



groups and repres e ntatio n theory ( Chari and Presslevl 



group s ana representatic 
I1994I: iGomez et al.lll996l) 



2. Fermions in ID - historical overview 

In the mid-60 's many physicists worked on extend- 
ing the results obtained by Lieb and Liniger (1963) 
and McGuire (1964) for ID bosons with delta-function 
i nteraction to the pro blem of ID fermions. McGuire 
(IMcGuird . Il965l 119661) solved the eigenvalue problem of 
N — 1 fermions of the same spin and one fermion of 
opposite spin and studied the low-lying excited states 
with repulsive and attractive potentials. The dynamics 
of this on e spin -down Fermi problem has been studied 
(|McGuird . Il990h . The problem of AT - 2 fermions of the 
same spin with two fermions of opposite spin was solved 
by Flicker and Lieb (1967). A further step came when 
Gaudin (1967a, 1967b) and Yang (1967) solved the gen- 
eral problem in terms of a nested BA for arbitrary spin 
population imbalance. 1 Gaudin derived the groundstate 



1 Missing phase factors for the spin sector in Eq. (4c) of Gaudin 
(1967a) were corrected in Gaudin (1967b). A thorough treatment 
of the ID Fermi pro blem can be fo und in Gaudin's book on the 
Bcthc wavefunction llGaudinl.fT983) . 



energy for the balanced (fully paired) case for attractive 
interaction, pointing out t hat the result is equiva lent to 
that for repulsive bosons (jLieb and Linigerl . Il963h . The 
delta-function interacting two-component Fermi gas is 
commonly referred to as the Gaudin-Yang model. 

Yang's concise solution of the problem had a profound 
impact. As already remarked above, a key point in the so- 
lution is that the matrix operators describing many-body 
scattering can be factorised into two-body scattering ma- 
trices, provided that a set of cubic equations - the Yang- 
Baxter equation - are satisfied by the two-body scatter- 
ing matrices. This in turn is equivalent to no diffraction 
in the outgoing waves in three-body scattering processes. 
In this sense Yang's solution completes McGuire's for- 
mulation of the scattering process in the context of the 
ID Bose gas. Indeed, the i?-matrix obtained for the ID 
Bose gas is known as the s implest nontr ivial solution of 
the Yang-Baxter equation ( Jimbol . [l989[) . 

The solution of the ID Fermi problem triggered a se- 
ries of further breakthroughs. Yang (1968) obtained the 
(S-matrix of the delta-function in teracting many-bod y 
problem for Boltzmann statistics (|Gu and YaneL 11989b . 
The exact solution of the ID Fermi gas with higher 
spin s ymme try was obtained by Sutherland ( Sutherland! . 



spin s ymme 

Il968l . Il975t ). The ID Hubbard model solved by Lieb 
and Wu (1968) is a fundamental model in the theory 
of strongly correlated electron systems. Its solution 
is a significant example of the factorization condition 
(the YBE) in which the quasimomenta of particles k 
is replaced by sinfc. The Lieb-Wu solution thus gives 
a similar set of integral equations as Yang's Fredholm 
equations for the continuum gas. This exactly solved 
model has been extensively studied in the literature. 
The exact results for the Lieb-Wu model not only pro- 
vide the essenti al physics of ID st r ong l y correlated elec- 
tronic systems (jEssler et all 120051: [Hal . Il996l: ITakahashil . 
Il999h . but also are relevant to phenomena in high 
T c superconductivity. Indeed, the ID Hubbard model 
is an archetypical many-body system featuring Fuldc- 
Ferrel-Larkin-Ovchinnikov (FFLO) pairing, universal 
Tomonaga-Luttinger liquid (TLL) physics , spin-charge 
separation and quantum entanglement (lEssler et all 
120051 ; iGu et all 12004 lLarsson and Johannessonl 120051). 

Alt hough further study (|Takahashil . Il970bt lYand . 
1970) of the ID Fermi gas was initiated soon 
after its solution, it was not until much later 
that this model began to receive more atten- 
tion (lAstrakharchik et al. . 2004 ; Batchelor et al I I2006al 



iFuchs et qZj , 12004 llida and Wadatil . 120051: iTokatlvt I2004D 

as a result of the brilliant experimental progress in ultra- 
cold atom physics. The fundamental physics of the model 
is determined by the set of the generalised Fredholm in- 
tegral equations obtained in the thermodynamic limit. 
Takahashi (1970a) discussed the analyticity of the Fred- 
holm equations in the vanishing interaction limit. A thor- 
ough study of the Fredholm equations for the Gaudin- 
Yang model with att ractive and repuls i ve interactions 
has been carried out ( Guan et all l2007t iGuan and" Mai 



2012] ; llida and Wadatil 120071 120081; IWadati and lidal 
20071 ; IZhou et all 120121) . The numerical solution of the 



Fredholm equations has also been discussed in the con 
text of harmonic traps (|Colome-Tatchel [2008 ; Hu et al 



20071; iKakashvili and Bolechl 120091 iMa and Yangl . I200S , 
2010aibl : IOrsol . l2007h . In particular, the eigenfunction has 



been obtained explicitly for the Fermi gas in the infinitely 
strong repulsion lim it by using the hard-core contact 
boundary condit ion (iGirardeaul 19601) and group the o- 
retical methods (jGuan L. et aUl2009l:lMa et all l2010l) . 



The next major advance with implications for the ID 
Fermi problem was the solution of the finite temperature 
problem for ID bosons. Yang and Yang (1969) showed 
that the thermodynamics of the Lieb-Liniger Bose gas 
can be determined from the minimisation conditions of 
the Gibbs free energy subject to the BA equations. Taka- 
hashi went on to make significant contributions to Yang 
and Yang's grand canonical appro ach to the therm o- 
dynamics of ID i ntegrable models (ITakahashil Il971al fbl. 
11971 119731 11974 iTakahashi and Suzuki I1972D . 



Takahashi gave the general name Thermodynamic 
Bethe Ansatz (TBA) equations to the Yang- Yang type 
of equations for the thermodynamics. He discovered 
spin strings patterns to the BA equations in addi- 
tion to those for the groundstate of the spin chain 
(ITakahashl Il971al) . Using a similar spin string hy- 
pothesis, Gaudin (1971) studied the thermodynamics 
of the Heisenberg-Ising chain. Lai (1971a, 1973) in- 
dependently derived the TBA equations for spin-1/2 
fcrmions in the repulsive regime. It turns out that 
Takahashi's spin string hypothesis allows one to study 
the grand canonical ensemble for many ID many-body 
systems wi th internal degrees of freedom, e.g., the ID 
Fermi gas (ITakahashil Il971bl). the ID H ubbard model 



(ITakahashl 119721 11974 lUsuki e^ a/J.ll990l). t he quantum 



sine-Gord on model dFowler and Zotos L 1 98 it) . th e Kondo 



problem ( Filvov et all Il98lfc iLowensteinl Il98ll ) among 



many other integrable models. 

Building on Takahashi's spin string hypothesis, 
Schlottmann derived the TBA equations for SU(N) 
fermions with repu l sive and attractive interactions 
(ISchlottmannl . Il993l 11994 . The Yang- Yang method 



has been revealed to be an elegant way to analyti- 
cally access not only the thermodynamics, but also 
correlation functions, quantum criticality and TLL 



physics for a wide ra nge of low-dimensio nal c 
many-body syst ems ( Essler et all 120051 ; 151 
ITakahashi Il999h 
of the ID Bose 



u antum 
I, fl99rj 

The Yang- Yang thermodynamics 
s has been tested in the recent 



experiments ( 


Armiio et al.. 2010L 201ll Armiio. 2011 


Jacqmin et al. 


. 2011; 


Kriiger et al.. 20101 Sagi et al. 


2012 


; Stimming et al. 


2010t van Amerongen et al. 


2008 


)■ 



Recently, numerical schemes have been developed to 
solve the TBA equations of the ID two-comp onent spinor 
Bose gas with delta-function interaction ( Caux et all 
120091 ; iKlauser and Cauxl . l201lh . The QTM method has 
also been applied to the thermodynamics of the ID 
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Bose and Fermi gases with repulsive delta-function in- 
terac tion ( KKimper and Patul . 120111 ; iPatu and Kliimpeit 
l2012h . The Canberra g roup and their col l abora- 
tors (|Guan and Batcheloii 2011 ; Guan and Ho , 2011 



Guan et all \20l& iHe et all |2010L 120111; IZhao et al . 
20091 ) developed an asymptotic method to calculate 



the thermodynamics of strongly interacting bosons and 
fermions in an analytic fashion using the polylog func- 
tion in the framework of the Yang- Yang and Takahashi 
methods. This approach does away with the need to nu- 
merically solve the TBA equations for these systems at 
quantum criticality, where the temperature is very low 
and the inter-particle interaction is strong. 



B. Renewed interest in ID fermions 

The renewed interest over the past decade in ID 
fermions has been on a number of related fronts. Here we 
give a brief introductory outline of these developments. 



of the Fermi surfaces. The FFLO state has since 
been studied by various meth ods: density-matrix 
renormaliza t ion g r oup fDMRG) dLiischer et al. 
iRizzi et all 120081 : ITezuka and Uedal. l2010h. quantum 
Mont e Carlo (QMC ) (iBatrouni et all . 120081 ; iBaur et all 



120101 : IWolak et all |2010|). me a n fiel d theor y and 



other method s (IChen and Gaol. 2012t iDatta 



2009 



2009: 



2010: 



Devreese et all 1201 it lEdge and Cooperl. 
Kaiala et q/J. |201 it iLiu et all I2007T 2008bl ; iParish et al 
20071 ; iPei et all\201& IZhao and Liul . l2008l) . 



Very recently, the asymptotic correlation func- 
tions and FFLO signature were analytically stud- 
ied using the dressed charge fo rmalism in th e con- 
text of the Gaudin-Yang model ( Lee and Guanl . 1201 it 
ISchlottmann and Zvvaginl . l2012bh . However, a con- 
vincing theoretical proof for the existence of the ID 
FFLO state is still rather elusiv e , see recent further de- 
velopments dBolech et all l2012t jPalmonte et all l2012t 
iLu et q/.l . l2012t) . 



1. Novel BCS pairing states 

Quantum matter at low temperatures has already been 
seen to exhibit some remarkable physical properties, such 
as BEC and superfluidity. Fermionic quantum matter 
with mismatched Fermi surfaces has long been expected 
to exhibit more exotic behaviour than seen in conven- 
tional materials. The two-component attractive Fermi 
gas is particularly interesting due to its connection with 
the exotic pairing phase - the FFLO state - involving 
BCS pairs with nonzero centre-of-mass momenta. In this 
phase, where the system is partially polarized, the Fermi 
energies of spin-up and spin-down electrons become un- 
equal. Originally, Fulde and Ferrell (1964) discovered 
that under a strong external field, superconducting elec- 
tron pairs have nonzero pairing momentum and spin po- 
larization. Larkin and Ovchinnikov (1965) found that the 
formation of pairs of electrons with different momenta, 
i.e., k and —k+q with non zero <f is energetically favoured 
over pairs of electrons with opposite momenta, i.e., k and 
— k, when the separation between Fermi surfaces is suffi- 
ciently large. Consequently, the density of spins and the 
superconducting order parameter become periodic func- 
tions of the spatial coordinates. 

Theoretical study of the FFLO state in ID inter- 
acting fermions was initiated by K. Yang (2001), who 
used bosonization to study the pairing correlations. The 
FFLO-likc pair correlations and spin correlations for the 
attractive Hubbard mod el were later investigated nu- 
merically by two group s (iFeiguin and Heidrich-Meisnerl 
120071 ; ITezuka and Uedal . l2008l ) . Both groups showed the 
power-law decay of the form n pair oc cos(/cffloM)/M q 
for the pair correlation, with spatial oscillations de- 
pending solely on the mismatch /cfflo = it(n\ — 
n^) of the Fermi surfaces. Thus the momen- 
tum pair distribution has peaks at the mismatch 



2. Large-spin ultracold atomic fermions 



It was shown (|Ho and Yipl . Il999l ; lYip and Hoi . Il999l ) 



that large-spin atomic fermions exhibit rich pairing 
structures and collective modes in low-energy physics. 
Further progress towards understanding many- body 
physics with large-sp in Fermi gases was made (jWul . 
120051 : IWu et al. I. l2003f ) on spin-3/2 systems which c an be 
realiz ed with 132 Cs, 9 Be and 135 Ba ultracold atoms (|Wul . 
120061 ). Such systems exhibit a generic SO (5) symmetry 
(isomorphically, Sp(4) symmetry). The spin-3/2 system 
with SU(A) symmetry can exhibit a quartet state 
(four-body bound state). More generally, ultracold 
atoms offer an exciting opportunity to investigate 



large-spin symmetry (ICherng et al. 


. 20071; Corboz et al. 


2011 


; Honor kamp and Hofstettcr, 


2004 Krauser et al. 


2012 


; Rapp et al., 2007; Szirmai and Lewenstein 


2011 


; Tu et al. 20071: Zhail. 120071 Zhao et al. 2006 


Zhou and Semenofll, 20061). The trimer state ("trions") 



consisting of fermionic 6 Li atoms in the three energeti- 
cally lowest substates has been reported ( Huckans et all 
120091 ; lLonrpe" et al. 1, 120101 : IWilliams et aZ.l . l2009h . 

On the other hand, fermionic alkaline-earth atoms dis- 
play an exact SU(N) spin sym metry with N = 21 + 
1 where J is the nuclear spin (|Cazalilla et all 120091: 
iGorshkov et al. I, 120101 : Ixl I2010D . Such fermionic sys- 
tems with enlarged SU (N) spin symmetry are expected 
to display a remarkable diversity of new quantum phases 
and quantum critical phenomena due to the existence of 
multiple charge bound states. De Salvo et al. (2010) 
have reported quantum degeneracy in a gas of ultracold 
fermionic 87 Sr atoms with / = 9/2 in an optical dipole 
trap. An experiment by Taie et al. (2010) dramatically 
realised the model of fermionic atoms with SU(2)(£)SU (6) 
symmetry where electron spin decouples from its nuclear 
spin I = 5/2 for 173 Yb together with atoms of its spin- 
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1/2 isotope. This group also successfully realized the 
SU(6) Mott-insulator state with ul tracold fermions of 
173 Yb atoms in a 3D optical lattice (jTaie et aZ.1 . 12012ft . 

In the context of large-spin ultracold atomic fermions, 
Lecheminant et al. (2005) considered ID ultracold 
atomic systems of fermions with general half-integer 
spins. The instability of the BCS pairing phase and 
molecular superfluid phase in these systems have been 
studied by a low-energy approach. The low-energy 
physics and competing orders in large-spin fermionic 
systems in a ID lattice were f urther investigated 
(lAzaria et all 120091 ; ICapponi et~al\ . 120071: iNonne et all 
I2010L 201 lh . In this scenario, a new class of integrable 
models of ultracold fermi ons and bo s ons w i th large-spin 
symmetries were found ( Cao et al\ . 120071 : I Jiang etai\ , 
l2009l . 120111) . They derived the BA solutions for spin- 
3/2 fermions with SO(5) symmetry and the Sp(2s + 1)- 
invariant model of fermions. 

From the integrable model perspective, the study 
of multi-component attractive fermions was initiated a 
long time ago by Yang (1970) and Takahashi (1970b). 
In the light of ultracold higher spin atoms, Con- 
trozzi and Tsvclik (2006) proposed an exact solu- 
tion of a model describing the low energy physics of 
spin-3/2 fermionic atoms in a ID lattice. The ex- 
act results obtained from ID many-body systems with 
higher spin symmetries have provided insi ght into un- 
derstanding the few-body physics of trio ns ( Guan et all 
l2008al : iHe et all l2010t iLiu et all l2008aD; quartet states 
(four-body charge bound states) ( Guan et 'all , 120091 : 
ISchlottmann and Zvvagini l2012al blldl: and an arbi- 



trary large s pin-neutral bound st a te of different sizes 
dGuan et all l2010t iLee and Guanl. l201ll [Schlottmannl . 



I1993L 11994 lYang and Youl . l201lt lYin et ad . l2011aD . The 

study of critical phenomena and universal TLL physics in 
low-dimensional ultracold atomic Fermi gases with large 
pseudo-spin symmetries is a rapidly developing frontier 
in ultracold atom physics. 



3. Quantum criticality of ultracold atoms 

Quantum criticality describes a V-shaped phase of 
quantum critical matter fanning out to finite temper- 
atures from the quantum critical point (QCP). It is 
associated with competition between the two distinct 
ground states near the QCP. Near a QCP, the quan- 
tum critical behaviour is characterized by the energy gap 
A ~ £ -2 and a diverging length scale £ ~ \fi — fi c \~" , 
where /i c is the critical value of the driving parame- 
ter /i. The universality class of quantum criticality is 
characterized by the dyna mic critical exponent z and 
the c o rrelation expon ent v ([Fisher et all Il989t ISachdevl . 
119991 : IWilsonl . Il975l) . The many-body system is ex- 
pected to show universal scaling behaviour in the ther- 
modynamic quantities at quantum criticality due to the 
collective nature of many-body effects. Thus a uni- 
versal and scale-invariant description of the system is 



expected through the power-law scaling of thermody- 
namic properties. However, understanding the various 
aspects of quantum criticality in quantum systems rep- 
resents a major challenge to our knowledge of many-body 
physics (IColeman and Schofieldl 20051 iGegenwart et~al. . 



20081: lLohnevsen et aZ.1. 120071: ISachdev and KeimerLl2011 



Sondhi et all Il997l : IVoital . l2003h 

Ultracold atoms have become the tool of choice to 
simulate and test universal quantum critical phenom- 
ena. The study of quantum criticality and finite-size 
scaling in trapped atomic systems is thus attract- 



ing considerable interest (Cam 


postrini and Vicari. 


20091. 


2010allbl: Ceccarelli et all 


2012t Fang et al. 


2011a: 


Hazzard and Muellerl 20111: Kato et al. 


, 2008: 


Pollet et al., 20101). The experimental study of critical 



behaviour in a trapped Bose gas was initiated by Donncr 
et al. (2007). In particular, significant experimental 
progress on quantum criticality and quantum phase 
transitions in 2 D Bos e atomic gase s has b e en ma de 
dGemelke et all 120091: iHuang et all I2010L l2011al lbl: 
IZhang et aZ.I . l201lL l2012h . 

Some of the remarkable features of criticality in general 
are the notions of universality class and symmetry. Using 
integrable quantum field theory, Zamolodchikov (1989) 
was able to show that the 2D Ising model in a mag- 
netic field, or equi valently the quantum Isin g chain with 
a transverse field ( Henkel and Saleurl Il989f ) displays E s 
symmetry close to the critical point. Such exotic quan- 
tum symmetry in the excitation spectrum was observed 
in a recent experiment in the traditional s etting of con- 
densed matter physics (jColdea et al 1 I2010D . 

Zhou and Ho (2010) have proposed a precise theoret- 
ical scheme for mapping out quantum criticality of ul- 
tracold atoms. In this framework, exactly solvable mod- 
els of ultracold atoms, exhibiting quantum phase transi- 
tions, provide a rigorous way to explore quantum crit- 
icality in many-body systems. The equation of state 
has been obtained for a number of key integrable mod- 
els, allowing the exploration of TLL physics and quan- 
tum criticality. These include the Gaudin - Yang Fermi 
gas jGuan and Hoi . 1201 ll: lYin et all |2011at IZhao et al 



2009 ). the Lieb-Liniger Bose ga s dGuan and Batchelor , 
20111 ). the Fermi-Bose mixture (|Yin et all l2011bl) and 



the spin-1 spinor Bose gas with antiferromagn etic spin- 
spin exchange interaction ( Kuhn et all l2012al Fbh . The 
exact results for the scaling forms of thermodynamic 
properties in these systems near the critical point illus- 
trate the physical origin of quantum criticality in many- 
body systems. 



4. Experiments with ultracold atoms in ID 

Many remarkable ID quantum phenomena have been 
experimentally observed due to recent rapid progress 
in material synthesis and tunable manipulation of ul- 
tracold atoms. These developments have provided 
a better understanding of significant quantum sta- 
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tistical effects and strong correlation effects in low- 
dimensional quantum many-body systems. The ob- 
served results to date are seen to be in excellent 
agreement with results obtained using the mathemat- 
ical methods and analysis of exactly solved models. 
These include th e experimental realizati o n of the Tonks 



2006) 



Paredes et al 



Girar deau gas dKinoshita et all 12004 ; 
2004 ). a quantum Newton's cradle (IKinoshita et al 
and quantu m correlations 
2011 



Endres et al 



Guarrera et all 



dBetz et all 12011 
20121 lHaller et al 



20111 ; IKinoshita et all 120051; iTolra et all l2004h . In par 



ticular, Haller et al. (2009) made an experimental 
breakthrough by realising a stable highly excited gas- 
like phase, called the super Tonks-Girardeau gas, in the 
strongly attractive regime of bosonic Cesium atoms. 
The Yang- Yang thermodynamics and thermal fluctu- 



further tested in a series of publications ( 


Armiio et al. 


20101. 2011; 


Armiiol. 2011; Jacamin et al. 


. 2011L 2012 


Kriiger et al. 


, 2010; Sagi et all 20121 Stimming et al. 


2010l: van Amerongen et all 20081). The universal low- 



energy physics was demonstrated as hostin g a TLL 
(|Blumenstein et aZ.I . [20lTl lHaller etal\ . l2010al) . 

The experimental research using ultracold Fermi 
gases to explore pairing phenom ena in a ID Ferm i 
first reported in 2005 dMoritz et all 120051) . 



gas was tirst reported in 2005 (|Moritz et al\ 
In a major breakthrough towards understanding the 
exotic pairing signature and quantum phase diagram 
of the attractive Fermi gas, Liao et al. (2010) mea- 
sured the finite temperature density profiles of trapped 
fermionic 6 Li atoms, see FigfT] They confirmed the 
key features of the T = phase diagram predicted 



from the exact solution 


(Batchelor et all 2006a 


Feiguin and Heidrich-Meisnerl 20071: Guan et al. 


20071: Hu et al.. 20071 


Iida and Wadatil. 2007 


Kakashvili and Bolechl 120091: 


Orsol 2007; Parish et al. 


2007). 



C. Outline of this review 

In light of these recent developments we review the BA 
solution of the Gaudin-Yang model in Sec. II and discuss 
the physical understanding of the solution in terms of 
BCS pairing, the polaron problem, molecule states and 
the super Tonks-Girardeau gas. In Sec. Ill we further dis- 
cuss many-body phenomena in the Gaudin-Yang model. 
Especially, we discuss ID fermions in a harmonic trap and 
review the universal features of ID interaction, including 
magnetism, FLLO-likc pairing, TLL physics, spin-charge 
separation, universal thermodynamics and quantum crit- 
icality. In Sec. IV we review recent progress on mixtures 
of ultracold atoms and the exact solution of the ID Fcrmi- 
Bose mixture. 

Sec. V reviews the exotic many-body physics of 
ID multi-component interacting fermions, including the 
three-component Fermi gas, the SU (N) invariant Fermi 
gases, and spin-3/2 fermions with SO(5) symmetry. The 



discussion in this Section covers magnetism for sys- 
tems of large-spin fermions, trions, molecular states of 
different sizes, multi-component TLL phases, universal 
low-temperature thermodynamics and critical behaviour 
caused by population imbalance. In Sec. VI, we focus on 
the asymptotics of various relevant correlation functions 
for the Gaudin-Yang model and multi-component Fermi 
gases. The characteristics of the FFLO-likc pairing corre- 
lations and spin-charge separation correlation functions 
are discussed in the framework of conformal field theory. 

The experimental breakthroughs with quasi-lD ultra- 
cold atoms and tests of ID many-body physics are re- 
viewed in Sec. VII. A brief conclusion and an outlook on 
future developments are given in Sec. VIII. 



II. THE GAUDIN-YANG MODEL 



The Hamiltonian 



h 2 d 2 



-hi [ [<f\{x)^{x) - 4(a# 4 (*)) dx (1) 



describes a ID (5-function interacting two-component 
(spin-i) Fermi gas of N fermions with mass m and an 
external magnetic field H constrained by periodic bound- 
ary conditions to a line of length L. The function V(x) 
is the trapping potential. The field operators (f>i and (j>\ 
describe the fermionic atoms in the states \\.) and |f), 
respectively. The 5-typc interaction between fermions 
with opposite hyperfine states preserves the spin states 
such that the Zeeman term in the Hamiltonian (JXJ) is a 
conserved quantity. 

The expe ri menta l realization ( Liao et all |2010| ; 
iMoritz et all 120051 ) of this system of interacting 
fermions is described in Sec. VI. The coupling constant 
,9id = h 2 c/m, wher e c = — 2/aip c an be tuned by 
Feshb ach resonance ( Bergeman et all 120031 : lOlshanii 
119981) . For repulsive interaction c > and for attractive 
interaction c < 0. Where appropriate, we use units of 
h = 2m = 1. A dimensionless coupling constant 7 = c/n 
is used to characterize physical regimes, i.e., 7 ^> 1 for 
the strong coupling regime and 7 <C 1 for the weak 
coupling regime. Here n is the linear number density. 



A. Bethe ansatz solution 

For a homogeneous gas, i.e., V{x) = 0, the eigenvalue 
problem for Hamiltonian (TTJ) reduces to the ID A^-body 
delta-function interaction problem 



fi 2 N Pi 1 

9m ^ — ' Hnr. 



dx 2 



91D £ 

l<i<j<N 



5{x t - Xj ) (2) 
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solved by Gaudin and Yang. Bethe's hypothesis states 
that the wavefunction of such a many-body system is a 
superposition of plane waves. In the domain < xqi < 
XQ2 < ■ ■ ■ < xqn < L, the wave function is given by 



of fermions in the hyperfine levels |f) and \\.) such that 
Nf > JVj,. This gives the second eigenvalue problem 



Xi(ki)A E (P\Q) = exp(ifc l L)^(P|Q) 



(7) 







P 



[P,Q]expi(k P1 XQi + . . . + kp N XQN) (3) 



where both P = Pi,...,Pjv and Q = Qi, ■ ■ ■ ,Qn are 
permutations of the integers {1,2,..., N}. The sum runs 
over all AH permutations P and the coefficients [P, Q] of 
the exponentials can be arranged as an JV! x JV! matrix. 
A plane wave repeatedly reflected from the hypcrplancs 
XQ i = XQ j gives total of AM plane waves. The idea in set- 
ting up such an ansatz is an attempt at a hypothetical so- 
lution followed by demonstrating that it gives the eigen- 
function of the many-body problem rather than solving 
the problem directly. 

The derivative of the wavefunction is discontinuous 
when two particles arc infinitesimally close to one an- 
other. This property can be derived by considering the 
eigenvalue problem Hip = Eip in the center of mass co- 
ordinate X = (xQ i + xq .)/2 and the relative coordinate 
Y = (xQ i — XQ j ) of the two adjacent particles involved. 
This discontinuity in the first derivative of the wave func- 
tion and the continuity of the wave function at XQ i = xq 
give a two-body scattering relation between the adjacent 
coefficients = Yp_ p .£...jj,..., where £ p denotes an AM 

dimensional vector in the AM x AT! matrix. 



The matrix operator Y°j b is defined by 



-i(ki 



kj ) Pab 



(I 



\(ki — kj) — c 



(4) 



where / is the identity and P a b is the permutation opera- 
tor acting on the vector Due to the Fermi statis- 
tics, P ab = — 1 for all a and b. Yang denoted the unequal 
indices a, b, c, in the three particle scattering process as 
the interchanges with coordinates x a , x& and x c under 
the permutation Q a , Qb, Qc- The consistency condition 
for factorizing the many-body scattering matrix into the 
product of two-body scattering matrices Y°j b leads to the 
celebrated YBE 



yabybcyab 
jk ik ij 



\rbcyrabybc 
ij ik jk 



(5) 



Defining the P-matrix by Ri 



where KfKf = 1. 

PijY?j and spectral paramters u = &2 — k± and v = k^—k2 
the YBE is often written in the form 

Ri2(u)R2 3 (u + v)R 12 (v) = R,23(v)R 12 (u + v)R 23 (u). (6) 

Returning to solving the problem of A^ particles in a 
periodic box of length L, the second step is to apply 
the periodic boundary condition ip(xi, . . . , Xi, . . . , Xn) = 
ip{x\, . . . , Xi+L, . . . ,xpf) on the wavefunction with period 
L for every 1 < i < N. For an irreducible representation 
Rt/, = [2 N ±, l N *~ N ±], where Nf and N± arc the numbers 



where Ae(P\Q) is an abbreviation of the amplitude of 
the wave function which provides the eigenvector of the 
N operators with i = 1, . . . , N 

9t<(/Si) = R l+1 ,i(k l+1 - hi) . . .R N ,i{ k N - hi) 

xR lti (ki - fa) . . .Ri-\,i{ki-\ - h). (8) 

Using Bethe's hypothesis again, Yang solved the eigen- 
value problem (J7J) by the ansatz 

Ae{P\Q) = J2 a Pl ...p M F(X Pl , Vl )... F{\ PM ,y M ), (9) 



where yi < y 2 ■ ■ ■ < yu 
down-spin fermions and Ai , . 
within the function F(X,y) 



are the coordinates of the 
, Am are the spin rapidities 
TTV-i fctzA+ic^ By the 

llj — 1 fcj + i— A— lc' J 

symmetry P = [2^, l N *~ N i] } the vector Ae describes 
a spin system with a number of iVj, spins on an A^-site 
lattice. 

The generalized ansatz ^ plays an important 
role in soly i ng m ulti-component many-body problems 
( Sutherland! . Il968l ) . Alternatively, the eigenvalue prob- 
lem can be worked out in a straightforward way 
in terms of the QISM, where the operator 9^(fc;) 
can be written i n term s of the quantum transfer ma- 
trix (I Jiang et all l2009t iKorepin et all Il993l ; iLi et all 
120031 ; iMal. Il993t lOelkers et all 12006^) 7 This approach 
was int roduced in the s t udy of the ID H u bbard 
model (lEssler et all l2005t iMartins and Ramosl . Il998t 
iRamos and Martins! . Il997f j. 

The energy eigenspectrum is given in terms of the 
quasimomcnta {ki} of the fermions via 



E 



h 2 N 

— Yk* 



.fc? 

2m ^ 3 



(10) 



subject to the BA equations which in terms of the func- 
tion eb(x) 



_ x+ibc/2 
x-ibc/2 



are 



exp(ikiL) = J ei (fci - A a ) , 



a=l 



\ ei (Ac - kj) 

3=1 



J] e 2 (Xa - X P ) , (11) 

0=1 



for i = 1, 2, . . . , N and a = 1, 2, ... , N±. All wave num- 
bers kj are distinct and u niquely define the wave function 
© (jGu and Yand . ll989h . 

The fundamental physics of the model is determined 
by the BA equations (TTTj) . For repulsive interaction, the 
quasimomcnta {ki} arc real, but the rapidities {A a } are 
real only for the groundstate. The complex roots A a 
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are the spin strings for excited states. In the thermody- 
namic limit, i.e., L,N — >• oo, where N/L is finite, the BA 
equations (jllj) can be written as generalized Frcdholm 
equations 

ri(fc) = f 2 K x {k-k')r 2 {k')dk', 

2tt J_ B2 

Si 



r 2 (k) = I Ki{k - k')n(k')dk 
'l 

,B 2 



K 2 {k-k')r 2 {k')dk' 



-B 2 



(12) 



where the integration boundaries B\ and B 2 are deter- 
mined by n = N/L = n(k)dk, = NJL = 
J^b 2 r 2 {k')dk' . In the above equations, the kernel func- 
tion Ke(x) = 27 (tc/ 2 ) ll +x'> ■ ^ ne Iunc tions r m (k) denote 
the Bethe root distributions, with r\(k) the quasimo- 
menta distribution function and r 2 (k) the spin rapidity 
parameter distribution function. The groundstatc energy 
per unit length is given by E = jfg k 2 ri(k)dk. 

For the attractive regime, the quasimomenta {ki} 
of fermions with different spins form two-body bound 
states, i.e., the wave numbers are c omplex with kj 



\'i ± i c /2 in the thermodynamic limit ( Takahashl . 11971b! 
lYand. Il970ft . Here i = 1, 



N±. The excess fermions 



,N-2N 



I- 



have real quasimomenta {kj} with j = 1, 
In the thermodynamic limit, the density of unpaired 
fermions pi(k) and the density of pairs p 2 (k) satisfy the 
Frcdholm equations 



Pi (*) 



+ [ ' Ki(k-k')p 2 (k')dk' 

2n J-A 2 



1 



+ / K^k-k^pitk^dk' 
K 2 {k-k')p 2 {k')dk'. 

-A 2 



(13) 



The distribution p 2 (k) coincides with the distribution 
function of the real parts of the bound states. The lin- 
ear densities are defined by N/L = 2 j M p 2 {k)dk + 

J- Al Pi( k ) dk and N il L = I-A 2 P2(k)dk. The ground- 
state energy per unit length is given by E = 
$ M M (2k 2 - c 2 /2) p 2 (k)dk + f\ k 2 Pl {k)dk. 

In the context of magnetism, the magnetization per 
unit length is defined by M z = (n— 2n^)/2. By definition, 
the groundstate energy can be expressed as a function of 
total particle density n and magnetization M z . In the 
grand canonical ensemble, the magnetic field H and the 
chemical potential p can be obtained via the relations 



H 



dE{n,M z ) 
1 dM~ z : 



dE(n, M z 
dn 



(14) 



which have been used to work out the phase dia- 
grams of the attractive Fermi gas ( Hu et al\ l2007t lOrsol 



20071) and the repuls ive Fermi gas (jColome-Tatchel l2008t 
Guan and We now turn to extracting such in- 



formation from both the discrete and continuum versions 
of the BA solution. 



B. Solutions to the discrete Bethe ansatz equations 

The ID Fermi gas (flj with spin population imbalance 
exhibits an unconventional pairing order that presents a 
major subtlety of many-body correlations in the Gaudin- 
Yang model. Starting with the discrete BA equations 
([lip , we will review how the exact solution enables us 
to understand precisely such subtle many-body physics 
driven by the interaction of quantum statistics and dy- 
namics. In particular, we will see that for the weakly 
and strongly attractive coupling regimes ID interacting 
fermions give significantly different phenomena: weakly 
bound BCS-likc pairs vs tightly bound molecules. 



1. BCS-like pairing and tightly bound molecules 

For weakly attractive interaction, i.e., L\c\ <C 1, two 
fermions with spin-up and spin-down form a weakly 
bound pair with a small binding energy £b = —h 2 \c\/mL, 
where the t wo-body binding energy is less than the ki- 
netic energy ( Batchelor et all . l2006al ) . The complex con- 
jugate pair leads to an exponential decay of the wave 
function with a factor e~^ Xi ~ Xj ^ 2 . Thus the balanced 
case has a BCS-like fully paired state where the size of 
the Cooper pairs is much larger than the mean average 
distance between the fermions. In this weakly attrac- 
tive regime, the energy gap separating the first triplet 
excited state from the groundstate is found to have an 
asymptotic behaviour A 2n 2 \/-!rh\ exp (— 7r 2 /(2H)) 
as H — > ( Fuchs et al\ . 120041 iKrivnov and OvchinnikovL 
Il975l ). In fact, the BA equations (fTTj) for weak attrac- 
tion give an explicit relation H tn 4^j- (2it 2 m z + 4|7|m z ) 
between the external field and magnetization in the ther- 
modynamic limit. The lower critical field gives the energy 
gap at m z = 0. This relation indicates a vanishing en- 
ergy gap A = H c — > for 7 — > 0. Here the mag netization 
is de f ined by m z =: M z In = (Nf - N i )/2N dHe et ali . 
120091 : llida and Wadatl . l2007h . 

For a polarized gas with weak attraction, Fermi statis- 
tics leads to segmentation in quasimomentum space, i.e., 
the excess fermions are located at the two outer wings in 
quasimomentum space, see Fig. [21 For a finite size system 
with arbitrary polarization P — (iV-j- — Nj/)/N, the BA 
equations (TTT]) determine N± weakly bound BCS pairs 
with fcg w Aq ± \\/\c\/L and N — 2N .\. unpaired fermions 
with real fc, ( Batchelor et all l2006ct) . In this case, {A Q } 



and {kj} axe symmetrically distributed around zero in 
the quasimomentum parameter space, see Figj2] 

Assuming that N± is odd and N is even, the first few 
leading orders of the positive roots {A a } and {kj} are 
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Weak attraction 



-• •- 



-• •- 



Strong attraction 
--t« k 



-9 •- 



-* * •- 



# Spin up 

Spin down 

FIG. 2 Schematic BA root configurations for pairing and de- 
pairing in the Gaudin-Yang model. For weakly attractive in- 
teraction, the unpaired roots sit in the outer wings due to 
Fermi statistics. For strongly attractive interaction, the un- 
paired roots can pe netrate into the centra l region, occupied 
by th e bound pairs (|Batchelor et al . 2006c; I ida and Wadatil . 
120071 ). 



determined by the equations 



A a 



2n n ir 



c 

Lk 



3c 



c 
L 



L 



(-M±-l)/2 

£ 

a=l 
a 0=1 



k 2 M 



2A„ 



X 2 a — X 2 



(JV-2MJ/2 

— y 

2L ^ 

3=1 



2X a 



(15) 



where 
1,2,.. 



and 



M±±l A/+3 JV-A/4,-1 
2 ' 2 ' ' ' ' ' 2 ' 

2 . These case a = j3 is excluded in Eq. (fl5|). 
The root patterns reveal the cooperative nature of 
many-body effects, i.e., an individual quasimomentum 
depends on that of all the particles. Here the momenta 
of unpaired fermions and bound pairs depend on the scat- 
tering energies between pair and between paired and un- 
paired fermions. This indicates that the quantum statis- 
tics of the weakly interactin g fermions is m utual accord- 
ing to exclusion statistics (jHaldanej . 1 1 9 9 lh . From Eq. 
( TBI, the ground s tate en ergy per unit length is given by 
(jBatchelor et all l2006ch 



I = ^n 3 t n 2 + in^ 2 + 2cn tH + 0(c 2 ). (16) 



Here the groundstate energy (|16p is also valid for weakly 
repulsive interaction, i.e., for Lc< 1. This leading order 
correction to the interaction energy indicates a mean field 
effect. 

On the other hand, for strong attraction, i.e., L\c\ 3> 1 
(or c 3> kp), the discrete BA equations (fTTj) give the 
root patterns k\ w A; ± i^c for bound pairs and real kj 
for unpaired fermions (|Takahashil . Il971bt lYand . Il970t) . 
Here i = 1, . . . , iVj. and j = l,...,Ni, the number of 

unpaired fermions JVi = N — 2iVj,. The binding energy 

2 

£b = — % is the largest energy scale than the kinetic en- 
ergies of pairs or excess fermions. For a strong attraction 
(up to order 0(l/c 3 )), N fermions have root patterns 



(jBatchelor et all |2006<J| 
^ (2n« + l)7r 



x (2n b + 1)tt 
a u , A« — ah, 



L 2L 
where the effective statistical parameters are given by 



1 / 2N-2NA 1 



1 



4^ 

L\c\ 



and n u = -N x /2, -JVi/2 + 1, . . . , Ni/2 - 1 with n b = 
-NJ2, — JV4./2+I, . . . , 2V4./2-I. In this strong attraction 
limit, the groundstate energy per unit length is given by 
E/L = Eq + 2Eq + n.j,£b, where the energies of excess 
fermions and pairs are given by 



2 a 2 



E b = \n\-K 2 al. 



(17) 



The bound states behave like hard-core bosons which 
can be viewed as ideal particles with fractional exclusion 
statistics. However, the bound pairs have tails and they 
interfere with each other. It is impossible to separate 
the intermolecular forces from the interference between 
molecules and single fermions. From this explicit form of 
the groundstate energy, we see that for n^^> x = rif— nj,, 
the single atoms are repelled by the molecules, i.e., 



E(n i7 x) - E w ^ n l^ 2 



4:X 



12x(x + nj.) 



> 0. (18) 



Here Eg is the groundstate energy per unit length of 
the balanced gas. This result indicates that the single 
atoms are repelled by the molecules on the tightly bound 
dimer limit. The atom-dimer scattering problem of three 
fer mions in a quasi-one-dimen sional trap has been stud- 
ied ()Mora et all 12004 12005H . 



2. Highly polarized fermions: polaron vs molecule 



In higher 
immersed in 



gives rise to the 



dimensions, a spin-down fcrmion 
a fully polarized spin-up Fermi 
quasiparticle phen o menon 



called Fermi polaron dBruun and Massignanl 
Combescot and Giraudl 20081: ICombescot et all 



Parish, 



2011 



Klawunn and Recati|. | 2011t iMathv et al 



Prokof 'ev and Svistunovl.' 



2010: 



2007 



_2011 
l2008al lb 



Schmidt and Ensa . 12011 ). The Fermi polaron is a 



spin-down impurity fermion dressed by the surrounding 
scattered fermions in the spin-up Fermi sea. In particu- 
lar, recent observations of Fermi Polaron s in a 3D or 2D 
tunable Fermi liquid of ultracold atoms dKohstall et al 



20121: iKoschorreck et all 120121: iNascimbene et all I200S ; 



Schirotzek et all 120091) provide insightful understanding 
of quasiparticle physics in many-body systems. For 
an attractive polaron, with increasing attraction, the 
single spin-down fermion undergoes a polaron-molecule 
transition in the ferm i onic m edium ( Nascimbene et all . 
120091: ISchirol zek et oil. 120091) . 
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For repulsive interaction, theoretical studies have 
suggested the existence of such novel quasipar- 
ticles repulsive polaron s (iMassignan and Bruun 



20111: iNgampruetikorn et all [2012t iPilati et all I201C : 



Schmidt and Enssl . 120111; ISchmidt et all 120121) . The 



properties of repulsive polarons, such as the energy, life 
time, and quasiparticlc residue, give a fundamental un- 
derstanding of the coherent nature of the quasiparti- 
cle. The repulsive polaron is metastable and eventu- 
ally decays to either a molecule state or an attractive 
polaron with particle-hole excitations in the majority 
Fermi sea. This quasiparticle phenomenon was experi- 
mentally observed by a magnetically tuned Feshbach res- 
onance on the BEC side with positive scatteri ng length 
( Kohstall et aE l2012tlKosch.orreck et al.l . 120121 ). 

So far, most studies concerning the first or- 
der nature of the pola ron- molecule transition in a 
3D fermionic medium 
Combescot and Giraudl. 



Mathv et all l201ll) 



2010: 



2007 



Bruun and Massignan, 
2008t ICombescot et all 
involve a variational ansatz with 
some approximati ons that are ultimately not justified 
in low dimensions dGiraud and Combescovl . 120091 : iParishl 
1201 ll ) . It is generally accepted that there actually do not 
exist quasiparticlc excitations in ID systems due to the 
collective nature of the ID many-body effect. The el- 
ementary excitations in ID are still eigenstates, where 
all particles are involved in a low energy nature. There- 
fore, we cannot find a simple operator, acting on the 
groundstate, to get a quasiparticlc excitation, unlike for 
higher dimensional systems. But this does not rule out 
a well-defined quasi-particlc like behaviour, e.g., a po- 
laron, which is a typical example of the collective na- 
ture of the ID many-body effect. The quantum impurity 
problem in ID trapped ultracold atoms ha s shed new in- 
sight on the collective n ature of particles ( Catani et all 
l2012HPalzer et all l2009ft . 

The BA solvable models are likely to provide a rig- 
orous trea tment of pol a ron-like phenomena in different 



mediums (1 Guanl . 2012t iLeskinen et al. I, l20ld: iLi et all 



l2012t iMcGuirej . Il966l ). In particular, ID attractive po 



laronic phenomenon does occur if one (or a few) spin- 
down fermion (fer mions) is (are) immersed into a large 
spin- up Fermi sea dGiraud and Combescovl l2009t iGuaiil 



2012: Klawunn and Recati. 2011 : Leskinen et al., 201C : 



Massel et all l2012t IMcGuirej . Il966t IParishl 120111) . The 

excitation energy of a system with one spin-down fermion 
has a certain momentum-dependent relation, which in- 
cludes a mean field attractive binding energy plus a clas- 
sical kinetic energy of polaron with effective mass m* . A 
cartoon picture is shown in Fig. [3] 



McGuire (1965,1966) studied the exact eigenvalue 
problem of TV — 1 fermions of the same spin and one 
fermion of the opposite spin. He calculated the en- 
ergy shift caused by this extra spin-down fermion by 
solving the equation azi + 1/ tan Zj = constant with 
a = 4/(gL) for the quasi momentum ki = 2zt/L with 
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FIG. 3 Schematic BA root configuration of polaron-molecule 
crossover in the ID attractive Fermi gas. In the weakly at- 
tractive limit, the single impurity fermion (blue) dressed by 
the surrounding scattered spin-up fermions (red) from the 
medium behave like a polaron (dashed oval) with a mean field 
binding energy and an effective mass m* « m. For strong at- 
traction, it binds with one spin-up fermion from the Fermi 
sea to form a tightly bound molecule (yellow circle) of a two- 
atom with a sole molecule binding energy and an effective 
mass m* « 2m. 



i = 1, . . . , TV. Here g > for an attractive interaction 
strength. McGuire found a hermitian conjugate pair 
Zi,2 — ± i/3 and TV — 2 real roots Zi. The energy is 

2 v^W 



given by E = jji J2i=i ^f-This single impurity prob- 
lem was recently studied ( Guanl . |2012[ ) by means of the 
BA equations (|11[) . For an attractive interaction, the 
quasimomenta fc^-j- = p ± i/3 of a pair and TV — 2 real 
roots {ki} with i = 1, . . . , TV — 2 are deter mined by the 
equations (fTTj) with TVj. = 1. It was found ( Guanl . 120121 ) 



that the imaginary part /3 in the pair is determined by 
the equation j3L = tanh~ 



P\c\ 



ff2.fe 2 /4 " F° r an excited state 
with total momentum of the system q, the spin-down 
fermion associated with the weakly bound pair in the 
fully-polarized Fermi sea thus has a nonzero momentum 



/ 



l-2|c| 



E 



V 



(19) 



which depends on all individual momenta of the spin- 
up fermions. This givens a collective signature of the 
ID many-body effect. Thus the energy shift is explicitly 
given by 



£7(g > JV > JVi = l)-Bf(JV t ,0) 



2 n 2 



£p-b 



2m* 



(20) 



which behaves like a polaron quasiparticle. Here, the at- 
tractive mean field binding energy of the polaron is given 
by e p _h w — -t-cf|7| for weak attraction. The Fermi en- 
ergy is ep = ^\ n2 K 2 - We see that this binding energy 
depends solely on the Fermi energy of the medium and 
interaction strength in ID. In Eq. ([2T)|) . the polaron-like 
state has an effective mass m* ~ m(l + 0(c 2 )) which 
is almost the same as the actual mass of the fermions 
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due to the decoupling from the bound pair in the weak 
coupling limit. We point out that the polaron-like state 
only occurs for few impurity fermions immersed into a 
fully-polarized Fermi sea. 

For a weakly repulsive interaction and in the thermo- 
dynamic limit, the low energy physics of the ID Fermi 
gas is described by a spin-charge separation theory. The 
spin rapidity parameters decouple from the quaismo- 
menta of the fermions. However, using the BA equations 
(|11[) . a single spin-down fermion immersed into the ID 
fully-polarized Fermi sea with weak repulsion can form 
a repulsive Fermi polaron, with energy of the form (|20[) 
and an effective mass m* w m(l + 0(c 2 )). But here 
the impurity fermion receives a positive mean field shift 



>e F \l 



from the fermionic medium. 



In the opposite limit, a spin-down fermion immersed 
into a fully polarized spin-up medium with strong attrac- 
tion, i.e., with L\c\ 3> 1, the bound pair has = pH/3 
and the N — 2 real ro ots {ki} with i = 1, . . . , N — 2 are 
given by (iGuanl . l2012h 



ki 



Lit 
L 



Ap 



1 - 



L\c\ 



(21) 



with rij = ±1, ±3, . . . , ±(iV-f — 1). For an excited 
state with total system momentum q, the relation be- 
tween the centre-of-mass pair quasimomcntum p and 
the total momentum of the system q is given by p ~ 



which is independent of the indi- 
vidual quasimomenta of the spin-up fermions. The en- 
ergy shift is given by AE = Em — (J- with the chemical 
potential p = n 2 7r 2 , where the molecule energy is given 
by 

H 2 q 2 



E 



M 



2m* 



The binding energy of the bound state is 
E b 



2m 



T_ 
2 



37 



(22) 



(23) 



which tends to the binding energy of a sole molecule Sb = 
— j— — in the strongly attractive regime L\c\ — > oo. The 
effective mass of the molecule 



2m 1 



(24) 



becomes twice the actual mass of the fermions in this 
limit. From the shift energies (|20[) and (|22|) . we see clearly 
that as the attractive interaction grows, the spin-down 
fermion binds only with one spin-up fermion from the 
medium to gradually form a tightly bound molecule. The 
polaron-molecule crossover is regarded as a change from a 
mean field attractive binding energy of a polaron with an 
effective mass m* = m to the binding energy of a single 
molecule with an effective mass m* = 2m as the attrac- 
tion grows from zero to infinity. The non-equilibrium 
dynamics of an impurity in a ID lattice within a har- 
monic trap have bee n studied using num erical methods 



and the BA solution (Mas scl et al 



i g num e: 
• 2012). 



C. Solutions in the thermodynamic limit 

In the last section we discussed the solutions to the 
discrete BA (|TT|) in the limits |c| — ^ 0, oo. They give rise 
to different phenomena in the two extreme limits. Usu- 
ally, the many-body phenomena of interest refers to the 
physics of the system in the thermodynamic limit, where 
N, L — > oo keeping N/ L finite. This entails considering 
the solutions to the two sets of Fredholm equations (fT2")l 
and (|13p for the repulsive and attractive regimes. 



1. BCS-BEC crossover and fermionic super Tonks-Girardeau 
gas 

In order to conceive the physical nature of the super 
Tonks-Girardeau gas, we first recall the Lieb-Liniger Bose 
gas with zero-range delta-function interaction, where 
the Tonks-Girardeau gas was deter mined by a Fermi- 
Bosc mapping to an ideal Fermi gas ( Girardeau! . Il960f) . 
For a strong attractive interaction, McGuire (1964) pre- 
dicted that the quasimomenta of the bosons are given 
in terms of a bound state of A-particles, with k±j ' 



±\c [N - 2j + 1] with j = 1, . . . , N/2. In this case, the 
wave function is approximately given by 



(25) 



l<i<j<N 



where N = U/(n — l)!/\/27r) \c\ i ~ n ~ 1 ^ 2 is a normaliza- 
tion constant. The energy of the McGuire cluster state 
is given by E = —j^c 2 N(N 2 — 1). However, if the 
interaction strength is abruptly changed from strongly 
repulsive to strongly attractive, the highly excited gas- 
like state may be metastable against this cluster-like 
state due to the Fermi pressure inherited from the re- 
pulsive Tonks-Girardeau gas. This gas-like state exhibits 
a more exclusive quantum statistics than the free Fermi 
gas, and is called super Tonks-Girardeau gas. The super 
Tonks-Girardeau gas-like state was first predicted by As- 
trakharchik et al. (2005) and further proved by Batchelor 
et al. (2005b) using the exact BA solution of the Lieb- 
Liniger Bose gas. Remarkably, such a highly excited state 
was r ealized in a breakthrough experiment ( Haller et all 
12009ft . 

The equally populated components in an attractive 
Fermi gas give rise to physics related to the cross-over 
between a BCS superfluid and a BEC (IChen et all 201C ; 



Fcig uin et all 



20051 : iTokatlvl . 



2004lWadati and IidaU2007ft . In this con 

N/2, 



20121; iFuchs et all |2004 Iida and Wadati 



text, for a balanced attractive Fermi gas with 
the discrete BA equations (fTTj) reduce to 



cxp (2iX a L) 



n 

0=1 



A« — \p + icF 

A Q — Xr — icp , 



(26) 
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with cp = —2/af d . This is equivalent to the equations 



exp (ikjL) 



N B 

n 

1=1 



kj — ki + ic B 
k,j — ki — \c B 



(27) 



with cb = ^/ a \d f° r the Lieb-Liniger gas in the su- 
per Tonks-Girardcau phase under the identificatio n cb = 
2c F , N R = N/ 2 and m B = 2m F (jChen et all \201& 
IWadati and Tidal . l2007h . Since the bound pair formed by 
two fermions with opposite spin has a mass ms = 2m f, 
the M bound pairs are equivalently described by the su- 
per Tonks-Girardeau phase of the interacting Bose gas 
with effective ID scattering length af D = ^af D . This re- 
lation is also obtained by an exact mapping based on the 
two-body scattering problem associated with BEC-BCS 
crossover ( Mora et al I l2005al) . 

For the balanced Fermi gas, the binding energy is sub- 
tracted from the energy that gives the energy of the 
bosonic pairs of a two-atom, with result Eq = E+N^eb = 
2^ J2ati The energy eigenvalues of the bosons are 

given by E — J2f=i kj ■ ^ n ^is regard, the ground- 
state of the balanced attractive Fermi gas c an be viewed 
as th e fermionic super Tonks-Girardeau gas ( Chen et all 
2010). The identification between the balanced attrac- 
tive Fermi gas and the attractive Lieb-Liniger Bose gas 
suggests an effective attraction between pairs. 

In the thermodynamic limit, the BA equations (|26|) 
give a particular Fredholm equation which can be de- 
duced from (|T3"|) . i.e., 




i r=- 
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4.414 
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17.78 
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q/D 

FIG. 4 The normalized pair quasimomentum distribution 
function f(q) for differ ent values of interaction strength 7 
(|lida and Wadatl . l2005f ). The analytic result for the distri- 
bution function matches the numerical solution. The quasi- 
momentum distribution indicates the fermionization from mu- 
tual statistics to non-mutual generalized exclusion statistics 
as 7 increases. From Iida and Wadati (2005). 



nance has been repo rted ( Mora et~al\ , l2005allbl ). see also 
dFeiguin et a/J . 120121) . 



P2(k) = 



1 



K 2 (k - k')p 2 {k')dk' (28) 



where the Fermi pair momentum Q 2 is determined by 
n = 2 j^o, p 2 (k) . It t u rns out that (|Chen et dj l2010l 
llida and Wadatl . 120051 : I Wadati and Mai . 120071 ) the re- 
duced Fredholm Eq. (|28|) maps to the Lieb-Liniger in- 
tegral equation for ID spinless bosons on identifying 
rriB = 2n%F 1 Nb = Nf/2 and 7/3 = ^f- For a weak 
attraction, the groundstate is the BCS-like pairing state 
with a pairing correlation length larger than the av- 
erage intcrparticle spacing and the energy is given by 
E = jLn 3 7r 2 — i« 2 |c| + O(o 2 ). In particular, for strong 
attraction, the groundstate of bound pairs determined 
by (|2"5|) can be regarded as a parti cular super Tonks- 
Girardeau gas of hardcore bosons ( Chen etall , |2010() . 
The distribution function of the pair density plotted in 
Fig. H] provides an understanding of the subtle BEC-BCS 
crossover in the balanced Fermi gas. In the weak cou- 
pling regime, the single quasimomentum essentially de- 
pends on that of the oth er particles. Th i s gives a signa- 
ture of mutual statistics (lAneziris et all 1 199 It iHaldand . 
Il99ll ; IWilczekl . Il982t IWul 11994) . However, in the limit 
I7I — > 00, the quasimomentum distribution becomes an 
equally-spaced separation. This indicates free Fermi non- 
mutual statistics. Further study on the dimer-dimcr 
scattering properties in the confinement induced rcso- 



Furthermore, i t was demonstrated ( Girardeau! 120 lot 
iGuan and Chenl . l201dh that another metastable highly 
excited gas-like state without bound pairs in the strongly 
attractive regime can be realized through a sudden switch 
of the interaction from strongly repulsive to strongly 
attractive. In the limit c — > —00, this gas- like state 
is still an eigenstate of the system with the energy 



per particle E 



7T 2 (1 + 



4 In 2 



12 In 2 



), but it is 



3'" " 1 | 7 | 
a highly excited state. From the experimental point 
of view, these different quantum states can be tested 
from measuring the frequencies of the lowest breathing 
mode from the mean square radius of the ID trapped 
gases in a harmonic po t ential ( Astrakharchik et q/.l . l2004l ; 
iMenotti and Stringaril 20021). e.g.. in the super Tonks- 
Girardeau Bose gas ( Haller et all 20091 ). The low breath- 
ing mode featuring different states of the strongly repul- 
sive and attractive Fermi gas can be analyzed via the lo- 
cal density approximation. The lowest breathing mode is 
given by the mean square radius of the trapped fermionic 
Tonks-Girardeau gas uj 2 = — 2{x 2 ) / (d{x 2 ) / duj 2 .) , sec 
Fig. [5] Here (x 2 ) = J p(x)x 2 dx/N . The frequency 
ratio ui 2 /uj 2 exhibits a peak which is a typical charac- 
teristic of the super Tonks-Girardeau phase. Further 
evidence for this fermionic super Tonks-Girardeau gas- 
like state has been seen in the experimental observa- 
tion of the fermion ization of two distinguishable fermions 
(IZiirn et q/.l . l2012D . 
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FIG. 5 The square of the lowest breathing mode frequency 
vs the groundstate energy per unit length vs the rescaled in- 
teraction strength aia/w^. The quantum gases are trapped 
in a ID harmonic potential V x = ^muj 2 x 2 . Here GS and 
FSTG stand for the frequency ratio ui 2 /uj 2 for the ground- 
state and fermionic super Tonks-Girardeau gas, respectively. 
From L.-M. Guan and Chen (2010). 



2. Solutions to the Fredholm equations and analyticity 



hyE/L 



E% 



EH 



El 



2E b 



(n t - n;) 3 7r 2 



where 



8n± 



15|d3 



(l27T 2 (n T - nj 2 



4807V 2 ; + 5?iJtt 2 ) 



(30) 



(3 



1 



2(2nt-n 4 ) 3(2n t -n+) 5 



-— i -J3 (l80n ; n 2 + 207r 2 n 3 . - 90n t n| - 22 7 r 2 n 3 . 
15|c| 

+15nJ - 120^ + 637r 2 ^ri t - 607r 2 n ; n^)] . (31) 



This energy is highly accurate for arbitrary polarization 
as can be seen in Fig. [6j The high precision of expansions 
for the groundst ate energy of the attractive Fermi gas 
were also studied (|He et all 120091 ; E da and Wadatil . [2007t 
IZhou et aZ.l . l2012t) . 

In contrast to the strong coupling case, it is more dif- 
ficult to proceed with asymptotic expansion for the two 
sets of Fredholm equations ([T2"j) and (fT3"|) at vanishing in- 
teraction strength. In terms of the polarization P, the 
gro undstate energy in weak attraction limit was found to 
be (|lida and Wadatil . l2007f ) 



Despite the two sets of Fredholm integral equations 
(|T2|) and ([TBI for the homogeneous gas bein g derived 
long ago (jTakahashl [l971ti lYand . Il967l Il970lh their an- 
alytical study are still restricted to particular regimes, 
e.g., 7 > 1, I7I < 1 and < -1- Tnc nrst f ew 
terms in the asymptotic expansions of the groundstate 
energy for the ID attractive Fermi gas for both the 
strong and weak cou pling cases has been calculated in 
terms of power s e ries (IGuan et al. l. l2007t lHe" et al. I. I2009t 
llida and Wadatil. I2005L 120071) and in terms of Legcndrc 



polynomials ( Zhou et ali l2012j ). The first few terms 



of the groundstate e nergy have been derived recently 
( Guan and Mal . [20T2T ) by an asymptotic expansion for (a) 
strong repulsion, (b) weak repulsion, (c) weak attraction 
and (d) strong attraction. 

For strong repulsion, the gro undstate energy of t he 
Gaudin-Yang model is given by ()Guan and MaL l2012ft 



frrV 2 

32(ln2) 



1 



4 In 2 1 

7 

2 C(3) 



12(ln2) 



7 

in 3 7r 2 



7° 

1 _ . 



4S/1T 



£ (256n 3 - f 7r 2 n 2 n; 



for P = 0, 



for P > 0.5. 



(29) 



Here £(z) is the Riemann zeta function. The leading or- 
der (I/7) correction was also found in (jBatchelor et all 
l2006d;lFuchs et aUl200l . Fig. ©shows that this ground- 
state energy is a good approximation for the balanced 
and imbalanced Fermi gas with a strongly repulsive in- 
teraction. 

For strong attraction, the groundstate energy is given 



E 
L 



12 



(1 + 3P 2 )-4 F (1-P 2 )| 7 |-P27 2 



(32) 

The coefficient B2 is a complicated function obtained 
from the power series expansions with respect to 7. How- 
ever, it contains di vergent sums and the c oefficients are 
singular as 7 — > (jlida and Wadatil . |2007[) . So far, only 
the leading order correction to the interaction energy is 
mathematically convincing and consistent with the result 
(fTrJ)) obtained from the discrete BA equations ([TT]) . Be- 
yond the mean field term, finding the next leading term 
in the groundstate energy is still an open problem. For 
zero polarization, Krivnov and Ovchinnikov (1975) found 

0(c 2 ) M -^(ln| 7 |) 2 obtained from the ID Hubbard 
model in dilute limit. Iida and Wadati (2007) found the 
term 0(c 2 ) ps — 7 2 n 3 /12. This difference reveals a sub- 
tlety of the vanishing interaction limit, i.e., the two limits 
c — > and the thermodynamic limit (N, L — > 00 with 
N/L finite) do not commute. In fact, the groundstate 
energy (|32[) only counts the density distributions away 
from the integration boundaries in the Fredholm equa- 
tions ([HI) and ([13]), i.e., \B t - k\ ~ c and |A; - k\ ~ c 
with i = 1,2. At the integration edges, these distribu- 
tion functions are no t analytically extractable as 7 — > 
(IGuan and Mai . I2012D . 

The analyticit y of the groundstate energy at 7 = i s 
still in question jGuan and Mai . 120121 rTakahashilll970cl) . 
Takahashi showed that (a) the groundstate energy func- 
tion f(n^,ny,c) is analytic on the real c axis when 
n-f =/= and (b) f(n^, nj,; c) is analytic on the real c axis 
except for c = when — n±. However, the Fredhom 
equations for weakly repulsive and attractive interactions 
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FIG. 6 The groundstate energy as a function of 7 = cL/N in 

units of 2m ^ 2 . The figure shows the comparison between the 
asymptotic solutions and numerical solutions of the Fredholm 
equations for different polarization. In the attractive regime, 
the binding energy et = — c 2 /2 was subtracted. From Guan 
and Ma (2012). 
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FIG. 7 Phase diagram for finite strong interaction with | — y | = 
10 and density n — 1. The dashed lines denote the two critical 
lines (|34[) and {55) . The coloured phases are obtained by 
numerical solution of the energy-magnetization ()33[) . From 
He et al. (2009). 



are identical as long as B\ > B2 and A\ > A2 where the 
integr ation boundari e s mat ch each other between the two 
sides (|Guan and Mai . l2012h . In this identical region, the 
asymptotic expansions of the energies of the repulsive and 
attractive fermions are identical to all orders as c — > 0. 
But the identity of the asymptotic expansions may not 
mean that the energy analytically connects due to the di- 
vergence of the Fredholm equations in the region c — > iO. 



III. MANY-BODY PHYSICS OF THE GAUDIN-YANG 
MODEL 

So far, we have only discussed the groundstate prop- 
erties of the Gaudin-Yang model. We now survey the 
wide range of fundamental many-body physics exhibited 
in the model. 



A. ID analog of the FFLO state and magnetism 

The particularly interesting feature of the attrac- 
tive Fermi gas is the exotic FFLO-like pairing, where 
the system is gaplcss with mismatched Fermi points 
between the two Fermi seas. In the gapped phase, 
it is well understood that the correlation function 
for the si ngle particle Green's function decays expo- 
nentially (|Bogoliubov and Korepml Il988l Il989l Il99dh 
(ipl gV'i.s) — > e~ x ^ with £ = vp/A and s =t, |, whereas 
the singlet pair correlation function decays as a power 

x~ e . Here A is 



of distance, i.e., (ipt^fpt^ip 1,^1,1) 
the energy gap, and the critical exponents £ and 9 are 
both greater than zero. However, once the external field 
exceeds the lower critical field, the system has a gaplcss 



phase where both of these correlation functions decay as a 
power of distance and the pairs lose their dominance. The 
molecule and excess fermions form the polarized FFLO- 
like pairing state, where the spacial oscillations of the 
pairing correlation are caused by an imbalance in the 
densities of spin-up and spin-down fermions, i.e., i%f — nj,. 
In section I VII we will further discuss the pair and spin 
correlations with the spacial oscillation signature in the 
context of conformal field theory. 

In terms of the polarization, the Gaudin-Yang model 
with attractive interaction exhibits three quantum phases 
at zero temperature: the fully paired phase which is a 
quasicondcnsate with zero polarization, the fully polar- 
ized normal Fermi gas with P = 1, and the partially 
polarized FFLO-like phase with polarization < P < 1, 
see Fig. [7] In this gappless phase, the magnetic prop- 
erties can be exactly described by the external field- 
magnetization relation 
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-e b + fi u -[i° 



(33) 



where [i b = \i + eb/2 and ^t" = /1 + H/2 are given by 
(|14|) . This relation reveals an important energy transfer 
relation among the binding energy, the variation of Fermi 
surfaces and the external field. 

For fixed density and strong attraction, the paired 
phase with magnetization M z = is stable when the 
field H < H c i , where the lower critical field is given by 
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When the external field exceeds the upper critical field 
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Magnetic Field H/e b 

FIG. 8 Magnetization vs the external field H/eb f° r c = — 10 
in the units 2m = H = 1 for different densities n. The dashed 
lines are plotted from the analytic result (|36[) . The solid 
curves are obtained from numerical solutions of the dressed 
energy equations. The inset shows similar comparison be- 
tween analytic and numerical results for the susceptibility vs 
external field H/e b . From He et al. (2009). 



showed that when the magnetic field H — >• H c , the 
magnetization m z ~ \J H — H c for a fixed chemical po- 
tential. However, for fixed density, the magnetization 
m z ~ (H - H c )/{-KV b N ) as H = H c + 0+. This leads 
to a finite onset susceptibility given by \ — l/( 7ru Ar) 
with the pair density stiffness v h N — i>f/4 in the strong 
attraction limit 7 — > 00. Here we further remark 
that for finitely strong attraction the onset susccptibil- 
ity X = # (6) /Mv) wh crc a (1 + ^ + ^) is 

the Tomonaga-Luttinger liquid parameter at the criti- 
cal point and v b N = ~f + f^i + 4^ J i s the stiffness of 
bound pairs in the limit H — >• i? c + + . 

The magnetization in the Hubbard model with a half- 
filled ba nd gives rise to t he square root dependence on 
the field (jTakahashil Il969l ) , where low density solitons ap- 
pearing in the spi n sector above the critical field behave 
like free fermions (iJaparidze and Nersesvanl . Il978l Il98lt 
iPokrovskv and Talapov . 1979T ). Rigorously speaking, the 
linear-field dependent magnetization is clearly seen from 
the energy transfer relation ([3"3"| . where the effective 
chemical potentials /j, u oc (2m z ) 2 and fi b oc (n — 2m z ) 2 . 
Thus the linear term m z in the relation ([55)) gives the 
finite susceptibility at the onset of magnetization. 



a phase transition from the FFLO-like phase into the 
normal gas phase occurs, see Fig. [7] The lower critical 
field gives the energy gap in the spin sector. 

The magnetization can be ob t ained from Q33I). see 
Fig. H It was fo u nd dGuan et all [2007t iHe et all 12009 



llida and 1 Wadati l2007t IWovnarovich and Pend . 119911 ) 
that in the vicinity of the critical fields H c i and H C 2, the 
system exhibits a linear field-dependent magnetization 



M z 



2 

b 



11 

2 7 2 



1 



H, 



i-H 



1 + H 



1 81-Tr" 

r 6h\ 3 

12 _ 16(^-6) 
T 1 3b7F 



(36) 



with a finite susceptibility. For fixed total number of 
particles, or say in a canonical ensemble, the magnetic 
field driven phase transitions in the ID Fermi gases 
with an attractive interaction are linear-field-dependcnt, 
which was also found in the SU (N) attractive Fermi gas 
dGuan et all l2010l : iLee et aLl . l20ll . 

The magnetism of the attr active Fermi g a ses has been 
discu ssed by Schlottmann (ISchlottmannl . 1 19931 1 1994 
119971 : ISchlott maim and Zvvaginil2012allbl ). However, the 
argument which was made o n the initial slope of the 
magn etization in these papers ( Schlottmannl . I"l993. 1994, 
Il997l ) does not appear to be correct for a fi xed total num- 
ber o f particles. The reason is discussed in ( Wovnarovichl 
Il99ll ): "the bound pairs which have to be broken up 
to yield the particles with uncompensated spins form a 
Fermi sea, their density of states is finite at the Fermi 
level, and that keeps the in i tial s usceptibility finite". It 
is also shown ( Vekua et all l2009h that the curvature of 
free dispersion at the Fermi points couples the spin and 
change modes and leads to a linear critical behaviour and 
finite susceptibility for a wide range of models. They 



B. Fermions in a ID harmonic trap 

In experiments, ID quantum atomic gases are pre- 
pared by loading ultracold atoms in an anisotropic 
harmonic trap with strong transverse confinement and 
weak longitudinal confinement. In general, interact- 
ing many-body systems trapped in an harmonic po- 
tential is a rather complicated problem. The prob- 
lem of the ID Hamiltonian ([T]) trapped in an harmonic 
pote ntial ^muxjx 2 has been studied by various meth- 
ods (IColome-Tatchel 120081: ICul [2012bt iGao and Asgari . 



20081; iGirardeaul. l2010t iGirardeau and Mingguzzil 12007 



Hu et all l2007nMa and YaneL 120091: lOrsol . l2007t 1 1 Yang 
2009l ; lYin erafl . hoilbl) . For N = 2, the eigenvalue prob- 



lem of the trapped gas has been studied analytic ally in 
dBusch et all Il998t lldziaszek and Calarcol 120061). The 



energy shift for N = 2 ( Busch et all Il998l) is given by 



r(-_E/2+i/4) — l/ a iD> where an> is a scattering length 



/77 IY-.E/2+3/4) _ 

v^ r (-j5/2+i/4) - 

and T(x) is the Euler gamma function (IBusch et all 
Il998h . The system of two for mions with arbitrary interac- 
tion in a ID harmonic potentia l has been experimentally 
investigated ( Ziirn et all , l2012f ). In this experiment, the 
Tonks-Girardeau state and the metastable super Tonks- 
Girardeau state have been observed. 

This problem for arbitrary number of particles was 
studied analytica lly ( Guan L. et all [20091 : iMa and Yangl . 
120091 : lYand . 120091) .' where the limiting cases c — > ±00 and 
c = have been studied using group theory. In partic- 
ular, Yang (2009) gave an analysis of the groundstate 
energy of fermions in a ID trap with 5-function inter- 
action. In the light of Yang's argument, for any value 
of interaction strength, the eigenvalue problems: (a) the 
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trapped Hamiltonian with symmetry Y = [N — M, M] 
in full oo^ space and (b) the Hamiltonian in region Ry 
with the boundary condition that the wave function van- 
ishes on its surface, are equivalent. Here the region Ry 
is bounded by C^_ M x Cj^ planes at which the wave 
function vanishes. For any value of g, the ground- 
state wave function for problem (b) has no zeros in the 
interior of Ry and is not degenerate. This thus suggests 
that the groundstate energy of the system with total spin 
J = N/2 — M increases monotonically and approaches to 
the energy E J=N / 2 - The Lieb and Mattis theorem fur- 
ther suggests Ej > Eji if J > J'. 

For c — > oo, the groundstate energy of the trapped gas 
with total spin J is given by Ej = J2n=o (5 + n ) = 
\N 2 , which is independent of the total spin J. For 
c = and J = N/2 — M, the energy is given by 
Ej = \ ([N/2 + J] 2 + [N/2 - J} 2 ) . Ma and Yang (2009) 

argued that Ej/N 2 



f'{M with 



1/2 fort -> 00 

fj(t) = { l/4+(J/A0 2 fort = 

- (1/2 - J/N)t 2 /A fort -> -00 



l2007f ). This type of approximate treatment is known as 
the local density approximation (LDA). 

To ensure the validity of the LDA, the correlation 
length £(z) should be much smaller than the character- 
istic inhomogeneity length £ inh = ld ^ z y dzl , i.e., £(2) < 
£inh- The two length scales £(z) and ^i n h are determined 
by the local chemical potential /i(z) and the local density 
n(z). Obviously, from the definition of £j n ft, the LDA be- 
comes invalid near the edge of an atomic cloud where the 
density drops rapidly. However, in real measurements, al- 
most all signal strengths are proportional to the density. 
Therefore, due to the very small density at the edge, the 
central region of large density dominates the measure- 
ment signals. For a large number of particles, N — > 00, 
the density profiles of the trapped gas can be precisely 
analyzed within the LDA. 

In an harmonic trap, the equation of state (|14p can 
be reformulated within the LDA by the replacement 
/i (x) = /.i (0) — ^muj 2 x 2 in which x is the position and 
uj x is the frequency within the trap. Using the LDA for 
the ID Bose gas in an harmonic trap, its g lobal chemical 
potential reads ( Kheruntsvan et a/1 120051) 



where t = -7=. In particular, for c — > 00, the exact wave 
function of the system VE' = tpAipJ where the spatial wave 
function and sym metric spin wave fun ction have been 
derived explicitly ( Guan L. et all [20091) 



det [<f>j(xi)] J i= i',...,N 

[N—M,M]i 



(N\) 

N\/((N-M)\M\) 

4kj = £ {yl N ~ M ' M] Q a }z a 



Hg = fMi[n(z)] - V(z) = (jlo[ti(z)] - -mcjfz 2 , (38) 

where the local chemical potential /io[^(z)] at position 
z is given by the chemical potential for a homogeneous 
system of an uniform density n = n(z). The total number 
of atoms N is given as N = J n(z)dz. 

Similarly, for a ID two-component Fermi gas 
in an harmonic trap, the global chemical poten- 
tial is dHeidrich-Meisner et all l2010at iHu et all l2007t 
iMa and Yand . |2010at lOrsol 120071) 



for the symmetry R 
with Qi = nLi Uf=M 

,[M,M] 



[N - M, M]. Here Q a = P a Q x 
[+1 sgn(x; — Xj). The basis tensor 

functio n yg VI ' JVJ> was constru cted explicitly from group 
theory (jGuan L. et aLl . l2009h . 



The fcrmion density distribution for ID interacting 



tions on top of a uniform density cloud (IGao et ai. 


20061: Gao and Asgaril. 20081: 


Guan L. et all 200£ 


Ma and Yane. 20091: Rieol et al. 


, 2003J). These os- 



^hom[n(x),P(x)} 



/i - ^muj 2 .x 2 



dilations can be described by an analytical form as ( | ao an sg. 
of the density distri bution dButts and RokhsarL Il997t 
iGleisberg et all l2000t ISoffing et adl20lil )~ 



where the chemical potential /J,hom[n, P] can be obtained 
from the homogenous gas (fT4]) . n(x) is the total lin- 
ear number density and P(x) is the local spin polar- 
ization. They can be determined from restriction on 
the total particle number N = n (x) dx and po- 
lari zation P = f^° r ^n u (x )dx/N which are rewritten 
(IGao and AsgariU 2008t IHu et all l2007t lOrsol . l2007t 



n(x) ss tlq{x) 



(—1) 2 cos{2k,F{x)x) 
ttLf 1 — x 2 /L 2 F 



(37) 



Nai D /4 



(Naj D ) P 



h (x) dx 

n u (x)dxxa 2 x . (39) 



for x < Lp, where the density cloud is given by the 



Thomas- Fermi profile, i.e., no(x) 



2ujL f 



^l-x 2 /L 2 



with a Thomas- Fermi radius Lp = J N/cu. If the longi- 
tudinal confinement is weak enough, the atomic density 
varies smoothly along the longitudinal direction and so 
the ato mic gases can be treated as locally homogeneous 
systems ( Hu et al\ . 120071 : kheruntsvan et al. l. l2005t [O rso. 



Here h (x) = l/\'y(x)\, a x = ^/h/(rnui x ) and n" (x) = 
n u (x)/\c\. 

Calculations for the integrable homogenous attractive 
gas confined to a ID trapping potential thus lead to a 
two-shell structure composed of a partially polarized ID 
FFLO-likc state in the trapping centre surrounded by 
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FIG. 9 Quantum phase diagram of the Gaudin-Yang model in 
the fi — H plane. The phase boundaries are obtained from Eq. 
(|14|l in terms of the numerical solution of the BA equations 
dHJ. From Orso (2007). 



wings composed of eit her a fully paired state or a fully 
polar i zed Fermi gas ( Gao and Asgaril . 120081 : iHu et all . 
120071 ; lOrsol 120071) . see Fig. M This prediction was ver- 
ified by Liao et al. (2010) by the observation of three 
distinct phases in experimental measurements of ultra- 
cold 6 Li atoms in an array of ID tubes. The analytical 
study of the phase diag ram of the ID attractive Fermi gas 
has been presented in (iGuan et al. 1. 120071 : lG uan and 
l201ll ; llida and Wadatil . |2007l ). 



C. Tomonaga-Luttinger liquids 



The TLL (|Luttingerl . [l96l iTomonagal Il950t ). describ- 
ing the collective motion of bosons, has played an im- 
portant role in the novel description of universal low 
energy phy s ics fo r low-dimension al many-body physics 
(jGiamarchil . 12004 iGogolinl Il998[) . In ID systems of in- 
teracting bosons, fermions or spin systems, the effect 
of quantum fluctuations is strong enough to yield strik- 
ing anomalous quantum phenomena. In this approach, 
for example, the low energy physics of a ID interacting 
fcrmion system can be described by a bilinear form of 
bosonic creation and annihilation operators. The TLL is 
henomenologically treated by bosonization techniques 



phenomcnologically trc 

(ICazalilla et all .120111: iGiamarchil . 12004 IGogolinl Il998t 
iTsvelik and Wiegmann . 19831) based on a linearization 



of the dispersion relation of the particles in the collective 
motion, i.e., uj(q) = v s \q\, here v s is the sound velocity 
of the collective motion. In contrast to the Fermi liq- 
uid, this thus leads to a power law density of states for 
the TLL at the Fermi energy Ep, i.e., \E — Ep\ a , where 
the exponent a = (K + 1/K — 2) /4 depending on the 
so-called TLL parameter K. 

In general, the correlation functions of such ID systems 
at zero temperature show a power-law decay determined 
by the TLL parameter K and the velocity v s . These crit- 



ical systems not only have global scale invariance but ex- 
hibit local conformal invariance. With the help of exact 
BA solutions, a wide class of ID interacting systems can 
be mapped onto TLLs in the low-energy limit, including 



spin-charge separati 


on (lEssler et all, 20051: Giamarchi 


2004; Kawakami and Yand. 19901: SchulzL |1990L 1991 


Solvoml. 


19791; IVoitl. 


1995). Moreover, progress in treat- 



ing such collective motion of particles beyond the low- 
energy limit was made by Imambekov and Glazman 
(2009a,2009b) and Imambekov et al. ( 2011). This 
method can be applied to a wide variety of ID systems 
with collective motion of particles. This generalized TLL 
theory could be possibly justified through exact BA re- 
sults for ID intcgrable models in ultracold atoms and 
correlated electronic systems. 

In contrast to the conventional quasiparticles carry- 
ing both spin and charge degrees, the elementary exci- 
tations form spin and cha rge waves that propagate with 
different velocities in ID ( Gogolinl . Il998f ). The relativis- 
tic dispersion relation for each one of these excitations is 

written as w^(p) = 



where A, 



i M is the energy 

gap and is the velocity. For a gapless excitation with 
vanishing energy gap i> M = d p 0J tl (p). For ID interacting 
systems, this gives a phonon dispersion that leads to con- 
formal invariance in the excitation spectrum. However, 
for a large energy gap, the dispersion can be rewritten as 
w M (p) = A M + U 2p 2 /(2A AI ) := A M + p 2 / (2m* ) which is 
the classical dispersion of a free particle with an effective 
mass m*. 

From the BA solution (fTTj) with TVf = N±, the charge 
and spin velocitie s are v r ._s = \vf (l ± 7/7 T 2 ) for the weak 
coup ling regime (jBatchelor et all . l2006allcl ; iFuchs et all . 
120041 ). Here the Fermi velocity vf = Hnn/m. For 
strong attraction, the charge and spin velocities are 
given by v c = \ v f(^ ~ I/7) and v„ = y / A(l — 2/7) 
(|Batchelor et a/J . l2006al lc1: IFuchs et all 120041 ) with an en- 



ergy gap A sa 7^^-. This gap increases with increasing 
interaction strength 7 so that the spin velocity is diver- 
gent in the strongly attractive limit. However, for strong 
repulsion the charge velocity v c = vf{1 — 4 In 2/7) tends 
to the Fermi velocity and the spin velocity goes to zero 
v s = 3 ^(1 - 6 In 2/7) (|Lee et aZ.I . [2012h due to suppres- 
sion of spin transportation due to the strong repulsion, 
see Fig. [TU] We shall discuss TLLs and spin-charge sep- 
aration phenomena in the attractive Fermi gas in the 
following two sections. 



D. Universal thermodynamics and Tomonaga-Luttinger 
liquids in attractive fermions 

The Yang- Yang formalism with its generalisation for 
the study o f therm odynamics of BA integrable systems 
( Takahashil . Il999h is a convenient tool for the study 
of universal thermodynamics and quantum criticality in 
the presence of external fields. At finite temperatures 
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FIG. 10 Charge and spin velocities vs dimensionless inter- 
action strength 7 for the ID balanced Fermi gas. The solid 
lines are the velocities obtained by numerically solving the 
BA equations (|11[) . The dotted lines denote the analytical 
result for the velocity in the weak coupling regimes. From 
Batchelor et at. (2006a). 



and in the thermodynamic limit, the densities in the 
Fredholm equations (|13p evolve into occupied and un- 
occupied roots in the whole parameter spaces, namely 
A\,A% — > 00. In particular, the roots in the spin sec- 
tor form complicated string patterns that characterize 
the spin excitations, i.e., spin wave bound states. The 
density distribution functions of pairs, unpaired fermions 
and spin strings involve the densities of 'particles' pi{k) 
and 'holes' Pi(k) (i = 1,2). Following the Yang- Yang 
grand canonical ensemble method, the grand partition 
function is written as Z = tr(e~ w / 7 ") = e~ G / T , in terms 
of the Gibbs free energy G = E - HM Z - pn - TS 
and the magnetic field H, the che mical potential p 
and the entropy S (|Takahashil Il999t) . In terms of the 
dressed energies e b (fc) := T \n(p2(k) / p 2 {k)) and e u (fc) :— 
T\n(pi(k)/pi(k)) for paired and unpaired fermions, the 
equilibrium states are determined by the minimization 
condition of the Gibbs free energy, which gives rise to 
a set of co upled nonlinear i ntegral equations - the TBA 
equations dTakahashil . ll999D . For the attractive Gaudin- 
Yang model, these equations are 



e\k) 
6 u (fc) 



2(fc 2 -p- X -<?) + TK 2 * ln(l + c- eb «/ T ) 



■TK 1 * ln(l + e" 



\k)/T 



2 - p - -H + TK X * ln(l + c 

OO 

-T^^*ln(l + r;7 1 (fc)), 



-e h (k)/T 



(40) 



\n m (X) = ^ + X, * m(l + e~ eU ( A >/ T ) 



gr^lnCl + T?" 1 ^)). (41) 



The function r\n (A) := ^(A)/^(A) is the ratio of the 
string densities. Here * denotes the convolution integral 
(/*<?)(A) = ,f.l/ (A-AW , T he function T lm {k) 
is given, e.g., in (|Guan et all 120071: iTakahashil . Il999t) . 
The Gibbs free energy per unit length is given by G — 
p b + p u where the effective pressures of the bound pairs 
and unpaired fermions are given by 



P b = 



P 



T 



dfcln(l + e- e W/ T ), 

1 

— / dfcln(l+e- £U W/ T ). (42) 
2tt J_ qo 



In the grand canonical ensemble, the total number of 
particles associated with the chemical potential p can 
be changed. The Fermi sea of unpaired fermions can be 
lifted by the external field. The entropy S is a mea- 
sure of the thermal disorder. The spin fluctuations (spin 
strings) are ferromagnetically coupled to the Fermi sea 
of unpaired fermions. The direct numerical computa- 
tion of the TBA equations was presented by Kakashvili 
and Bolech (2009). The TBA equations for this model 
involve an infinite number of coupled nonlinear integral 
equations that impose a number of challenges to access- 
ing the physics of the model. 

For zero external field, the lowest excitations split into 
collective excitations carrying charge and collective exci- 
tations carrying spin. This leads to the phenomenon of 
spin-charge separation. The charge excitations are de- 
scribed by sound modes with a linear dispersion. How- 
ever, for the external field in excesses of the lower criti- 
cal field the spin gap vanishes. In contrast to the spin- 
charge separation formalism, the spin-charge coupling 
drastically changes the critical behaviour in the attrac- 
tive regime of the Fermi gas. The TBA equations (|4Tj) 
indicate that the spin fluctuations (the spin wave bound 
states) are ferromagn etically coupled t o the Fermi sea 
of unpaired fermions ( Zhao et all . I2009T I . In contrast to 
the antiferromagnetic coupling Jaf = — j|[P u (^\ H) for 
repulsive regime (jGuan et HB. l2008bh . the spin-spin ex- 



change interaction in the spin sector is described by an 
effective spin- 1/2 ferromagnetic chain with a coupling 
constant Jp « ||y p u (T, H) > in the strong coupling 
regime |c| 3> 1. The ferromagnetic spin wave fluctuations 
are produced due to the thermal fluctuation in the Fermi 
sea of unpaired fermions. However, Jp tends to zero 
for 7 — > 00. Therefore the spin transportation becomes 
weaker and weaker until it vanishes as I7I — > 00. At zero 
temperature all unpaired fermions are polarized and spin 
strings arc fully suppressed. In this gaplcss phase, exci- 
tations involve particle-hole excitations and spin-string 
excitations. The TBA equations (|4Tj) can be greatly sim- 
plified in the strong coupling regime due to the suppres- 
sion of spin fluctuations, where 77^" 1 ~ e -£H/T _s. q as 
T — > 0. Thus one can extract the universal TLL physics 
using Sommcrfeld expansion for temperatures less than 
chemical potential and magnetic field. 

In fact, in this spinless phase, the spin fluctuation is 
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suppressed in the limit T — s- and | — ^ | ^> 1. Thus 
the bound pairs and unpaired fermions form a two- 
component TLL. Conformal invariance predicts that the 
energy per unit length has a universal finite-size scaling 
form that is characterized by the dimensionless number 
C, whic h is the centra l charge of the und erlying Virasoro 
algebra (|Affleckl . Il986l iBlote et all Il986l) . The finite-size 
corrections to the groundstate ener gy have been analyt- 
ically derived ( Lee and Guar! [20 111) 



£n - 



Ctt 
6L* 



a — u , b 



(43) 



where C = 1 with v u and Vb the velocities of unpaired 
fermions and bound pairs, respectively. For strong inter- 
action, they are given explicitly by 
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(44) 



where A\ = 4ri2, A 2 = 2ni+n 2 and n 2 = nj,. We will de- 
scribe universal behaviour of the macroscopic properties 
of this Fermi gas in the following subsections. 

Although a phase transition in ID many-body sys- 
tems at finite temperatures does not exist, the system 
does exhibit universal crossover from relativistic disper- 
sions to quadratic dispersions. Thus at low tempera- 
tures, the bound pairs, normal Fermi gas, and the FFLO 
phase become relativistic TLLs of bound pairs (TLLp), 
unpaired fermions (TLLp), and a two-component TLL 
(TLLpp), respective ly, see Fig. [TT1 A detailed discus- 
sion has been given (|Yi et all l2012t I Zhao et~cd\ , 120091) . 
For the temperature ksT <C Ep, where Ep is the Fermi 
energy, the leading low temperature correction to the 
free energy of the polarized gas can be calculated ex- 
plicitly using Sommerfeld exp a nsion with the pressures 
(1421). namely (iBatchelor et all \201& iGuan et all l2007t 
iHe et all 1200a IZhao et all l2009h 



F(T,H) 




irCk 



6h 

irCk B T 

nCk 2 B T 2 i 
6h VF ' 



for TLLpp 
tor TLLp 
tor TLL F 

(45) 

which belongs to the universality class of the Gaussian 
model with central charge C = 1. For strong attraction, 
the velocities are given in (|4"4")l . In the above equation, 
the groundstate energy E (H) is as given in Sec. IH.B.ll 
In fact, from the TBA equations (|4"Tj) . the universal ther- 
modynamics (|45[) can be shown to be valid for arbitrary 
interaction strength. 

The two branches of gapless excitations in the ID 
FFLO-likc phase form collective motions of particles. 
The low energy (long wavelength) physics is described 




FIG. 11 Quantum phase diagram of the Gaudin-Yang model 
in the T — H plane showing a contour plot of the entropy 
in the strong interaction regime. The dashed lines are de- 
termined from the deviation from linear temperature depen- 
dent entropy obtained from the result (|45|) . The universal 
crossover temperatures separate the TLLs from quantum crit- 
ical regimes. From Yi et al. (2012). 



by an effective Hamiltonian 



H cS = ■£ [(d x u ) 2 + (d x 9 u ) 2 ] + -f [(dzfo) 2 + (d x 



h d x (j) u _ (d x <j) u + 2d x (/) b ) 

2 ^ y/H 



(46) 



as long as the spin fluct uation is frozen out (jVekua et all 
l2009l ; lzhao et aZj . l2009h . Here the fields d x <f>i,d x 6i with 
i = b,u are the density and current fluctuations for 
the pairs and unpaired fermions. However, in the spin 
gapped phase, i.e., for H < H c i, the energy gap in 
the spin sector leads to an exponential decay of spin 
correlations, whereas the singlet pair correlation and 
charge density wave correlations have a power-law decay 
(|Cazalilla et all 120051 ; iGao et all l2007f ). Thus the sys- 
tem in the spi n gapped phase forms a so-called Luther- 
Emery liquid ( Luther and Emervl Il974l ) . 



E. Quantum criticality and universal scaling 

As we have seen, the ID attractive Fermi gas exhibits 
various phases of strongly correlated quantum liquids and 
is thus particularly valuable to investigate quantum crit- 
icality. Near a quantum critical point, the many-body 
system is expected to show universal scaling behaviour in 
the thermodynamic quantities due to the collective na- 
ture of the many-body effects. In the framework of Yang- 
Yang TBA thermodynamics, exactly solvable models of 
ultracold atoms, exhibiting quantum phase transitions, 
provide a rigorous way to treat quantum criticality in 
archetypical quantum many-body systems, such as the 
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Gaudin-Yang Fer mi gas dGuan and Hoi. 1201 ll) . the Lieb- 
Liniger Bose gas dGuan and Batchelorl l201l[) . a mixture 
of bosons and fermions ( Yin et ai . 20121 ). and the spin-1 
Bose gas with both delta- function interaction and anti- 
ferromagnetic interaction ( Kuhn et all l2012al lbl) . 

At zero temperature, the quantum phase diagram 
in the grand canonical ensemble can be analytically 
deter mined from the so-c a lled d ressed energy equa- 
tions dGuan and Batchelorl 1201 ll iGuan and Hoi 1201 it 
iTakahashil . Il999f ) 



'(A) 



A 2 



M 7 



K 2 (A- A')e b (A')dA' 



K ± (A- k)e u (k)dk, 



-A! 



-A 2 



K x (k~ A)e b (A)dA, (47) 



which are obtained from the TBA equations in the limit 
T — > 0. The integration boundaries ^2 and charac- 
terize the Fermi surfaces for bound pairs and unpaired 
fermions, respectively. It is convenient to use dimension- 
less quantities where energy and length are measured in 
units of binding energy Eb and c _1 respectively. In terms 
of the dimensionless quantities fl := yu/et, h := H/eb, 
t := T/sb, n ■= n/\c\ = 7 -1 , p := P/\ceb\, the phase 
bou ndaries have been de termined analytically from (|47|) , 
see (jGuan and Hoi . l201ll) . There are four phases denoted 
by vacuum (V), fully paired phase (P), ferromagnetic 
phase (F), and partially paired (PP) or (FFLO-\ikc) 
phase presenting the same phase diagram as in Fig. [9j 

The low density and strong coupling limits are partic- 
ularly important to study quantum criticality. In fact, 
the TBA equations (|4"Tj) can be converted into a dimen- 
sionless form wi th the above reseati ng. Following the 
notation used in ( Guan and Hoi . [20 111 ), the phase bound- 
aries between V — F, V - P, F - PP and P - PP are 
denoted by [i c \ to {i c ± respectively. The closed forms of 
the critical fields 



Mel 



Mc3 



Mc2 



= -2 + 3^ (1 -^ + 2^ (1 - /l)2 - (48) 
are needed to determine scaling functions of thermody- 
namic properties. Here [i c \^c2 applies to all regimes 
and fj, C 3,fJ-c4 are expressions in the strongly interacting 
regime. The above critical fields /i C 3,/^ C 4 correspond to 
the upper and lower critical fields in the h — n plane , 
which were foun d in dGuan et all 120071 ; iHe et all l2009t 
llida and Wadatil . l2007l) .~ 

The TBA equations (gO]) and (gl]) encode the micro- 
scopic roles of each single particle that lead to a global 
coherent state - quantum criticality. The quantum crit- 
icality is manifested by universal scaling of thermody- 
namic properties near the critical points. The key input 



to obtain critical scaling behaviour is to derive the form of 
the equation of state which takes full thermal and quan- 
tum fluctuations at low temperatur es into account. The 
dimensionless form of the pressure ()Guan and Hol . l201lh 



p(t,fi,h) := p/(|c|e 6 ) =p 



(49) 



serves as equation of state, where the pressures of the 
bound pairs and unpaired fermions are given by 
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In these equations the functions F% 



1 n ' 



ft and /« 

arc defined by F^ u := Li„ (- e x ".-/*) and fy u := 
Li„ (— e Vb ' u ' *) , with the notation Vb — 2ji + 1, v u = 
jl + h/2. The function Li s (z) = Y2kLi z /k s is the poly- 
log function and Iq{x) = J2T=o (Mp^f ) 2k ■ Despite the 
equation of state (j4"9")l having only a few leading terms 
in expansions with respect to the interaction strength, 
it contains thermal fluctuations in contrast to the TLL 
thermodynamics (|45l) . 

The TLL thermodynamics (|4"5")) has been derived from 
low temperature expansion along T <C \fx — fx c \. This 
universal thermodynamics is a consequence of th e lin- 
early dispersing phonon modes ( Maeda et aZ.I . [2007h . i.e., 
the long wavelength density fluctuations of two weakly 
coupled gases or a gas of bound pairs or single fermions. 
The quantum critical regime lies beyond T » \fx — fx c \. 
In this limit, the equation of state (HHJ) provides closed 
forms for the scaling functions of thermodynamic quanti- 
ties, such as density, magnetization and the compressibil- 
ity. Near the critical point, the thermodyna mic functions 
can be cast into a u niversal scaling form ( Fisher et all 
ll989HSachdevlll999h . The explicit universal scaling form 
of the density for T >• \fi — /U c | is 
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(51) 



Here the constants n 3 and n 4 are the background den- 
sities near the critical points /13 and ^4. These constants, 
together with a and b , are known explicitly in terms of h 
(IGuan and Hol . l201lh . 

The universal scaling form for the compressibility is 
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(F-PP) K 
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wher e k 3 , k D 4 , A4 , A5 are also known ( Guan and Hoi . 

Houl). 

In the Gaudin-Yang model, the above density and com- 
pressibility can be cast into the universal scaling forms 

n(M,T,x) = no + ^ +1 -^( M ^ Mc ), (53) 

k{^T,x) = Ko + T^-Zf ^W-^ y (54) 

with dimensionality d = 1. Here the scaling functions 
are Q{x) = A Q Li!/ 2 (— c x ) and J-(x) = A ( gLi„ 1 / 2 (— e x ) 
from which one can read off the dynamical critical ex- 
ponent 2 = 2 and correlation length exponent v = 1/2 
for different phases of the spin states. Such results illus- 
trate the microscopic origin of the quantum criticality of 
different spin states, i.e., the singular parts in (|53[) and 
(|54"]) characterize sudden changes of the density of state 
of either excess fermions or bound pairs. The TLL is 
maintained below the cross-over temperature T* which 
indicates a universal crossover from a r elativistic dis- 
persi on into a nonrelativistic dispersion ( Maeda et all 
120071 ). The quantum criticality driven by the external 
field H gives rise to the same universality class. 

Using the LDA presented in ()39p . quantum criticality 
of the bulk system can be mapped out through finite tem- 
perature density profiles in the trapped gas. For small 
polarization, the chemical potential passes the lower crit- 
ical point /i C 2 = — \ from the vacuum into the fully paired 
phase then passes the upper critical point /i C 4 from the 
fully paired phase into the FFLO-like phase. At finite 
temperatures, quantum criticality of the Fermi gas can be 
seen clearly from contour plots of entropy in T — /i plane 
see Fig. IT27 a). The typical U-shape crossover temper- 
ature T* separates the quantum critical regimes where 
T* oc — fi c \. The crossover temperatures are deter- 
mined by minimums or maximums of the magnetization 
or by the breakdown of the linear-temperature-dependent 
entropy (jZhao et aZ.I , l2009lb Fig.[n](b) and (c) show that 
the unpaired density curves for different temperatures in- 
tersect at the critical points /x C 2 and /z C 4, respectively. 

The phase boundary separating the fully paired phase 
from the FFLO-like phase can be mapped out from the 
density profiles of unpaired fermions in the trapped gas 
at finite temperatures. As the temperature decreases, 
the compressibility evolves a round peak sitting in the 
phase of the higher density of state. It diverges at zero 
temperature. Similarly, for the high polarization case, 
the density profiles of unpaired and paired atoms can 
be used to map out the phase boun daries /x r i ( V — > F 
) and fi c3 (F PP), respectively (|Yin et all koiibl) . 
This signature can be used to confirm the quantum crit- 
i cal law, as per th e recent experimental measurements 
(|Zhang et all l2012t) . 




FIG. 12 Quantum criticality of the ID Fermi gas in an har- 
monic trap for low polarization: (a) contour plot of the en- 
tropy in the t — /i plane. Symbols indicate the crossover 
temperature T* separating the quantum critical regimes from 
vacuum, single component TLLp of paired fermions and two- 
component TLLpp of FLLO-like states. The intersection of 
the density curves at different temperatures can map out the 
critical points (b) /j, C 2 and (c) ^ C 4. The corresponding com- 
pressibility curves (d) and (e) intersect at the same critical 
points after a subtraction of the background compressibility. 
From Yin et al. (2011). 



evide n t in th e sc aling analysis dCampostrini and Vicari . 
20091. l2010al lbl: ICeccarelli et al\ . l2012t IZhou and Ho . 
20101) . It has been proved that quantum criticality of 
the bulk system can be revealed from the singular part 
of a thermodynamic quantity near the trapping centre 
x = 0. The scaling behaviour exists in the limit of large 
trapping size. E.g., under a scale change b, the singular 
part of the density below the c ritical dimension d c can 
be written (jZhou and Hoi . l2010h 



The universal scaling behaviour of the homogenous 
system can be mapped out through the density profiles 
of the trapped gas at finite temperatures. However, in- 
homogencity caused by the finite-size scaling effect is 



n((i,T,u; 2 ,x) = b-^+^^g (pb L / v ,Tb*,w !i V',x/b\ 

with p = fi — fi c and y = 2 + l/v. Choosing 
Tb z = 1, the scaling function with finite-size trapping 
is Q(p,, T\D, x) = G (j2/T 1 ^" z \D,xT 1 ^ z ) with D = 
w 2 /T y / z . The scaling behaviour is revealed through plot- 
ting g against fj, at x = for different temperatures 
with a fixed D. This means that the scaling behaviour 
of the homogeneous system could be extracted from th e 
trapping centre in a small window ( Zhou and Hoi . l2010h . 
Nevertheless, the finite- size error lies with in the current 
experimental accuracy ( Zhang et aZ.I . [20T2I ) . One can ei- 
ther lower the temperature or increase the interaction 
strength such that all data curves for the physical proper- 
ties at different temperatures collapse into a single curve 
with a proper scaling in the trapped gas. 
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F. Spin-charge separation in repulsive fermions 

Landau's Fermi liquid theory provides a universal de- 
scription of low energy physics of interacting electron sys- 
tems in higher dimensions, where interactions only lead 
to finit e renorm alizations of physical properties ( Landaul 
Il957allblll958h . The quasiparticle excitations have a di- 
vergent lifetime when the excitation energy goes to zero. 
Renormalization of individual quasiparticles leads to a 
similar Fermi liquid behaviour, e.g., a finite density of 
states and a step-like singularity in momentum distribu- 
tion at zero temperature. The deviations from the values 
of physical properties of noninteracting systems present 
the interaction effect. However, the low energy physics 
of ID interacting many-body systems does not have such 
quasiparticle-type excitations. In ID many-body sys- 
tems, all particles participate in the low energy excita- 
tions and form collective motions of the charge and spin 
densities with different velocities (ICazalilla et all 2011 



Giamarchl [2004 iGoeolinl . 119981: [Tsvelik and Wiegmann 
19831) . Thus the low energy physics only depends on 



the TLL parameter and the velocities of collective charge 
and spin oscillations. Introducing charge and spin Boson 
fields C , CT = ± /V2, n c , a = (n t ± ILJ /V2, the 
low energy physics of the ID spin- 1/2 repulsi ve Fermi gas 
can be described by an effective Hamiltonian (Jg iamarchil . 
12004 ISchulzl . Il99lh 



H = H C + H a 



2<?i 



(27ra) 2 



dx cos( 



(55) 



The fields <p v and n„ obey the standard Bose commu- 
nication relations 11^] = \5 VIjL 5{x — y) with u, (j, = c, u. 
The parameter a is a short-distance cutoff. The last 
term in the effective Hamiltonian (|55[) characterizes the 
backscattcring process, i.e, it corresponds to a 2kp scat- 
tering. The ID interacting Fermi gas separates into 
charge and spin parts 



H,, = I dx 



■kv v K v ^ 2 



2itK, 



(56) 



The coefficients for different proce sses are given phe- 

chl, 12004) . The coefficient 



nomenologically (jGiamaref 
v c /K c is the energy cost for changing the particle density 
while v„ / K a determines the energy for creating a nonzero 
spin polarization. The compressibility and susceptibility 
are given by k = 2K c /{ttv c ) and x = 2A' (T /(7rv cr ), respec- 
tively. 

At fixed point g\ = the specific heat is given by 
a linear temperature dependent relation c = 7 C T where 
IcHo = { v f/v c + Vp/v a ) /2. Here 70 is the specific heat 
coefficient of noninteracting fermions. The susceptibility 
is given by x/xo = Vf/v<t where the noninteracting sus- 
ceptibility reads xo = ^/{^vp). Thus the Wilson ratio at 
the fixed point is given by 



Rw 



2v r 



V c + Vcr 



(57) 



This ratio gives a universal feature of collective mo- 
tions of particles. The TLL parameter K v and charge 
and spin velocities v v determine universal behaviour 
of correlation functions (ICheianov and Zvonare 3, 12004 
iGiamarchiL 12004 ISchuld . Il991ri . 

On the other hand, in terms of the TBA formalism, the 
physical quantities are described by the set of nonline ar 
integral equations ([Lai Il97l[ Il973f iTakahashil Il971bf ) 



e{k) = k 2 -n-*L-TY, A,*ln(l + e-^W/ T ) 
bj (A) = jH - TKj ;* In (l + e - £(A)/T ) 

OO 

+T^r jm *ln(l+c-^W/ T 



m—1 



(58) 



with j = 1, . . . ,00. Here Tj m {k) is given in Takahashi 
(1999). The free energy per unit length F and the pres- 
sure P follow as 

rp />00 

Fttfin-P, P=tt ln (! + c~ £{k)/T )dk, (59) 

where n denotes the particle density. The spin wave 
bound states give an antiferromagnetic ordering at low 
temperatures. The population imbalance associated with 
the Fredholm equations (fT2j) lead to three distinguishable 
phases: the spin singlet groundstatc with magnetization 
m z = (zero magnetic field); a magnetic phase with fi- 
nite magnetization (0 < H < H s ) and a fully polarized 
phase with m z = 1/2 (H > H s ). The critical field, at 
which the density of down-spin atoms is zero, is given by 
dLee et al\ . \2Qm 
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2^ 



27m 2 I tan 1 



27m 



(60) 



The phase diagram in the chemical potential-magnetic 
field plane is shown in Fig. [T51 

Spin charge separation is a hallmark of the TLL physics 
for ID interacting fermions. For arbitrary repulsive cou- 
pling c > in arbitra ry magnet i c field H < H Cl the TBA 
equations (f5"8"]) yield ( Lee et all . I2012T ) 



F = E n 



1 1 
\ _ + _ 

6 \ V s V c 



(61) 



where v c and v s are, respectively, the holon and spinon 
excitation velocities 



e' c (k ) 
2irp c (k ) ' 
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(62) 



However, the velocities v c and v s can be analytically 
calculated only for strong and weak interactions. The 
numerical solutions of spin a nd charge veloci t ies hav e 
been discussed in the literature (IBatchelor et aZl. l2006allcl : 
iFuchs et all 12004 iRecati et all |2003[ ). see Fig. [TO] In 
an harmonic trap, there is an imbalanccd mixture of 
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FIG. 13 The zero temperature phase diagram of the Gaudin- 
Yang model in the repulsive regime. Here the chemical po- 
tential and magnetic field are rescaled by €>, = . PP 
denotes partially polarized phase with finite magnetization. 
FP denotes the fully polarized phase. 
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FIG. 14 A short laser pulse near the trap centre could ex- 
cite the charge density and spin density wave packets in an 
harmonic trapped Fermi gas. Charge and spin velocities can 
be manifested in a spatial separation of the spin (solid curve) 
and density (dashed curve) wave packets. From Recati et al. 
(2003). 



two-component fermions in the trapping center and fully 
polarized fermions at the edge (lAbedinpoure£ al. , l2007t 
IColome-Tatchl 120081: IrJa and Yand . I2009L l2010al) . The 
exact analytical solution of quasi-one-dimensional spin- 
1/2 fermions with infinite repul sion for an arbi t rary c on- 
fining potential was present in ( Guan L. et all [20091) . 

In the context of exact solutions, considerable work 
has been done to derive low temperature analytic re- 
sults for BA solvable models. These include the work 
on the free energy of spin c hains at low temperatures un- 
der a small magnetic field (IMezincescu and Nepomechiel . 
Il993l : iMezincescu et al. |, 119931) and the calculation of 
the leading temperature dependent terms in the free 
energy of the massive Hciscnberg model by Johnson 
and McCoy (1972). Later, Filyov et al. (1981) de- 
rived an exact solution to the s-d exchange model ex- 
pressed as a series in term s of the temperature. Follow- 
ing t h e method proposed (IMezincescu and Nepomechiel . 
119931: M ezincescu et al. |, 119931) for small external field 



H < 1, Lee et al, (|Lee et all l2012h solved the TBA 
equations (|58p by the Weiner-Hopf method, where the 
dressed energy potential is simplified as e(k) w k 2 — A 
with the cutoff chemical potential A := fj, + 2Pln2/c + 
cH 2 /Ai: 2 P + cT 2 /6P. The terms in the function A pro- 
vide insights on spin and charge density fluctuations. 
With t he help of this potential and using Sommerfeld's 
lemma (jPathrial . Il996l) , the free energy of the system de- 
fined by Eq. (|59[) is given by 
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(63) 



which gives the universal low temperature form of spin- 
charge separation theory (|61[) with the excitation veloci- 
ties 



27m 1 - 



41n2 

7 



1 - 



61n2 



7 



For the external field approaching the saturation field 
H s , the charge and spin velocities can be derived from 
the relations (|()2"j) . The leading terms in the velocities are 
then found to be 
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The susceptibility values for different values of the chem- 
ical potential are consistent with the field theory pre- 
diction xv s = 6/-K with 6 = 1/2 (|Giamarchil 120041 ) for 
strong repulsion. On the other hand, the TBA equa- 
tions (|58[) show that the spin-spin interaction for the 
system of the polarized fermions with strong coupling 
7 3> 1 can be effectively described by the isotropic 
spin-1/2 Heisenberg chain with a weak antiferromagnetic 
coupling J = — ||yp u (T, H). For the isotropic spin-1/2 
Heisenberg chain, the susceptibility at H = is given 
by X = l/(^ 7r2 ) whi ch coincides with the field theory 
prediction x v s =9 /it ( Lee et all . l2012h . 

Furthermore, for strong repulsion, the pres sure of the 
gas w ith a weak magnetic field is given by ( Lee et all . 
I2OI2D 



P 



2irh 2 



T2 Lia 



(64) 



This result provides the low temperature thermodynam- 
ics which extends beyond the range covered by spin- 
charge separation theory. With the low temperature ex- 
pansion, the pressure (|64[) significantly evolves into the 
thermodynamics of two free Gaussian fields at criticality. 
The result for the free energy at low temperature gives a 
universal signature of TLLs where the leading low tem- 
perature contributions are solely dependent on the spin 
and charge velocities. This opens a way to experimen- 
tally explore how the low temperature thermodynam- 
ics of a ID many-body system naturally separates into 
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two free Gaussian field theories. This may possibly be 
used to test the spin and charge velocities in experiment 
with an ultracold atomic ID Fermi gas in an harmoni c 
trap ( Cheianov and Zvonarevl |2004 Iflecati et al. . 2003). 
where the spin and charge of the ultracold atoms refer 
to two internal atomic hyperfine states and the atomic 
mass density, respectively, as per the theoretical scheme 
to explore spin-charge separation waves in Fig |14l This 
phenomenon has already been experimentally observed 
in ele ctron liquids ( Deshpande et aLl . l2010HJompol et all , 
I2009D . 



IV. FERMI-BOSE MIXTURES IN ID 

The experimental advances in trapping and cool- 
ing ultracold quantum gases have also led to real- 
izations of degenerate Fermi-Bose mixtures with var- 
ious combinations of fermionic and bosonic atoms 
(iTruscott et al 



(IChen et all 120101 iGirardeau and Mingguzzil, 120 071 : 



iHu et all l2006t llmambekov and Demleri. l2006allbt). 
magnetism (Batchelor et all 2005aTGuan et all 2008bt 



llmambekov and Demler. 2006b), correlation func 



tions (iFang et all l2011bt iFrahm and Palacioi. l2005t 



llmambekov and Demler 2006bl) . thermod ynamics 
quantum criticality (|Yin et al 020091120121) . 



and 



such as 6 Li-'Li 



12001 



i Hadzibabic et oil 120021: ] Stan et all , l 200^i). 

( Inouve et all 12004 lOsoelkaus et all 120061: iRoati eTal\ . [ + b ' 



6 Li- 23 Na 
40 K- 87 Rb 



A. Groundstate 

The ID (5-function interacting mixture of Mb spinlcss 
bosons and Mi fermions with spin-up and M2 fermions 
with spin-down is described by the Hamiltonian ([2]) with 
Zeeman term |if(Mj — M 2 ). H is an external mag- 
netic field. The Bethc wave function for this model re- 
quires symmetry under exchange of spatial and internal 
spin coordinates between two bosons or fermions with 
different spin states and antisymmetry for two fermions 
with the same spin state. The BA equations for this 
Fermi-Bose mixture with an irreducible re presen tation 

19711) 



iMi-Af 2 



are (Lai and Yan 



12002ft. 6 Li- 87 Rb (ISilber eta ll |2005l), 173 Yb- 174 Yb 



dFukuhara et all 120091) and 6 Li- 17 4 Yb ([Hansen et al 
20 111) and ^LT 133 Cs (jRepp et all 12013 : iTung et al 



A I 



20131 ). This experimental success opens up a further 



gateway for exploring striking quantum many-body phe- 
nomena though tuning interactions between inter- and 
intra-species of atoms. 

In light of the experimental realizations of Fermi-Bose 
mixtures, various theoretical methods have been used 
to study phases of superfluids and Mott insulators, 
instabilities of collapse and demixing and quantum 
correlations of the ID Fer mi-Bo se mixtures, such 
as t he mean-field ap p roach ( Dad, 
orv (ICazalilla and Hoi 120031: 
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Mathey et all 12004 IMathev 



20101 : iRizzf and Imambekovl 



2003), TLL the- 



Lewenstein et al 
20071: 



methods (|Pollet et al 



20081) 



2004; 



Qrignac et al 



IVarnev et al 



2008 



and numerical 
20081: iTakeuchi and Mori 120071: 



Zuiey et al 



200£ 



The TLL field 



theory (ICazalilla and Hoi . 120031: IMathev et all 12004 



iRizzi et al . 20081 ) predicts that the binary mixtures 



of bosons and spin-polarized fermions with population 
imbalance present competing ordering - i) strong at- 
traction between the two species leads to collapse; ii) a 
strong repulsion leads to demixing; iii) subtle tuning of 
the intra-and inter-particle scattering leads to pairing 
and two-component TLLs. 

On the other hand, the ID Fermi-Bose mixture with 
equal masses of bosons and spin-polarized fermions and 
with the same strength of delta-function interaction 
between boson-boson, boson-fermion and fer mions 
with different spins was so lved a long time ago (|Lail . 
Il974ri lLai and Yangl . Il97ll ). The exactly solved model 
provides a benchmark towards understanding var- 
ious quantum many-body effects in ID Fermi and 
Bose mixtures. Recently, particular theoretical in- 
terest has been paid to the groundstate properties 



where M = M 2 + Mf,. In these equations kj, with j = 
1, . . . , N, are the quasimomenta of the particles and A a , 
with a = 1, ...,M, are parameters for fermions with 
spin-down and the bosons. At with b = 1, ...,Mj, are 
the parameters for the bosons. 

Lai and Yang (1971) showed by numerically solving 
the set of coupled integral equations that the energy 
of the system for the mixture of bosons and polarized 
fermions is a monotonic decreasing function with respect 
to the ratio of bosons and the total number of particles 
in the system. The groundstat e properties of the sys- 
tem have been further stu died ( Batchelor et all l2005al : 
IFrahm and Palaciosl 120051) . 

The magnetic properties of the mixture of bosons and 
polarized fermions provide further insight into the com- 
peting ordering. Application of the external magnetic 
field to the polarized fermions causes Zeeman splitting 
of the spin-up and spin-down fermions into different en- 
ergy levels. The groundstate can only accommodate 
fermions that are in the lower energy level. Therefore 
it is expected that when the direction of the magnetic 
field is along the spin-up (H > 0) direction, spin-down 
fermions can no longer populate the groundstate. The 
free energy of the strongly coupled mixture is given in 
the form F = —(n — ni{,)H/2 + Eo w here the groundstate 
energy per unit length is given by (|Guan et all l2008bl: 



26 



llmambekov and Demleit l2006al lbl) 
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Here rib is the boson number density. If the ex- 
ternal field exceeds the critical field H c = 8po/c, 
where the pressure per unit length is given by Pq ~ 
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the system enters a 

phase of fully-polarized fcrmions. The susceptibility in 
the vicinity of the critical field H c diverges as 
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(67) 



This van Hove type of singularity is subtly differ- 
ent from the linear field-dependent magnetization in 
the two-component attractive Fermi gas with polariza- 
tion, where the effective interaction between the bosonic 
pairs is weakly at tractive in the Tonks-Girardcau limit 
(|Chen et all l2010f ). 

In contrast, for weak repulsion the groundstate en- 
ergy presents a mean field effect, i.e., the energy per unit 
length is 
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(68) 



The magnetization r 
ccptibility \ ^ 



z « t-(V2H + ^) and the sus- 

47T V v 7T / 

follow from this equation. These 
results show that in the weak coupling regime the square- 
root field-dependent behavior of magnetization emerges 
for finite external field. 

Furthermore, using C FT, Frahm and Palacios 



(|Frahm and Palaciost 120051 ) lave computed the asymp- 
totics of boson and fermion Green's functions for a mix- 
tures of bosons and polarized and unpolarized fermions. 
The Fourier transform of equal time boson Green's func- 
tion has a form n b (k) ~ |fc — fcol^ 6 . The response 
function of fermions follows a form n a (k) ~ sin(fc — 
ko)\k — ko\" f . The exponents rib and Vf are determined 
by finite-size energies and momenta. The presence of 
spin population imbalance of fcrmions gives rise to dif- 
ferent critical exponents. In fact, the exact result indi- 
cates that no demixing occurs in the Fermi-Bose mix- 
ture with a repulsive interaction. A thorough study of 
the ground state properties, including density distribu- 
tions and the single particle correla t ion fun ctions was re- 
ported in ( Imambekov and Demlerl . l2006bl ) . Using exact 
solution with local density approximation in a harmonic 
trap, the density profiles of the Fermi-Bose mixture were 
predicted. In the weakly interacting regime, bosons can 
condense to the trapping centre whereas fermions spread 
out due to the Fermi pressure. For strong repulsion, the 
boson density distribution is extended in a wider region 
and the fermion density shows strong non-monotonous 
behavior, see FigfTSl 



fermions 



bosons 




FIG. 15 Density profile of Fermi-Bose mixture in strong cou- 
pling regime 7 > 1. At zero temperature, fermions sit at 
the wings while the trapping centre comprises of a mixture of 
fermions and bosons. From Imambekov and Demler (2006b). 



B. Universal thermodynamics 

The fully-polarized Fermi-Bose mixture is described by 
the Hamiltonian 



n = 



dx 



2m 



2m f 1 



^,T,t,T,t 



^l^l^b^b + g bf ^l^\^f^ b - yU/*/*/ ~ fJ-b^b 

(69) 

where Vt 7 ;,, 'F/ are boson and fermion field operators, mb, 
mf are the masses, fj,b, /!/ are chemical potentials of 
bosons and fermions, and g b b, 9b f are boson-boson and 
boson-fcrmion interaction strengths, respectively. 

For bosons and fermions of equal mass and equal in- 
teraction strength gbb — gbf, a different set of BA equa- 
tions from (|65|) were derived ( Imambekov and Demlerl 
l2006allbh with 



N 



exp(ifcji) = ei(kj - A a ), JJei(fcj - A Q ) = 1, (70) 

a=l i=l 

where j = 1, . . . , N and a = 1, . . . , M. Here Mb is the 
number of bosons. The BA equations (|65|) with a par- 
ticular choice of fermion spin polarization arc physically 
equivalent to (|70|) . This can be seen from the fact that 
the groundstate energy of the model (|69|) is given by the 
same expression as (|66p . where the spin-down fermions 
are gapfull. 

At finite temperatures, the equilibrium states become 
degenerate. The thermodynamics of the model are deter- 
mine d from the integral equations (jLail . Il974al lYin et all . 
2009) 



e(k) = k 2 -fi f -T Ki (A - k) In (l + e ~^ A ^ T \ dA, 
<p(A) = n f -Hb-T I K\ (k - A) In (l + C - £ W/ T ) dk, 

J — OO 

(71) 
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For fixed temperature T and chemical potential [if, fib, 
the pressure is given by P = In (l + cT' i ^/ T ^ dk. 

In this grand canonical ensemble, it is convenient to 
use the chemical potential /i and the chemical bias H = 
[if — fib to di scuss the zero tem perature phase diagram 
of the model (|Yin et all l2012h . see Fig. OH The phase 
boundaries were determined by analysing the dressed en- 
ergy equations obtained from the TBA equations (|7II) in 
the limit T — s- 0. The critical line 
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(72) 



separates the mixed phase and the pure fermion phase. 
Here dimcnsionless units have been used, i.e., H = H/eq 
and flf = /i//eo with eo = c 2 . In a harmonic trap, 
this phase diagram can be presented within the LDA, 
e., fj,l(x) + mu^i 2 / 2 = M°(0) and fjP f (x) + muj 2 f x 2 /2 



^f(O). It was found ( Imambekov and Demlerl . l2006bl : 
lYin et all I2009D that for both strong and weak inter- 
actions bosons and fermions coexist in the central part 
and fermions sit at the wings. For weak interaction, 
bosons can condense into the centre while the fermions 
spread out in the whole trapping space due to Fermi pres- 
sure. However, for the strong interaction regime, the 
Fermi density shows strong non - monotonous beh aviour 
(|lmambekov and Demlerl . l2006bl lYin et all l2009h . The 
significant feature of correlation functions of the mix- 
ture in the Tonks-Girardeau l imit w ere discussed in de- 
tail ( Imambekov and Demler . 2006b). In particular, the 
Fourier transform of the Bose-Bose correlation function 
is governed by the TLL parameter Kb via n b (k) ~ 
\k\~ 1+1 /(' 2Kb ^ for k — >• and the Fourier transform of 
the Fermi-Fermi correlation function has singularities at 
k = kf,kt + 2kb. The discontinuity at kf + 2kb indicates 
an interaction effect. 

The TBA equations (|71j) have been used to explore 
scaling behaviour of the thermodynamics in the Fermi- 
Bose mixture. The universal l eading order tem perature 
corrections to the free energy (|Yin et all I20H 
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indicate a collective TLL signature at low temperatures. 
In the strongly repulsive regime for H <C 1 
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As already remarked, the TLL description is incapable 
of describing quantum criticality since it does not con- 
tain the right fluctuations in the critical regime. For the 
physical regime, i.e., c> 1, or T/eq <^ 1, the pressure is 
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FIG. 16 Phase diagram in the fi — H plane. Three distin- 
guished phases result from varying the chemical potential and 
the chemical potential difference H — Hf—fn: the pure boson 
phase for H < and /i > H/2; the pure fermion phase below 
the phase boundary H c in the region /i > —H/2 and the mix- 
ture of bosons and fermions above the phase boundary H c in 
the region H > 0. From Yin et al. (2012). 
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FIG. 17 Contour plot of entropy S vs chemical potential 
from the exact solution, showing quantum criticality driven 
by chemical potential for H = O.leo. The crossover tem- 
peratures (white squares and triangles) separate the vacuum, 
TLLf and TLLm from the quantum critical regimes. From 
Yin et al. (2012). 



Here the functions (3 and A are determined by 
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The function tp(A ) can be obtained from ([71]) by iteration 
(|Yin et q/J . l20ia . 

The equation of state (|74^) can be used to explore 
the critical behaviour of the model. E.g., the entropy 
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in Fig. [17] shows that the TLL is maintained below a 
crossover temperature at which the linear temperature- 
dependent entropy breaks down. In Fig. \T7\ TLLp de- 
notes a TLL of fully-polarized fermions and TLLm de- 
notes a a two-component TLL of the mixture of bosons 
and fermions described by ([73| . As the temperature is 
tuned over the crossover temperatures the scaling func- 
tion of thermodynamic properties give rise to the free 
fcrmion universality class of criticality that entirely de- 
pends on the symmetry excitation spectrum and dimen- 
sionality of the system. The equation of state pressure, 
determined by (|74p . contains universal scaling functions 
which control the thermodynamic properties in the quan- 
tum critical regimes. Near the critical point, the thermal 
dynamical properties can be cast into universal scaling 
forms, e.g., (|53"|) and (|54[) for a free Fermi theory of crit- 
icality, i.e., with d = 1, z — 2 and v = 1/2. The detailed 
analysis of quantum criticality of the Fermi-Bose mixture 
has been presented in Yin et al. 2012. As for the Fermi 
system, with the help of the exact solutions the critical 
properties of the bulk system can be mapped out from 
the density profiles of the trapped Fermi-Bose mixture 
gas at finite temperatures. 



V. MULTI-COMPONENT FERMI GASES OF 
ULTRACOLD ATOMS 

A. Pairs and trions in three-component systems 

A pseudo spin- 1/2 system of interacting atomic 
fermions has been experimentally realized by loading 
atoms within two lowest hyperfine levels, i.e., states 
|1) and 1 2). The interaction strength can be con- 
trolled by tuning the scattering length through Fesh- 
bach resonance. Weakly bound molecules exist in the 
phase where the scattering length is small and posi- 
tive. These molecules can form a BEC. The tunabil- 
ity of the scattering length across the Fcshbach reso- 
nance leads to divergent scattering length. As a re- 
sult the interactions can be effectively enhanced. In 
the strong interacting limit, these molecules may con- 
tinuously transform into BCS pairs such that the sys- 
tem reaches the BEC-BCS crossover (|Bartenstein et 



12004 iRegal et all 12004 IZwierlein et ctZ.I . l2005h . Univer 



sal many-body be haviour is expected in the crossover 
(unitarity) regime (jHeiselberd . l200lb IHoI l200l . 



It is of a great interest that the third pseudo spin 
state 1 3) is added to t h e two-component Fe r mi ga s 
dBartenstein et all 120051: iModawi and Leggettl . Il997lh 
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FIG. 18 Three-component fermionic atoms in an optical lat- 
tice. The colour pairing phase occurs for on-site attraction 
U < U c , whereas the trionic state occurs for U > U c , From 
Rapp et al. (2007). 



three low sublevels denoted by |1), \2) and |3) can form 
three possible pai r s |1) + |2), \2) + |3) and |1) + |3) 
( Bartenstein et o*l l2005h . ' Degenerate three-component 
fermions with tunable interparticlc scattering lengths 
«12, 0-23 and ai3 open up the possibility of novel 
many-body phenomena. 

Specifically, strongly attractive three-component 
atomic fermions can form spin-neutral three-body bound 
states called 'trions'. A degenerate Fer mi gas of atoms in 
three different hyperfine states of 6 Li ( Ottenstein et all 
120081) has been created at a temperature T = 0.37Tp, 
where Tp is the Fermi temperature. The spin state 
mixture of ultracold fermionic atoms has been further 
used to study subtle physics of three-body recom- 
bination, atom-dimer scattering, the atomic Efimov 
trimer, associati on and disassoc ia tion o f the atom-dimer 
etc dBraaten et all 2009 : Ferlaino et 



collis ions 
20091: 



Huckans et all 12009 : Ottenstein et al 



Pollack et all 120091: ISpiegelhalder et al. 



2008: 



2009; 



Wenz et all 120091: IZaccanti et all |2009h . In partic 



As a consequence, BCS pairing can be favored by 
anisotropics in three different ways, namely atoms in 



ular, the three-component interacting fermions have 
received considerable interest in the study of the 
quantum phas e transition from a colour super fl uid t o 
singlet trions (jCherng et all l2007t IRapp et all 120071) . 
see Fig. [TSJ This study also sheds light on colour 
superconductivity of quark matter in nuclear physics. 



1. Colour pairing and trions 

Loading three-component fermions with contact inter- 
action on a ID optical lattice, the system, i.e., the ID 
multi -component Hubbard mod el, is no longer BA solv- 
able (|Chov and Haldand . Il982h . In terms of bosoniza- 
tion approach , it ha s been found ( Azaria et all 120091: 
ICapponi et all |2008|) that the low-energy physics of the 
ID three-component Hubbard model shows the forma- 
tion of three-atom bound states, i.e., a quantum phase 
transition from a colour superfluid to the singlet trion 
state. Motivated by such exotic phases, the ID integrablc 
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three- component Fermi gas w ith delta-function interac- 
tions ([Sutherland I1968L Il975t iTakahashil Il970bt lYaneL 
119701 ) has been studied to give a prec ise understanding of 
colour pairing and the trionic st ate ( Guan et all l2008al 
iHe et aU l20ia iLiu et all l2008al) . 

The first quantized many-body Hamiltonian of the 
three-component <5-function interacting Fermi gas is as 
defined in Eq. ([2"]l. with an additional Zeeman energy 
term E z — J2%=i N v e l z (tf B , B). In this system, the 
fermions can occupy three possible hyperfine levels (|1), 
1 2) an d |3)) with particle number A^ 1 , N 2 and -/V 3 , respec- 
tively (ISutherlandl . Il968t iTakahashl Il970bt lYand . 1 19701) . 
The Zeeman energy levels e l z are determined by the mag- 
netic moments \x B and the magnetic field B. By con- 
vention, particle numbers in each hyperfine states satisfy 
the relation TV 1 > N 2 > iV 3 . Thus the particle numbers 
of unpaired fermions, pairs, and trions are respectively 
given by Ni = N 1 - iV 2 , N 2 = N 2 - N 3 and N 3 = N 3 
for the attractive regime. The unequally spaced Zeeman 
splitting in the three hyperfine levels can be characterized 
by two independent parameters Hi = 1 — e z (fi B ,B) and 

Hi = e z (fj, B ,B) ~ e. Here e = J2a=i e z/3 i s an average 
Zeeman energy 

For the ID three-component Fermi gas, the energy 
eigenspectrum is given in terms of the quasimomenta 
{ki} of th e fermions v i a Eq. (JTDJ) , which satisfy the BA 
equations ( Sutherland! . 1 19681 ) 
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(75) 



Here j = 1, . . . , N, a = 1, . . . , Mi and m = 1, . . . , M2- 
The parameters {A Q ,A m } are the rapidities for the in- 
ternal hyperfine spin degrees of freedom. It is assumed 
that there are M2 fermions in state |3), Mi — M% 
fermions in state |2) and N — Mi fermions in state 
|1). For the irreducible representation [3^32^21^1] , a 
three-column Young tableau encodes the numbers of un- 
paired fermions, bound pairs and trions given by N\ = 
N - 2Mi + M 2 , N 2 = Mi - 2M 2 and N 3 = M 2 . 

In the attractive regime, the BA equations (|75|) 
admit complex string solutions for kj, i.e., three- 
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For arbitrary spin polarization, (i) there are A^3 spin- 



neutral trions in the quasimomentum k space accompa- 
nied by A^3 spin bound states in the A-parameter space 
and A^ 3 real roots in the A-parameter space, (ii) N% BCS 
bound pairs in k space accompanied by N 2 real roots in 
A space and (iii) Ni unpaired fermions in k space. These 
root patterns are depicted in Fig. [HA 
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FIG. 19 Schematic root pattern for 3 trions, 4 colour pairs 
and 6 excess single fermions in the groundstate. For strongly 
attractive interaction, the unpaired and paired quasimomenta 
can penetrate into the central region occupied by tightly 
bound trions. From Lee et al. (2011). 



For weak attraction, we redefine the polarizations pi = 
N [ /N with i = 1,2,3 and the energy EL 2 /N 3 = e (-y) 
with dimensionless parameter 7 = cL/N. The ground- 
state ener gy fol lows from the BA equations (|75|) as 
dLee et al.U201ll) 



eo(7) 



loo Ion loo 
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+27 [pip 2 + pip 3 + p 2 p 3 ] + 0(7 2 ). (76) 

This r esult was al s o obt ained from the Fredholm equa- 
tions ( Guan et al 1 120121) . The leading order of the in- 
teraction energy gives a two-body mean field interaction 
energy between the particles with different internal spin 
states. The kinetic energy part indicates three free-Fermi 
gases of different species. 

For strong attraction (I7I 3> 1) these charge bound 
states are stable and the system is strongly correlated. 
The corresponding binding energies of the charge bound 
states are given by e r = h 2 c 2 r{r 2 — l)/(24m). Explicitly, 
the BCS pair binding energy is £b = h 2 c 2 /(4m) and the 
binding energy for a trion is et = h 2 c 2 /m. Without loss 
of generality, the numbers of trions N3 , pairs N 2 and un- 
paired fermions Ni are assumed to be even. The explicit 
root patterns to the BA equations l[75)) for trions, pairs 
and single fermions arc then 
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where i = 1, . . . , JV3, j = 1, . . . , N 2 and i = 1, . . . , JVi, 
see Fig. The real parts of the bound states have a 
hard-core signature, explicitly, 

(2nW + lk 



(2n< 3 > + 1)tt 

W 3L " 3 ' 

(2n ( V + 1)tt 
ke ~ 7 ai, 



A; 



2L 



-«2, 



30 



where tjM = -N r /2, -N r /2 + 1, . . . , N r /2 ~ 1 with r = 
1, 2, 3. Here a r = (1 + A r + A 2 ) with the function 

A -V^V" 4n. t fl(r-2) ^ 4Ttjg(K-r-l) 

3=1 i=j v J ' i=r+l v ' 

where k = 3 and is the step function. 

This result points to the interference effects among 
the molecule states and excess single fermions. From 
these roots, the groundstate energ y in the thermody- 
namic limit is given explicitly by dGuan et all l2008at 
iKuhn and Foersterl . 12013 : iLee and GuarJ . l201ll ) 
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with k — 3. Here n r = N r /L and N r is the number 
of r-atom molecule states. In general the ID interacting 
Fermi gases with higher spin symmetries have two distin- 
guishing features: (i) mean field theory for the weak cou- 
pling regime, and (ii) strongly correlated molecule states 
of different sizes in the strongly attractive regime. The 
fundamental physics of the model is determined by the 
set of equations (|75|) which can be transformed to gen- 
eralised Fredholm types of equations in the thermody- 
namic limit. The asymptotic expansion solution of the 
Fredholm equations for arbit rary component F ermi gas 
has been thoroughly studied ( Guan et all |2012[ ) . 



2. Quantum phase transitions and phase diagrams 

The groundstate energies (|?6")) and (fT9"| provide full 
phase diagrams in the H\ — Hi plane for the weak and 
strong coupling regimes, see Fig. [501 The fields Hi and 
Hi are determined through the relations ( Guan et all 
20081 IKuhn and Foersterl . [20121 ) 
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with the constraint condition n = ri\ +2n2 + 3n3. For the 
strong coupling regime in the absence of Zeeman split- 
ting, i.e., H\ — Hi = 0, trions form a singlet groundstate. 
However, the Zeeman splitting can lift the SU(3) degen- 
eracy and drive the system into different phases. For 
small Hi , a transition from a trionic state into a mixture 
of trions and pairs occurs as H 2 exceeds the lower critical 
value -fff 1 . When H 2 is greater than the upper critical 
value -fff 2 , a pure pairing phase takes place. Trions and 
BCS pairs coexist when H^ 1 < H 2 < H% 2 . These criti- 



cal fields are derived from the rel ations (|80|) ( Guan et all 
l2008at IKuhn and Foersterl . l2012h 
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FIG. 20 Ground state energy vs Zeeman splitting for (a) 
strong interaction |c| = 10 and (b) weak interaction |c| = 0.5 
with n = 1. The figure reveals a trion phase C, a pairing 
phase B, an unpaired phase A and four different mixtures of 
these states. Good agreement is found between the analytical 
critical fields (black lines) and the numerical solutions (white 
lines) of the TBA equations (|82|l . The pure trion phase C 
is present in the strong coupling regime, whereas the colour 
pairing phase is favoured in the weak coupling limit. From 
Kuhn and Foerster (2012). 



The phase transitions from the colour paired phase B to 
the mixed phase A + B of pairs and unpaired fermions 
induced by increasing Hi are reminiscent of those for 
the two-component system discussed in Sec. IIII.AI This 
mixed phase consisting of the BCS-pairs and unpaired 
fermions is referred to the FFLO phase, see Fig. I2TH 

For small Hi, a phase transition from a trionic into 
a mixture of trions and unpaired fermions occur as Hi 
increases. Using the relations (|8"tj)) , the trionic phase with 
zero polarization forms a singlet groundstate of trions 
when the field H < iff 1 , whereas when Hi is greater 
than the upper critical value H^ 2 all trions are broken 
and the state becomes a normal Fermi liquid, see Fig. 1201 
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In addition, for a certain regime of Hi and H 2 , there is a 
phase transition from the trionic state into the mixture 
of trions, pairs and unpaired fermions, see Fig. 1201 

In the weak coupling regime, the critical fields can be 
obtained from the relations 
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We see that either a mixture of BCS pairs and unpaired 
fermions or a mixture of trions and unpaired fermions or 
a mixture of trions, pairs and unpaired fermions popu- 
lates the groundstate for certain values of H\ and H 2 . 
These asymptotic results indeed agree well with the 
full phase diagram determ ined from numerical solutions 
(|Kuhn and Foersterl . 120121) . see Fig. HOI It is interest- 
ing to note that the pure paired phase can be sustained 
under certain Zeeman splittings. In this phase, the two 
lowest levels are almost degenerate for certain tuning of 
H\ and H2 ■ The persistence of this colour pairing phase 
is relevant for the study of phase transition between the 
colour BCS-pairing phase and the state of trions. 

In the thermodynamic limit, the grand partition func- 
tion Z = ti-(c~ H / T ) = cr G l T is given in terms of the 
Gibbs free energy G = E-fiN-HiNi-H 2 N 2 -TS where 
the chemical potential fi, the Zeeman energy Ez and the 
entropy S are given in terms of the densities of unpaired 
fermions, charge bound states, trions and spin-strings, 
which are all subject to the BA equations ([75)1 . The equi- 
librium states are determined by minimizing the Gibbs 
free energy, which gives rise to a set of coupled nonlinear 
integral equations - the TBA equations for the dressed 
energies £ a (a = 1. 2 3) (|He et alL l201li E ee and Guanl . 
120111: ISchlottmannl . [l993i Il994f) . In the thermodynamic 
limit, the pressure p is defined in terms of the Gibbs en- 



ergy by p = —(dG/dL), including three parts, p^\ p^ 
and p^ 3 \ for the pressure of unpaired fermions, pairs and 
trions, respectively, with 
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Here we have set the Boltzmann constant ks = 1. 

The quantum phase transitions in the model with Zee- 
man splitting may be analyzed via the dressed energy 
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FIG. 21 The fj, — H phase diagrams for (a) equally spaced 
Zeeman splitting and unequally-spaced Zeeman splitting (b), 
(c) and (d). From He et al. (2010). 



equations, which are obtained from the TBA equations 
in the limit T — > as 



e< 2 )(A) 
e«(fc) 



3A 2 - 2c 2 



3/i- A' 2 *e (1) (A) 
- [A'i + A 3 ] * e (2) (A) - [A 2 + A 4 ] * e< 3 >(A), 
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2A 
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-K 2 * e 2 (A) - [Ax + A 3 ] * e< 3 >(A), 

k 2 - fj, - H 1 - Ai * e (2) (fc) - A 2 * e (3) (fcj(82) 



Here we denote A„ * e^ a \x) = J_q K n (x — y)e^ a \y)dy . 
The integration boundaries Q a characterize the "Fermi 
surfaces" at e^ a '(Q a ) = 0. The chemical potential and 
magnetization are determined by Hi, H 2 , gir> and n 
through the relations — ^ = n, — -§§^ = n>i, —-§§^ = 
n 2 . The dressed energy equations (|82[) indicate effective 
interactions among trions, pairs and single fermions. If 
we denote effective chemical potentials /x* = \i + et/3 
for trions, /i b = fx + £b/2 + A 2 /2 for bound pairs and 
/i u = fx + Hi for unpaired fermions, the energy transfer 
relations among the binding energy, the Zeeman energy 
and the variation of che mical potentials betw een different 
Fermi seas are given by (jGuan et all l2008al) 



Hi = 2c 2 /3 + (/i u - ft), H 2 = 5c 2 /6 + 2(// b - 

Hi - H 2 /2 = c 2 /4: + ( At u - M b ). (83) 

These equations determine the full phase diagram and 
the critical fields triggered by the Zeeman splitting Hi 
and H 2 . Indeed, the relations ([55)1 give the results (|T6) 
and (|79[) in the weak and strong coupling regimes. In 
the phase diagrams Fig. [5U the phase boundaries can 
be also analytically determined by analyzing the band 
fillings through the dressed energy equations 
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3. Universal thermodynamics of the three-component fermions 

At low temperature regimes, the TBA equations pro- 
vide universal thermod ynamics of the attractive three- 
component Fermi gas ( He et all . |2010[ ). In the grand 
canonical ensemble, the unequal spaced Zeeman splitting 
lead to various phases in the fx — H plane, see Fig. [21] 
This fi — H phase diagram can be determined by the 
dressed energy equations (|82|) . The quantum phases at 
zero temperature persist due to the nature of collective 
motion (forming TLL phases) for a certain temperature 
range. Although there is no quantum phase transition in 
ID many-body systems at finite temperatures, quantum 
criticality leads to a crossover from a relativistic disper- 
sion to a nonrelativistic dispersion between the TLL and 
the quantum critical regime. The nature of the TLL 
physics is revealed from the universal leading tempera- 
ture corrections to the free energy in the different phases 
demonstrated in Fig. [20l i.e., 



F f=s En — 



ttT 2 
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(84) 



where the sum involves a term for each cluster component 
in the phase. E.g., + — for phase B + C. For strong 



attraction, v r 



h7rn r 



(I + ^A r + j,A 2 r ) with r = 1,2,3 
are the velocities for unpaired fermions, pairs and trions, 
respectively. Here the function A r are given by (|78[) . 

On the other hand, the pressure p^ of trions, pairs 
and excess fermions can be obtained in an analytical 
manner using the polylog function in the strong attrac- 
tive regime, with 



(85) 



for a = 1, 2, 3. Up to a few leading order terms, the 
functions A^ are 
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A (3) = 3/i + 2c 2 - 



,(2) - 



,(3) 



3|c| 



where J a = 2p^/(a\c\). 

The total pressure p — 52 a =i provides a high preci- 
sion equation of state through iterations with ([87]) . The 
thermodynamics and critical behaviour can be worked 
out in a straightforward manner in terms of the polylog 
function. The equation of state ((85)) covers not only the 
zero temperature result of the three-component strongly 
attractive Fermi gas but also the TLL thermodynamics. 
This result opens up further study of quantum criticality 



with respect to the phase transitions when the param- 
eters drive the system across the phase boundaries of 

Fig. mi 

For the repulsive regime, the thermodynamics of the 
three-component Fe rmi gas is determined by another set 
of TB A equations (jHe et all . 1201 it ISchlottmannl Il993l 
|l994f ). In this regime, spin-charge separation is a hall- 
mark of the ID three-component Fermi gas. In the 
low-lying excitations, interacting particles "split" into 
spins and charges as the temperature tends to absolute 
zero. The collective motion of fermions with only spin or 
charge, called spinons and chargons/holons (the antipar- 
ticle of a chargon), which have different velocities. The 
three-component Fermi gas has U(l) x SU(3) symme- 
try that l e ads t o two sets of spin waves. It is shown 
( He et all . 1201 lh that the low temperature thermody- 
namics of such a gas naturally separates into free Gaus- 
sian field theories for the U(l) charge degree of freedom 
and two spin rapidities. The free energy gives a universal 
low temperature TLL behaviour, namely 
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ttT 2 (C, 



(87) 



6 \ v s v c 

The spin and charge velocities can be derived (jHe et all , 
l201lh from the relations v c = e' (fen) /2np c (ko) and v s = 
<f>i ' (An) /2np s (An). For three-component fermions, 
there are two spin velocities v s \ and v S 2, where v s i = v S 2 
for pure Zeeman splitting. The central charge for the spin 
part is C s = 2 and for the charge part C c = 1. The rea- 
son for the value C s = 2 is because the 5T/(3)-invariant 
fermion model has two spin "Fermi seas" whose depen- 
dence on H are equal, i.e., we considered the case where 
Hi = Hi = H . 

Following the Wiener-Hopf method developed for 
the study of the thermodynamics of Heisenberg 
spin chains with SU(2) and S U(3) symmetries 
dMezincescu and Nepomechid . Il993t iMezincescu et all , 
11993ft . the groundstate energy of the gas in a weak mag- 
netic field is given by E = ±n 3 n 2 (l - ^ - 2n_Ln3^j _ 

! ' . Explicitly, the spin and charge 



4n 2 7r 4 

velocities 
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are derived in the strong repulsive regime. The spin ve- 
locity tends to zero while the charge velocity tends to 
the Fermi velocity as the interaction strength c — > 00. 
The low temperature properties of the equation of state 
follow from 




,A/T 



(88) 



where the potential function A = p,+2nP y 6 ^ c + 5^ 

The thermodynamics obtained from this 
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pressure covers the universal thermodynamics of the TLL 
given by ([57]). 



B. Ultracold fermions with higher spin symmetries 

1. Bosonization for spin-3/2 fermions with SO(5) symmetry 

Spinor Bose gases with spin-independent short-range 
interactions have a ferromagnetic groundstatc, i.e., 
the groundstatc is always fully polarized. In con- 
trast to the spinless Bose gas, the spinor Bose gases 
with spin-exchange interactions can display a differ- 
ent groundstate, i.e., either a ferromagnetic or an an- 
tiferromagnetic groundsta te solely depending on the 
spin-exchange inter a ction (|Hol . Il998t I Ho and Yid . l2000t 
lOhmi and Machidal . Il988l) . In this regard, the ID 
spinor Fermi gases with a short-range delta-function in- 
teraction and spin-spin exchange interaction are par- 
ticularly interesting due to the existence of various 
BCS-like pairi ng phases of quan t um liquids associat ed 
with the BA (lEssler et all 120091: iKuhn et all |2012allbl : 
iLee et all l"2009t IShlvapnikov and Tsvelikl . 1201 11 ). Large 
spin fermionic systems also exhibit many new phases of 
matter which do not appear in the usual spin- 1/2 sys- 
tems. In particular, spin-3/2 systems possessing a generic 



SO (5) symmetry have attracted much theoretical inter- 
est (ILecheminant et ~al\. 120051: IWul . l2005ll2006HWu et all 
12003 ). In these systems, s-wave scattering acquires inter- 
action in total spin singlet and quintet channels. The 
two interacting channels present a hidden 50(5) symme- 
try without fine tuning. The spin-3/2 systems exhibit a 
quintet pairing phase with total spin-2 and quartetting 
order as a four-fermion counterpart of the Cooper pair- 
ing. 

The Ham ilton i an describing these 5Q(5)-in variant sys- 
tems reads (|Wul . I2005L 120061: IWu et al. I. I2003D 



H 



^( r )(-_A*-, 



+ . go P o yr)P ,o(r) + g 2 £ P\ , f (r)P 2 , £ (r) \ (89) 

t=±1,±\S) J 

where d is the dimensionality and fj, is the chemical po- 
tential. The operators Pq and P\ m denote the spin sin- 
glet and quintet pairing operators given by Pp m {r) = 

E^(|f;^H||;^)V'l(r)^(r) with F = 1,2 and 
to = —F, —F + 1, . . . , F. Using the Dirac m atrices, the 
5*0(5 ) algebra has been constructed explicitly ( Wu et all . 
The spin singlet pair and quintet pair have been 
constructed in terms of the 50(5) scalar and vector op- 
erators. The 50(5) generators commute with the Hamil- 
tonian (|89|) . The models restore 5/7(4) symmetry when 
fjo = .92- The Fermi liquid theory of 50(5) symmetry 
describes different competing orders in the spin-3/2 con- 
tinuum and lattice models. The lattice version of the 




3:g = 3g, 



FIG. 22 Phase diagram of the ID spin-3/2 Fermi gas with 
SO(5) symmetry in terms of the singlet and quintet interac- 
tion channel parameters go and g2- Various competing orders 
of singlet and CDW pairing, as well as QROS and CDW of 
quartets are formed. From Wu (2005). 



Hamiltonian (|89|) can be viewed as the one-band gen- 
eralized Hubbard model exhibiting 50(5) symmet ry at 
arbitrary filling and 5 Q(7) symmetry at half-filling (jWul . 
120051 : IWu et all 120031) . 



The 50(5) invariance can apply equally well in the 
one-dimensional continuum model and the lattice model 
( Lecheminant et al 1 120051: IWuL [20051 120061) . Writing the 
left and right moving currents in terms of t he 50(5) 
scalar, vector and tensor currents (jWul . l2005ri . the low 
energy physics can be described by an effective Hamil- 
tonian density, which can be treated by bosonization. 
The result indicates two spin gap phases and various or- 
der parameters, see Fig. [22] The TLL phase exists in 
the repulsive region go > g 2 > where K c < 1. The 
quartetting phase (B) has two orders - the quasi-long 
range ordered superfluidity (QROS) [B.l) and charge 
density wave (CDW) of quartets (B.2). The pairing 
phase (C) separates into a spin singlet pairing phase 
(0.2) and dimenerization of spin Peierls order (Ol). 
Here the competition between the quartetting and pair- 
ing phases is characterized by the Ising duality. In 
this scenario, the low energy physics of ID arbitrary 
half-spin fermio ns of ultracold atoms has been studied 
by bosonization (ILecheminant et al. . 2005: N onne et all 
120101 : Is zirmai and Lewensteinl . 120111 ). See also a recent 
study on competing orders in ID h alf-filed fermi o nic ul - 
tracold atoms in an optical lattice ( Nonne et all 1201 ll ) . 
Two different superfluid orders, a confined BCS pairing 
phase and a confined molecular superfluid, were found 
for F = N — 1/2 fermionic ultraccold atoms. 
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2. Integrable spin-3/2 fermions with SO(5) symmetry 

Controzzi and Tsvelik (2006) write that "although high 
symmetries do not occur frequently in nature, they de- 
serve attention since every new symmetry brings with it- 
self a possibility of new physics" indicates a perspective 
of large spin fermionic atoms. Fortunately the ID realiza- 
tion of spin-3/2 fermions with 5Q (5) symmetry is exactly 
solve d under a suit a ble co ndition ( Controzzi and Tsvelikl . 



120061 : I Jiang et~al\ . l2009f) . In particular, Jiang et al. 
(2009) found that the Hamiltonian ([89]) reduces to an 
integrable many-body Hamiltonian 
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N 



3 = 1 J 
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(en + c 2 Sj-Sf) 5(xj 



xt) (90) 



on a line Co/c 2 = —3/4. Here the coupling constants are 
given by c = (go + 2<? 2 )/3 and c 2 = (g 2 - go) /3. The 
integrable Hamiltonian (|M|) possesses 5*0(5) symmetry. 
But the individual spin components are no longer con- 
served. However, 7i = Ns + Ns + Ni + Ni + N_s , 

2 2 2 2 2 

In = Ns — N_ 3 and 7 3 = Ni — N i are three indepen- 

2 2 _ I 2 

dent conserved quantities. The integrability is guaran- 
teed by the two-body scattering matrix 



t, - h - i# 



kj — k[ + i 



P + P. 



kj k[ 



3c" jl 
2 



jl kj - k[ 



which satisfies the YBE 



Si 2 (fci - k 2 )S 13 (k 1 - fc 3 )5 23 (fc 2 - fc 3 ) = 
5 23 (fc 2 - fcs)5i 3 (fci - fc 3 )5i 2 (fc 2 - fc 3 ). (91) 

In the above equation, PJJ 1 is the projection operator onto 
the spin m channel. 

The solut i on ha s been derived in terms of the BA 
( Jiang et q/J . l2009t l. In the repulsive regime, the low en- 
ergy physics can be described by the spin-charge separa- 
tion theory. The elementary spin excitations, including a 
spin-3/2 spinon, a neutral spinon an d a spin-1/2 spinon, 
has been studied (jJiang et al ., 2009), see Fig. [23] Here 
we see that the charge excitations (a) indicate a particle- 
hole type. The spin excitations (b) and (c) involve two 
real A and \i holes, respectively. The spin excitations (d) 
give two string-2 hole excitations. These spin excitations 
arc different from the case of SU(4) symmetric fermions. 
For an attractive interaction, competing pairing orders 
lead to quantum phases of pairs and quartets. However, 
the BA equations give very complicated root patterns. 
The study of the attractive integrable SO(5) symmetric 
model still remains an open problem. 

Although integrable spin systems with higher symme- 
tries have been extensively studied in the literature, ex- 
actly solved models of fermionic ultracold atoms with 
higher spin s ymme try are still restri cted to 577 (2s + 1 ) 
((Sutherland Il968l ) and Sp{2s + 1) (jjiang et all l201ll) . 
For spin-dependent interaction, the spin exchange inter- 
action between the i- and j-th particles can be written 
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FIG. 23 Charge and spin excitations: (a) charge particle-hole 
excitations; (b) spin excitations of two real A-holes in the spin- 
A sector; (c) spin excitations of two real /i-holes in the spin-/i 
sector and (d) spin excitations of two string-2 A holes in the 
spin- A sector. From Jiang et al. (2009). 



as a summation of spin projection operators P™ in the 
channels with even total spin m = 0, 2, . . . , 2s. The in- 
tegrable Sp(2s + l)-invariant models of atomic fermions 
are artificially written as 



JV 



i:) 2 



i=l 1 i^Lj 



where the interaction potential reads 



(92) 



^' = (-l) 2s+1 [ S +l/2-(-l) 2s ] cP t 



m=2,4,...2s 



The BA equations for the model (|92|) have been de- 
rived (jjiang et al\ . 1 2 llh . In addition, Melo and Martins 
(2007) derived the BA solution of a many-body prob- 
lem of interacting spin-s particles with an arbitrary U{1) 
factorized 5-matrix. 

The stability of spinor Fermi gases in tight waveguides 
was discussed in del Campo et al. (2007). 



C. Unified results for SU(k) Fermi gases 

Fermionic alkaline-earth atoms provide unique oppor- 
tunities to study exotic many-body physics in the con- 
text of higher spin symmetry, e. g., SU(k) with k = 
27 + 1, where 7 i s the nuclear spin (ICazalilla et al. . 2009: 
iGorshkov etaH l2010t iHermele e£ al.1 . 12009T) . 2 De Salvo 



2 In this section we use the symbol k to avoid any confusion with 
the number of particles N . 
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et al. (2010) have created a degenerate gas of ultra- 
cold fermionic atoms 87 Sr(F = I = 9/2) in an opti- 
cal trap. Taic et al. (2010) have reported the real- 
ization of a degenerate Fermi mixture of two isotopes 
of ytterbium atoms m Yb (1=1/2) and 173 Yb (1=5/2) 
with SU(2 ) x SU (6) symmetry. More recently they 
( Taie et all . 120121) have successfully realized the 5*1/(6) 
symmetric Mott-insulator state with the atomic Fermi 
gas of 173 Yb in a 3D optical lattice. It was found that 
loading fermions adiabatically into a higher symmetry 
Mott insulating state can achieve lower temperatures 
than the SU(2) symmetry state owing to differences in 
the entropy carried by isolated hyperfine spins. 



These alkaline-earth atoms have a particular filled elec- 
tron shell structure such that their nuclear spins decouple 
from the electronic angular momentum J in these two 
states. This decoupling implies that the nuclear spins 
give the hyperfine spins, i.e., F = I. For these atomic 
systems, the 21+1 hyperfine levels are likely to display 
SU (21+1) symmetry where the s-wave scattering lengths 
are independent of the nuclear spins. Such fermionic sys- 
tems with enlarged symmetries are mo tivated to simu- 
late quantum many- body phenomena ( Gorshkov et all 
120101 ; iHermele et all 120091: IXul |2010|) which may shed 
light on physics of strongly correlated transition-metal 
oxides, heavy-fermion materials and spin-liquids phases. 
In contrast to the weaker quantum spin effect for a com- 
posite object with larger spin in solid state, large hy- 
perfine spins can be e s sential to the quantum magnetic 
state s dCazalilla et al. 



120101 : IWu et all 12003 



sentia l to the quantum mag r 
20091: IGorshkov et all [201of [Wul. 
Xul. |2010D because spin fluctua- 



tions are accommodated in a large number of hyper- 
fine spin states. Large-hypcrfine fermions may also be 
used to st udy Cooper pairing phenomena within hypcr- 
fine spins (Irlo and Yirj. 19991) and quantu m information 
([Daley et all 120081 : IGorshkov et aLl . l2009h . 



The Hamiltonian for the ID TV-bod y (5-fun ction in- 
teracting fermion problem ( Sutherland! . Il968l ) is again 
as defined in Eq. There are now k possible hy- 

perfine states |1), |2), . . . , \k) that the fermions can oc- 
cupy. This system has SU(k) spin symmetry and 
U(l) charge symmetry. For an irreducible represen- 
tation [k n ",(k- l) N »- 1 ,...,2 N *,l Nl ], the Young dia- 
gram has k columns with the quantum numbers iVj = 
N l — iV' +1 . Here iV 1 is the number of fermions at 
the i-th hyperfine level such that TV 1 > iV 2 > . . . > 
N K . The groundstate properties and thermodynamics 
of ID 5'i7(«;)-invariant Fermi gases have been studied 
by means of the BA solutiom e.g., the t hree-component 
Fermi gas (|Guan et aU !2008ailHe et all l2O10t iLiu et al 



2008a! ) . 5 , ;7(4)-invariant spin-3/2 f ermion s ([Guan et al 
20091 : ISchlottmann and Zvvagin 12012 a b c ), and 



comp onent fermions (IGuan et all 120101: iLee and Guan 
20 111: ISchlottmann! . Il993l Il994t lYang and Youl . 12011 



Yin aZ.l . l2011ah 



1. Groundstate energy 

For the groundstate of the S'C/(K)-invariant Fermi gas, 
the generalized Fredh olm equations for c > are given 
by ([Sutherland! . Il968l) 



J-Bi 
m+1 „B a 

r m {k) = Y, / K 1+Sam {k-\)r a (X)d\ (93) 



a—m — l 



— B 



where 1 < m < re — 1 and (3q = l/(27r), ro(k) is the parti- 
cle quasimomentum distribution function, whereas r m (k) 
with m > 1 are the distribution functions for the K — 1 
spin rapidities. The groundstate energy E per unit length 
is given by E = J B g k 2 ro(k)dk. For the balanced case, 
the integration boundaries B m with m > 1 are infinitely 



large . Thus the distributions are given by ([Guan et al 
Soil) r m (A) = 27 JZvrmi^e-^du, where 



fo in (w) sinh — m)|w|c] 
sinh [-i«|o;|c] 



(94) 



with m = 1, . . . , k — 1. 

For strong repulsion, i.e., cL/N ^ 1, the groundstate 
energy of the balan ced re-component Fermi gas follows as 
(IGuan et all 12012!) 
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(95) 



with the constants Z\ = — - [tp(-) + C] and Z 3 = 
re" 3 [C(3, i) - C(3)]. Here ((z,q) and C{z) arc the Ric- 
mann zeta functions, ip(p) denotes the Euler psi function 
and C is the Euler constant. For k — > 00, it is seen that 
lim K _ i . 00 Z\ = limK^yoo Z% = 1. This indicates an insight 
into the hyperfine spin effect the groundstate energy 
([9"5[) reduces to the energy of the spinless Bose gas as 
c — > 00. This result was first noticed by Yang and You 
(2011) by means of the Fredholm equations ([9"3"[) . The 
suppression of the hyperfine spin effect is further man- 
ifest in the groundstate energy per unit length for the 
highly polarized case, which reads 
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(96) 



where m x = M x /L with A'h = YT^Zl N ' i+1 < N - This 
result shows that spin variation does not play an essen- 
tial role in the groundstate due to the strong repulsion. 
However, for small polarization, i.e., a small external field 
lifting the SU(n) symmetry, the integral boundaries B m 
with m > 1 are very large. In this case, logarithmic 
singularities arise in the zero temperature susceptibility. 
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This configuration is drastically changed as the interac- 
tion is decreased. 

For the weak coupling limit, the mean field result for 
the groundstate energy per unit length is 

E = oE^ 2 + 2c E E PiP3 + 0(c 2 ). (97) 



1=1 



i=l j=i+l 



Here pi = N*/N with i = 1,2, ...,k— 1 denote the 
polarizations, and N l is the number of fermions in the 
ith level. The first part is the kinetic energy of the k- 
componcnt fermions whereas the second parts is the in- 
teraction energy. This result is valid for arbitrary spin 
imbalance in the weakly repulsive and attractive regime. 
The higher order corrections have not been obtained. 
For the balanced case, i.e., Ni = N/k, the energy is 
7r 2 n 3 

E = — — - — h c(k — l)n 2 /n + 0(c 2 ) which is the same 

3fy 

as for spinless bosons with a weak repulsion as n — > oo. 
The groundstate properties for the limits c = + and 
c — > oo have been discussed by Schlottmann (1997). 

In the attractive regime, the complex string solutions 
for kj form m-atom bound states up to length 2, . . . , n 

with the binding energy for a bound state = £(£ 2 — 
l)c 2 /12. Such bound states were studied by Gu and Yang 
(1989). A bound state in quasimomentum space of length 

to takes on the form fc™ J = Ai™ -1 ^ + i(m + 1 — 2j')|c'| + 
0(exp(— where j = 1, . . . ,m. The number of bound 
states with length 1 < m < n is denoted by N m . The 
real part is Aq™ 1 " 1 " 1 . The unpaired atoms have real quasi- 
momenta ki. Takahashi (1970b) derived the Fredhohn 
equations for the attractive Fermi gas with an arbitrary 
number of components 

p m (X) = ™A)+EE / #s+m-2r(A- A)p s (A)dA 
r=l s=r J-Qs 

K i-Qs 

+ E / Ks- m (X-A) Ps (A)dA, (98) 

s=m + l J-Qb 

where pi(k) is the density distribution function of sin- 
gle fermions, p m {k) is the density distribution func- 
tion for the m-atom bound state with 1 < to < 
k. The total number of fermions is given by TV = 
ELi m ^m- The integration boundaries Q m , char- 
acterizing the Fermi points in each Fermi sea, are 
determined by n m := = J_q p m {k)dk. The 

groundstate energy per unit length is given by E = 



Em=l f- 



Qn 



mk 
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p m (k)dk. 



For weak attraction \c\L/N <C 1, the two sets of the 
Frcdholm equations ([93]) and (|98| for the repulsive and 
attractive regimes preserve the symmetry 



Pi, 

n<3™ 



Q, 
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c — >• — c. 



(99) 



Indeed, the groundstate energy of the /t-component Fermi 
gas with weak attraction has the same closed form as 
(|97l) with the replacement c — > — c. This means that the 
groundstate energy of the /t-component gas with arbi- 
trary polarization continuously connects at c = 0. 

For strong attraction \c\L/N ^ 1, the bound states 
of different sizes form tightly bound molecules of differ- 
ent sizes. In this regime, all the Fermi momenta of the 
molecules arc finite, i.e., |c| ^> Q m with to = 1, ...,k. 
Here Q\ characterizes the Fermi momentum of the sin- 
gle spin-aligned atoms. Therefore, in this regime the 
strong coupling condition allows one to expand the Fred- 
holm equations (|98| in powers of l/|c|. The ground- 
state energy per unit length of the gas with arbitrary 
polarization in the strong attractive regime is given by 
E = 52m=i( E m ~ n ™ £ b" l) )- Th e energy Ep of the cluster 



state of an ^-atom is given by (jGuan et all 120121 ) 



E n 



n 2 N 3 
3mL 3 
256 



1 



8 



48 



to 3 L 3 |c| 3 



mL\c\ m 2 L 2 \c\ 2 
p3 _ 



5 L 3 I 



_K (ioo) 

[ — G m + G m /15] 



The coefficients F m , G m and G m can be found in Guan 
et al. (2012). 

This closed form for the groundstate energy with ar- 
bitrary polarization is very accurate for a finitely strong 
attraction (for |c| > 5). This high precision groundstate 
energy has bee n given for the two-component a ttrac- 
tive Fermi gas dHe et all 120091 : llida and Wadatil . l2007t 
IZhou et all . |2012[ ). Up to order I/7 2 , the above ground- 
state energy of ID K-component fermions with arbitrary 
population imbalance can be further simplified to the 
general form given in (|79[) . 

From the result (|100p one can obtain full phase dia- 
grams and magnetism at zero temper ature. In particular , 
for the three-component F e rmi gases dGuan et all l2008al : 
iKuhn and Foersteri . 120121 : iLee and Guanl. 1201 ll) spin- 
3/2 Fermi gas with SU(3) s ymmet ry ( Guan et all , 120091 : 
ISchlott mann and Zvvaginl . l2012al). spin-5/2, 7/2 and 
9/2 at tractive Fermi gases (jSchlottmann and Zvvagin , 
2012bl 



2010 



d), and SU ( k) Fermi gases 



Lee and Guanl [20111: Schlottmannl 



Yang and Youl . l201ll : lYhT et aLl l2011af l 



Guan et al 



19931 . 119941: 



The SU(k) symmetry requires each hypcrfinc spin 
state to be conserved. Therefore, there are k chemical 
potentials associated with each spin state. For conve- 
nience in analyzing the quantum phases of the multi- 
component attractive Fermi gas, the Zeemann energy is 
chosen as E z /L = —^£=1 n tHt where Hi is an effec- 
tive external field for the cluster state of size £-atom. 
Here Hi denotes the chemical potential for unpaired 
fermions. The particle numbers arc changed by pN, 
with p the total chemical potential. These effective 
fields result in unequally spacing Zeeman splitting via 
Aj+i^ = —Hi-i + 2Hi — H i+ i. Here H K = because of 
the spin singlet state. Defining the effective chemical po- 
tential of the m-atom bound state by pe = p+(Hi+el)/£, 
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FIG. 24 Phase diagram of the ID integrable spin-3/2 Fermi 
gas with an attractive interaction and pure Zeeman splitting. 
Quantum phases of single excess fermions (I), colour BCS 
pairs (II), trions (III) and quartets (IV) are displayed in the 
/i — H plane for |c| = 1. The shaded area corresponds to the 
vacuum. From Schlottmann and Zvyagin (2012a). 



the result /i£ 



B 

L 



ELiV?) with i = l, 



can be obtained from the groundstate energy. 

The energy-field transfer relati on between H m and the 



effective chemical potentials /x m ( Guan et aLl . l20ld ) 

nn £ . . 

He = £ (fie - + 



(101) 



determines full phase diagrams of the system in terms of 
the effective fields H m and chemical potentials. In the 
special case of pure Zeeman splitting where A^+i.f = 
A for all £, the system has three distinct magnetic 
phases for a strong attractive regime. For weak cou- 
pling, even for pure Zeeman splitting, the phase dia- 
gram is very sophisticated, see Fig. [24] Some subtle 
phase diagrams for the spin-3/2, 5/2, 7/2 and 9/2 a ttrac- 
tive Fermi gases have been stud ied (jGuan et al 1 120091: 
ISchlottmann and Zvvaginl l2012bl lcl). 

Furthermore, for the polarized phase, the cluster states 
form multi-component TLLs in the low energy physics. 
The low-lying excitations are described by the linear dis- 
persion relations uj r (k) = v r (k — k 7 F ) 1 where the velocities 
can be calculated by 



L 1 



mn r r 



d 2 E r 



d 2 L 



(102) 



for a system featuring Galilean invariance. For a 
strong attractive interaction, the velocities for unpaired 
fermions and charge bound state of r-fermions are given 
by v r « ^(1 + ±A r + jtA 2 .) with r = 1, . . . , k. The 
function A r is given by (|75|) . These dispersion relations 
naturally lead to the universal form f or finite-size correc- 
tions to the groundstate energy (|Guan et aUl2010l) . 



E(L, N) w LEI 



nhC 
~~6L 



r=l 



(103) 



The central charge C = 1 for f/(l) symmetry. Here the 
universal finite-size corrections (|103j) indicate the TLL 
signature of the many-body physics. In this case the 
low energy excitations of the system are described by the 
CFT of the Gaussian model. In the next subsection, we 
shall discuss the thermodynamics of these cluster states. 



2. Universal thermodynamics of S'C/(«:)-invariant fermions 

Schlottmann (1993) derived the TBA equations for 
SU(k) fermions with repulsive and attractive interac- 
tion. Lee et al. (2011b) derived a different set of 
TBA equations which are more convenient for analysis 
of thermodynamics and phase transitions for the attrac- 
tive Fermi gas. To understand the thermodynamics of 
this model, it is crucial to separate different physical 
regimes, i.e., 1) groundstate for T — s- 0; 2) TLL phases, 
T < \n — fi c \ and T < \H — H c \; 3) quantum critical- 
ity, T > \fi - fi c \ and T > \H — H c \; 4) high tempera- 
tures, T 3> c 2 . The TBA equations involve an infinite 
number of coupled nonlinear integral equations. At zero 
temperature, the phase diagrams can be analytically or 
numerically obtained through the dressed energy equa- 



the limit T -)■ dGuan et all 


20101 Lee and Guan. 2011 


Schlottmannl. 1993, 1994; 


Schlottmann and Zvvagin 


2012allb: Yane and You. 12011 


Yin et al, 2011a). 



In the strongly attractive regime, the effective fer- 
romagnetic spin-spin coupling constants are given by 
j(r) ~ 7RP (r) for r = 1, 2, . . . , k - 1. is the pres- 
sure for charge r-atom bound states. In this sense, we 
may simply view the non-neutral charge bound state as 
a molecule with spin s = — r, which could flip its 
spin to form the spin wave bound states (spin strings) 
due to thermal fluctuations. However, in the physi- 
cally interesting regime where T « t r , T « Ai+i, 
and 7 3> 1 the breaking of charge bound states and 
spin wave fluctuations are strongly suppressed. The 
spin string contributions to thermal fluctuations in this 
regime can be asymptotically calculated from the TBA 

(r) m + jW T , j{T) , 



i.e., 

(r) 



f: 



It is obvious that 
becomes exponentially small as T — > 0. Thus 



each dressed energy can be written in a single parti- 
cle form e r (k) = ti 2 rk 2 /2m — p,^ + 0(l/j 3 ), where the 
marginal scattering energies among composites and un- 
paired fermions as well as spin- wave thermal fluctuations 
are considered in the chemical potentials p,^. 

With the help of this simplification, the th ermodynam- 
ics at finite temperatures have been given as ( Guan et all 
l2010h 



P 



(r) 



rm „3 
T5 Lis 



-M ^ 



E E 



i ^ 2 j - r 



i{i + r-2j)\c\ 



(104) 
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for r = 1, . . . , N. The total pressure of the system is given 
by p = Yl r =iP^ m Here Li s (x) is the standard polylog 
function. Furthermore, the suppression of spin fluctua- 
tions leads to a universality class of a multi-component 
TLL in each gapless phase, where the charge bound states 
of r-atoms form hard-core composite particles, i.e., the 
leading low temperature corrections to the free energy 
reads 



f ~ -f° FF~ ^2 — 

on z — ' v r 

r—l 



(105) 



The velocity v r of the r-atom bound state is given by 
(|102[) . This result is consistent with the finite-size cor- 
rection result (|103[) . In the above equation fo = — 



X^] 1 n r H r . This result proves the existence of TLL 
phases in ID gapped systems at low temperatures. The 
existence of the TLL leads to a crossover from a relativis- 
tic dispersion to a nonrelativistic di spersion between dif - 
ferent regimes at low temperatures ( Maeda et all |2007|) . 
Linear Zccman splitting may result in a two-component 
Luttinger liquid in a large portion of Zeeman parameter 
space at low temperatures. 



VI. CORRELATION FUNCTIONS 

It is well established that the universality classes of 
critical behaviour of two-dimensional systems are de- 
scribed by a rational conformal field theory (CFT), where 
the Hilbcrt space of states contains a direct sum of irre- 
ducable representations of a Virasoro algebra. Using the 
CFT one can calculate the critical exponents that char- 
acterize p ower law decay of correlation func t ions at larg e 



acterize power law dec ay ot c orrelation functions at larg e 
distance (|Essler et all 120051: iHenkel Il999t IVoitl . Il995l) . 



In general, CFT predicts that the two-point correlation 
function for primary fields with conformal dimension A ± 
is of the form 



(</»(r,y)</»(0,0)) = 



cxp(2iAD£;F2/) 



(vt + iy) (vt — iy) 



2A- 



(106) 



Here r is the Euclidean time (— oo < r < oo, — L < y < 
0) and v is the velocity of light. The conformal dimen- 
sions A^ can be read off from the finite-size corrections 
of low-lying excitations of ID systems via 

Eq-E = ^( X + N++N~), 

Pq-Po = ^{s + N+ - N-) +2ADk F , (107) 

where x = A + + A~ is the conformal dimension and 
s = A + — A~ is the conformal spin. The non-negative 
integers N + and iV - label the level of the descendant. 
AD represents the number of particles backscattered. 

Moreover, it was shown by Affleck (1986) that confor- 
mal invariance gives a universal form for the finite tem- 
perature effects in the free energy by replacing l/L with 



T in the conformal map z = exp(27rw/L). Considering 
a conformal mapping of the complex plane without the 
origin (corresponding to T = 0) onto a strip of width 1 /T 
in the imaginary time direction, the two point correlation 
function at finite temperatures reads 



(0(t,2/)<XO,O)) 



^2iADk F y f+ 



f+(y-ivT)f-(y+ivT) (108) 



, 2A d 



with f±(x) = (ttT/v) 2A / [sinh ^fx\ . For such crit- 
ical phenomena, the critical Hamiltonian of the gapless 
sector exhibits not only global scale invariance but also 
local conformal invariance, i.e., a local version of the scale 
invariance. Thus the critical Hamiltonian can be approxi- 
mately described by the conformal Hamiltonian which is 
described by the generators of the underlying Virasoro 
algebra with central charge C. 

The QISM provides a systematic way to calculate the 
critical exponents for BA integrablc models. Bogoliubov 
et al. (1986) derived explicit critical exponents for the 
XXX and XXZ chains. In this approach, the dressed 
charge matrix Z a p formalism presents a unified frame- 
work to calculate the conformal towers through the finite- 
size correctio ns to the eigenspect rum of multi-component 
BA systems ( Izergin et all U~989) . The universality class 
of critical exponents is uniquely determined by the sym- 
metry of the models associated with the quantum R- 
matrix. Consequently, the asymptotic behaviour of cor- 
relation functions for integrabl e mod els, such as 



penetrable 
metric t 



Bose gas 
J model 



Its et all 119901). the 



Kawakami and Yan: 



the Hubbard model llFrahm and Korepi 



IWovnarovichl Il989t IWovnarovich and Ecklei Il987ft . has 



supc rsym- 
19911) and 



19901. Il99lt 



been studied in the framework of the QISM. In addition, 
exact results for the form factors have been derived for 
several cases dKitanine et all Il999t iKoiima et all Il997t 
ISlavnovt [19891 Il990l). Further developments in the theo- 
retical study of correlation functions have been reported 
(ICalabrese and Cauxl. 120071: ICaux and Calabresel 12006 : 



Caux et all 120071: iGangardt and Shlvapnikovl l2003al lb 



Kitanine et all 120091: iKormos et qZj . |2009h . 



A. Correlation functions and the nature of FFLO pairing 

As remarked in Sec. I, for a system with partial polar- 
ization, the Fermi energies of spin-up and spin-down elec- 
trons are unlikely to match. This leads to a non-standard 
form of pairing known as the FFLO state. The power-law 
decay of the pair correlation n pair oc cos(fcFFLoM)/M Q 
with the spatial oscillations depending solely on the mis- 
match of the Fermi surfaces fcpFLO = 7T i n 1-~ n i) nas been 
identified numerically. This spatial oscillation is a typical 
characteristic of the FFLO pairing. 

The FFLO-like pair correlations and spin corre- 
lation for the attractive Hubbard model were in- 
vestigated by se veral groups via various m e thods , 



such as DMRG (IFeiguin and Heidrich-MeisnerL I2007 
iLiischer e^"aZl . l2008l : lRizzi et a/.l . l2008l:lTezuka and Ueda 
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2008L l20ld), QMC (iBatrouni et all 120081: iBaur et all 



2010l; IWolak et al\. |2010|). me a n fiel d theor y and 



Kaiala et al 



other method s ( Chen and Gaol. 20121: iDatt 



2011; 



Devreese et~ail 1201 it lEdge and Cooperl 120091 



Liu et al 



2009: 



2010; 



2007] . I2008bl ; iParish et al 



20071 : iPei et all 1201(1 IZhao and Lhl \20m . The FFLO 
signature is displayed by the on-site pair correlation func- 
tion O on - s i te (Zi, Zj) := (\l /o|c| 1 c| t c,- 1 1) i n the ID at- 
tractive Hubbard model (jTezuka and Uedal . 120081 ) where 
|\&o) is the eigenstate for the groundstate of the system. 

In Fig. [25l for polarization P = 0, the on-site pair cor- 
relation indicates a power-law decay without changing 
the sign (i.e., spatial oscillation). For low polarization 
P = 0.1, the pair correlation function still has a power- 
law decay but periodically changes its sign. For P = 0.4, 
the oscillation frequency of the pair correlation becomes 
faster due to a large mismatch of the two Fermi energies. 
The wave vector q of oscillations q = AkpFLO- Feiguin 
and Hcidrich-Mcisner (2007) showed that the pair mo- 
mentum distribution function has a peak at the mismatch 
of the Fermi surfaces k = Ak^FLOi i- e -: 



pair 



(109) 



where the pair correlation pf^ 1 oc | cos(Afcp_FLo|^ — 
?n|)/|Z — m\ A ( p \ see Fig. [26l The correlation exponent 
A(P) depends on both the polarization P and inter- 
action strength. Such a FFLO pairing wave number 
was also confirmed by the occurrence of a peak in the 
pair momentum distribution corresponding to the differ- 
ence between th e Ferm i momenta of individual species 
( Batrouni et all 20081: Feiguin and Heidrich-Meisnerl 



l2009t iHeidrich-Meisner et aUboiObt iRizzi et aZ.U2008j) . 

On the other hand, the critical behaviour of ID many- 
body systems with linear dispersion in the vicinities of 
their Fermi points can be described by conformal field 
theory. The critical behaviour of the Hubbard model 
with attractive interaction was investigated by Bogoli- 
ubov and Korepin (1988,1989,1990,1992). As demon- 
strated in previous Sections, the low-energy physics of 
the homogeneous ID Fermi gases with polarization is de- 
scribed by the TLLs of bound pairs and excess unpaired 
fcrmions in the charge sector and f erromagnetic spin- 
spin intera c tions in the spin sector ( Guan et all |20 lot 
IZhao et~ai. I, I2009D . This paves a way to study asymp- 
totic behaviour of various correlation functions by us- 
ing CFT. The method used to study correlation func- 
tions of the spin- 1/2 Fermi gas with attractive inter- 
actio n follows closely the method set out in the litera- 
ture lEssler et al\. 2005 ; Frahm and Korepin . |l990lll99ll 



iKawakami and Yang . 19911 : Wovnarovichl 1989l )~~Conse- 
quently, the asymptotic correlation functions and FFLO 
signature of the ID attractive spin-1/2 Fermi gas have 
been analytically studied by the dressed charge formal- 
ism (|Lee and Guanl . l201ll) . This study can be carried out 
naturally for ID attractive multi-c omponent Fermi gases 
(|Schlottmann and Zvvagml l2012bl) . 
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FIG. 25 Left panels: the on-site pair correlation function 
O on -sitc(zi, Zj) in the Zi — Zj plane for polarization (A^, N{) = 
(a) (20, 20), (b) (22, 18) and (c)(28, 12). Right panels: the axial 
density distribution profiles in a harmonic trap. Numerical 
setting with interacting strength U/t = —4. From Tezuka 
and Ueda (2008). 



In the gaplcss phase, the bound pairs and excess un- 
paired fermions form two Fermi seas which can be de- 
scribed by a two-component TLL, where the spin fluctu- 
ations are strongly suppressed at low temperatures. The 
conformal dimensions of two-point correlation functions 
of the ID attractive spin-1/2 Fermi gas can be calcu- 
lated from the elements of the dressed charge matrix Z 
describing the finite-size corrections for the low-lying ex- 
citations. The long distance asymptotics of various corre- 
lation functions are then examined through the dressed 
charge formalism at the T = 0. The finite-size correc- 
tions to the groundstate en ergy were computed from the 



tions to the groundstate energy were computed 1 
BA equations, with result ( Lee and Guanl . [201 if ) 



7T 

6I2 



E 



(110) 



where v u and Vb are the velocities of unpaired fermions 
and bound pairs, respectively. 

Three types of low-lying excitations are considered in 
the calculations of finite-size corrections: a Type 1 exci- 
tation is characterized by moving a particle close to the 
right or left Fermi points outside the Fermi sea; a Type 
2 excitation arises from the change in total number of 
unpaired fcrmions or bound pairs; a Type 3 excitation is 
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FIG. 26 (a) Momentum distribution function of pairs in the 
ID attractive Hubbard model with U/t — —8. (b) Power-law 
decay behaviour in real space and comparison of the DMRG 
result (symbols) to the bosonization (lines) . (c) The positions 
of the peaks in the momentum distribution function n£ alr . 
From Feiguin and Heidrich-Meisner (2007). 



caused by moving a particle from the left Fermi point to 
the right Fermi point and vice versa. Such kinds of exci- 
tation are also known as backscattering. All three types 
of excitations can be unified in the following form of the 
finite-size corr ections for the energy and total momentum 
of the system ( Lee and Guanl . l201lh 



AE = T (l *( AiV )'( z_1 ) vz_lAiV 

+ t {AD)ZVZ t AD + + 



a — u , b 



ap = 



2tt 

T 



l ANAD + N U AD U + N b AD b 



(111) 



a— u,b 

with the notation 



of which has a central charge C = 1. The finite-size 
scaling form of the energy (|111[) determines the critical 
exponents of two-point correlation functions between pri- 
mary fields (O' (x,t)0(x' ,t')). At zero temperature, the 
two-point correlation functions take the form 

(O(x,t)O(0,0)) = (113) 

exp(-2i{N u AD u + N b AD b )^x) 

(x - iv u t) 2A t ( x + h; M i) 2A " (x - iv b t) 2A t (x + iv b t) 2A > ' 



The conformal dimensions are given by 

Z bb AN u 



2A 



Z b „ AD b ± 



Z, lb ANi 



2 dot Z 



-2N1 



2A± = [Z ub AD u + Z bb AD b ± 



-2N±. 



Z„,,AN h - Zfr,,.AN,, 



2dctZ 



Here {ot = u, b) characterize the descendent fields 
from the primary fields. 

In this way the quantum numbers for the low-lying ex- 
citations completely determine the nature of the asymp- 
totic behaviour of these correlations. The exponential 
oscillating term in the asymptotic behaviour comes from 
the backscattering process. Various correlation func- 
tions, e.g., the single particle Green's function G^(x,t), 
charge density correlation function G nn (x,t), spin cor- 
relation function G z (x,t) and pair correlation func- 
tion G p (x,t), can be derived based on the choice of 
(AN Ul AN b ) which d efine the quantum numbers (JTT3J) 
( Lee and Guanl . 120 111) . We now discuss these correlation 
functions in detail. 

The asymptotic form of the single particle Green's 
function G-\(x, t) = {ipAx, t)ipf(0, 0)) is given by 



G t (x,t) 



Af^i cos (jr(n^ — 2n\ r )x) 
\x + w u t\ 6 ^\x + \v b t\ i 
A+. 2 cos (im^x) 
\x + iv u t\ e3 \x + iv b t\ e -i ' 
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(Qu) Z ub (Q b ) 
Z bu (Qu) Z bb (Q b ) 



Here the quantum numbers AN,, b are characte r ized b y 
the change in quantum numbers ( Lee and Guanl . 1 2 llh 



AD n 



AN U + AN b 



AN, 



(mod 1), AD b := (mod 1). 

^ 2 

(112) 

The dressed charges Z UU (Q U ), Z ub (Q b ), Z bu (Q u ) and 
Z bb (Q b ) satisfy a s et of dressed charge equations 
dLee and Guanl . 1201 llh 

The finite-size spectrum is described by a critical the- 
ory based on the product of two Virasoro algebras each 



where the critical exponents for the strong coupling 
regime arc given in terms of the polarization by 
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1 - 
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The constants Af i and A^^ cannot be derived from the 
finite-size corrections for low-lying excitations. 

The asymptotic form of the pair correlation function 

G p (x,t) = (tlj\(x,t)^l(x,t)tP t (0,0)^4 (0,0)) is given by 



G p (x , t) 



A pA cos (7r(n t - n±)x) 
\x + iv u t\ 6l \x + iv b t\ e * 
| A Pt2 cos (7r(n T - 3n^)x) 
\x + iv u t\ 6s \x + iv b t\ 84 - 



(115) 
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where the critical exponents in the strong coupling 
regime are 



71 ~ 21 

3 ~ 1 I 2(1-P) 

*J ~ 2 + hT> 



^ ~ 2 + "2M - ' 

3 ~ 5 _ (19P-3) 
; 4 ~ 2 2| 7 | ■ 



It is clearly apparent that the leading order for the 
long distance asymptotics of the pair correlation func- 
tion G p (x,t) oscillates with wave number Akp, where 
Ak F = 7r(nt - 

The Cooper pair correlatio n has been further discussed 
with a connection to CFT ( Schlottmann and Zvvaginl 
2012b). Moreover, the leading order for the charge den- 
sity correlation function G nn (x,t) = (n(x, i)n(0, 0)) and 
the spin correlation function G z {x,t) oscillates twice as 
fast with wave number 2/S.kp 1 namely 



G nn (x, t) r? + 



A nn ,i cos {2n{n^ ~ n^x) 
\x + \v u t\ 9 i 



A nn ,2 cos {2-kN^x) A n7li3 cos (2ir(n t - 2n l )x) 



\x + k'bt\ 
G z (x,t) « (m z Y + 



iM 6 



Z x2 cos (2n(n^ — ni)x) 



x 



iv u t\ 



A Z}2 cos (2tto;x) A^ 3 cos (27r(rt T - 2rt ; )x-) 
|a; + iu h t| e = |x + i^| e 3|x + iu 6 t| e * ( j 

Here A Z: , and A nn are constants and the correlation ex- 
ponents are given by 



4(1-P) 
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2(1-P) , 16P 
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The oscillations in G p (x,t), G nn (x,t) and G z (x,t) 
are caused by the mismatch in Fermi surfaces be- 
tween both species of fermions. These spatial 
oscillations give a novel si gnature of the Larkin- 
Ovch inikov pairing phase ( Larkin and Ovchinnikovl 
Il965h . This finding is consistent with the numerical re- 
sults from DMRG djfeiguin and Heidrich-Meisneii 12007 : 
Liischer et oZ] . l2008tlRizzi et all 120081: iTezuka and Ueda , 



2008U2010D and QMC (IBatrouni et all 120081: iBaur et at , 
20ialwolak et aZj . l201oF 7 It is remarkable to see that the 



asymptotics of the spatial oscillation terms in the pair 
and spin correlations are a consequence of Type 3 exci- 
tations, i.e., backscattcring for bound pairs and unpaired 
fermions. The asymptotic behaviour of the Fermi field, 
Cooper pair and charge density wave correlation func- 
tions for the ID attractive Hubbard model were com- 
puted for t he magnetic field H — > H r ^ + + at the half- 
filled band ([Bogoliubov and Korepinl Il990h . 

Furthermore, the correlation functions in momentum 
space can be given by taking the Fourier transform of 
their counterparts in position space, i.e., the Fourier 
transform of equal-time correlation functions reads 



g(x,t = + ) 



exp(ikox) 



(x--i0) 2A+ (a;-(-i0) 



2A- 



(117) 



where A* = Aj + A^ is given by 

g(k « k ) ~ [sign(fc - k )] 2s \k ~ k \ v . (118) 

Here the conformal spin of the operator is s = A + — A - 
and the exponent v is expressed in terms of the con- 
formal dimensions by v — 2(A + + A~) — 1. However, 
experimental observation of this momentum distribu- 
tion remains very difficult because the transverse expan- 
sion dominates the expansion along the axial direction 
( Bolech et al 1 l2012h . In this regard, it is practicable to 
consider the long-time behaviour of the distribution dur- 
ing the sudden expansion of spin-imbalanced ultracold 
lattice fermions with an attraction after turning off the 
longitudinal confining potential. However, the existence 
of the F FLO signature in the expand i ng ga s is still in 
question (|Bolech et all 120121 : iLu et all l2012h . 

For fields above the lower critical field H c i, these cor- 
relators decrease as power-laws. Thus the pairs lose 
their dominance, i.e., three types of ordering coexist - 
superconductivity, charge density waves and spin den- 
sity waves. Indeed the FFLO correlation is more ro- 
bust in the ID homogenous attractive Fermi gas. In 
the experimental setting discussed in Sec. VII, the quasi- 
1D system of ultracold fermionic atoms is formed by 
loading into a bunch of quasi-lD tubes created by two 
counter-propagating laser beams. By tuning the inten- 
sity of the lasers one could adiabatically tune the tun- 
nelling between two neighbour tubes. Such Josephson- 
like tunnelling is likely to lead to a 3D long-range or- 
der of the superconducting phase. A quasi-lD model 
consisting of at tractive Hubbard chains with in terchain 
tunnelling t± I Bogoliubov and Korepinl Il989l ). where 
t± is much less than the energy gap A. For a sin- 
glet pair state, the charge density wave is suppressed, 
thus the system shows an anomalous average value, i.e., 
(^j-j-^ij.) 7^ 0. This means that a superconductive cur- 
rent exists. Through analysis of the instability of the 
normal state (refer to the TLL phase) to the supercon- 
ductor transition, the critical temperature is estimated to 
be T c ~ A[g(i_ L /A) 2 ] 1 /( 2 ~' )/ ), where the 7' is the critical 
pair correlation exp onent of the ID homogeneous a ttrac- 
tive Hubbard chain dBogoliubov and Korepinl . ll989t) . For 
T > T c , the gapless excitation spectrum of the Cooper 
pair leads to power-law behaviour of pair correlation 
function. 



B. ID two-component repulsive fermions 



In Sec. IIII.FI we saw that for ID repulsive spin- 1/2 
fermions, the low energy physics of the model can be 
reformulated as two massless bosonic theories for the 
charge and spin degrees of freedom with dispersions 
uj(q) = v C:S \q\. Based on this spin-charge separation, the 
bosonization approach has been used to co mpute corre- 
lation functions of the ID Hubb a rd mod el dGiamarchi , 
20041: iRen and Anderson! . 119931: ISchulzl Il990l 11991: 
Tsvelikl . ll995ir Due to the fact that the spin and charge 
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excitations are independent of each other, the descrip- 
tion of critical phenomena in the repulsive Hubbard 
model h as been mad e by the conformal field theory ap- 
proach (lEssler et all l2005t iKawakami and Yantd . |'l99ll 
IWovnarovichl 19891 ) . Using the critical theory based 
on the product of two Virasoro algebras with central 
charge (7=1, one can systematically compute asymp- 
totics of correlation functions from the finite-size spec- 
trum of low-lying exc i tations in te rms of the BA solution 
(|Frahm and Korepinl . Il990l Il99ll) . \ n the same fashion, 
the long distance asymptotics of various correlation func- 
tions of the spin- 1/2 Fermi gas have been in vestigated in 
the strong repulsive regime ( Lee et all . l2012l ). The model 
is gapless and thus critical at zero temperature. At T = 
the correlation functions decay as some power of distance 
governed by the critical exponents. For T > the decay 
is exponential. 

In the context of spin-charge separation, the general 
two-point correlation function for primary fields <p with 
conformal dimensions Af s at T = and T > are given 
by 



(<pOmmo,o)) = 

and 



e -2iD a (k F ^ + k F i)x e ~2iD 3 k F ix 

Ua=c > - *Vat) 2A " (X + iv a t)^ 

(119) 



(<p(x,t)<p(0,0)) T = (120) 

(nT/Va) 2 ^ +A ° ' e ~2iD c (k Ft +k Fi )x e -2iD s k Fi x 



n„ 



sinh 



2A2 



sinh 



2A' 



^(x + iv a t) 



where kp\^ are the Fermi momenta, < x < L and 
— oo < t < oo is Euclidean time. The conformal dimen- 
sions of the fields can be written in terms of the elements 
of the dressed charge matrix as 



/ Z SS AN C - Z CS AN S 

2Af = Z CC D C + Z SC D S ±- SS 



2dctZ 



-2m 



2detZ 



+2A±. 



Here the dressed charge matrix is denoted by 

(Qc) Z cs (Qs) 



z = 



Z SC (Q C ) Z ss (Qc) 



which can be obtained fro m the dressed charge equa- 
tions, see ( Lee et al\ . 120121 ). Q c s are the Fermi bound- 
aries for charge and spin degrees of freedom. The non- 
negative integers A N a , and the parameter D a where 
a = c, s characterize the three types of low-lying excita- 
tions. AN a denotes the change in the number of down- 
spin fermions. characterizes particle-hole excitations 
where iV+ (N~) is the number of particles at the right 



(left) Fermi level jumps to. D a represents fermions that 
are backscattered from one Fermi point to the other, i.e., 



D r = 



AN S + AN S 



(modi), D s = 



AN C 



(modi). 



With the expressions (|119[) and (|120|) the general two- 
point correlation functions for the operator 0(x,t), 
namely (0(x, t)Cft (0, 0)) , can be written as a linear com- 
bination of primary fields with conformal dimensions A f s 
and their descendent fields from the finite-size spectra of 
the model. 

Various correlation functions of operators can be writ- 
ten in terms of the field operators ip a (x,t) and ip a (x,t) 
where a =t> 4-- E.g., 

(i) One particle Green's function: 

G CT (M) = WvOM)^(0,0)). 

(ii) Charge density correlation function: 

G nn (x,t) = (n(x,t)n(0,0)) 

where n(x,t) = n^(x,t) + n±(x,t) and n a (x,t) = 
ipl(x,t)^ a {x,t). 

(iii) Longitudinal spin-spin correlation function: 

G z (x,t) = (S z (x,t)S z (0,0)) 

where S z (x,t) = ^(n^(x,t) — n±(x,t)). 

(iv) Transverse spin-spin correlation function: 

G ± {x,t) = (S + (x,t)S~ (0,0)} (121) 

where S + (x,t) = i/j1(x, t)tpi(x,t), and S~(x,t) = 

il>l(z,t)il>t(x,t). 

(v) Pair correlation function: 

G p {x,t) = (^(x,t)^(x, t)V{(0, 0)V|(0,0)>. (122) 

The critical exponents for each of the above correlation 
functions are determined by the values of the quantum 



state with 






G t (x,t) : 


(AN C — 


1,AA S 


Gi(x,t) : 


{AN C = 


l,AN s 


G nn (x, i) . 


(AA C = 


0,AN s 


G z {x,t) : 


(AA C = 


0,AN s 


G ± (x,t) : 


(AA C = 


0,AN s 


G p (x,t) : 


(AA C = 


2,AA S 



0, D C 6 Z + ±,D S 6 Z- 

1, D c eZ,D s e2 + i) 
0,D C eZ,D s e Z) 



1, D c G Z • 



D, e z) 
D s e Z) 



with N^ s G Z>o for every case. The explicit results for 
these correlation functions for H <C 1 and H 



H r are 



given by solving the dressed charge equations ( Lee et all , 
2012). In the zero field limit H -C 1 the dressed charge 
equations were solved by the Wiener-Hopf method, ex- 
plicitly, 



Zss(Qs) 
Zsc(Qc) 
Zcs(Qs) 
ZcciQc) 



V2 ( 



&n±H 



1 

27x4. In 2 



2 1 c 
2V2H 

2 1 2 In 2 

7 



1 i\n(H /H) 
2H 
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Here Hq = \J ir 3 /(2e)H c with the critical field value H c w 
8n 3 7r 2 /(3c) for strong repulsion. 

For the approach to the critical field h —> H c the 
dressed ch arge equations can be solved by a symptotic 
expansion (jEssler et al ], 120051 iLee et aZ.1 . 120121) . with re- 
sult 



Z SS (Q S ) 

Z S c{Qc) 
ZcsiQs) 
Z C c(Qc) 



1 



1 



H 



8 

77 1 



1 



H 

H c 



1 - 



H 
H € - 



±1-1 



1 - 

1 H € 



1 + 

7T7 



1-^- 

He ' 



A few terms of the asymptotic expansion for the dressed 
charges together with the above values for the low-lying 
excitations determine the asymptotics of the correlation 
functions. As an illustration of this approach, in the 
zero field limit, the asymptotics of the density-density 
correlation function are given by 



G nn (x, t) ~ n 2 + 



Ai cos(2kp^x) 



(123) 



x + iv c t\ e ^\x + ivstf' 1 
A2Cos(2kp^x) A 3 cos(2(kpi + kp^)x) 



iv c t\ 



ivM 



iv c ty 



where the critical exponents are 



7 C 1 
9c3 
%2 



2 In 2 

7 

2 In 2 

7 

8 In 2 



2\n(H a /H) 
1 

2\n(H a /H) 



H " 


\ _ 8 In 2 




Hc < 




7T 2 7 


I He 


H 

H Cj 


) + 


8 In 2 

7T 2 7 


{h c 




' H^ 






u 


H c , 







The presence of the external field does not change the 
form of the correlator (|123j) . In the la rge field limit, i. e., 
H — >• H c , the exponents are given by (|Lee et ad 120121) 



,1 


= 2- 




/l- 


4 


- 32 


C2 


= 2 + 


32 




H 




,1 


= 2- 




/l- 


H c ^ 


32 


,2 


= 2- 


16 

7 






ff 

He 



1 - JL 

1 H c 



H 

H c 



(124) 



64 

7T7 



1 - 



H 

H, 



The constants A; with i = 1,2,3 depend on the in- 
teraction. For kp^ = kp^ = nn/2 = hp, we see 
that the density-density correlation contains 2kp and 
4/cF-oscillations. However, in the strong coupling limit, 
the system behaves like non-interacting spinless free 
fermions, thus the 2kp terms should vanish, i.e., Ai = 
A2 ~ 0. The leading contributions to the asymptotics of 
the density-density correlation function are from the 4k p- 
oscillation term because this oscillation is a consequence 
of int eractions ( Essler et all 120051 : iFrahm and Korepinl . 
Il990h . 



Furthermore, the large distance behaviour of the two- 
point correlation function determines the singularities of 
spectral functions near u sa ±v c ^ s ([k — kp). The correla- 
tion functions in momentum space can be determined by 
Fourier transforming the asymptotics of the above corre- 
lators. Of particular interest are the dynamical response 
functions such as the spectral function A(to, k) which can 
be obtained from the imaginary part of the retarded sin- 
gle particle Green's function. The interacting spectral 
function often has a non-zero width. There are singular- 
i ties in the sp ectral functi on Aju ), kp + q) for u — > v C}S q 
(lEssler et all 120051: lEsslerl . l2010h , with 



A(w, kp + q) 



(lu — v c q)^2~ forw 

(lo + v c q) iovui 

(u> — v s q) ai ~z forw 

(oj + v s q) ai ioruj 



v c q, 

- y c1> (125) 

v s q, 

—Wag- 



Here the exponent ai, which can be obtained from 
Fourier transformation, is always greater than zero. Us- 
ing TLL theory, the Fourier transforms of the zero- 
temperatu re single particle Green's function have bee n 
computed (|Meden and Schonhammerl Il992t IVoitl . Il993l) . 
The Fourier transform of the 2k p Luttinger liquid 
density correlatio n function was calculated recently in 
( Iucci et aZl . l2007t ). In general, the explicit calculation of 
the spectral function is muc h more i nvolved (lEssler et al 



2005 



2005) 



IFrahm and Korepinl 1 19901 : IFrahm and Palacios , 



C. ID multi-component fermions 

In general, SU(k) Wess-Zumino-Witten (WZW) the- 
ory can be used to capture the low energy be- 
haviour of a family of cr itical quantum spin chains 
([Affleck and Haldand . Il987h . The SU{ k) WZW mod- 
els of level £ describe integrable higher spin models with 
critical points characterised by the central charge C = 
£(k 2 ~ 1) / (£+k) and the scaling dim ensions of the primary 
fields (|Affleck and Haldand . Il987t ). The U(l) ® SU(k) 
symmetry interacting fermions in ID have k branches of 
states characterised by one charge degree of freedom and 
K — 1 spin rapidities, see Eq. (|93|) . The low- lying exci- 
tations are described by the linear dispersion relations 
u r (k) = v r (k — kp) with the charge velocity v c and spin 
velocit ies v r for r = 1, . . . , k — 1, respectively. The BA 
result (jGuan et all |2010( ) predicts that the energy has a 
universal finite-size scaling in the low energy excitation 
spectrum, i.e., E 0>L - E a ^ = -f£Er T V + Oi 1 / 1 ^ 2 )- 
The low energy spectrum can be interpreted in terms of 
a product of k independent Virasoro algebras with the 
same central charge C = 1. For vanishing magnetic field, 
the spin velocities are equal. Thus we have a spin-charge 
separation of the C = 1 Gaussian field theory (in the 
charge sector with U(l) symmetry) and C = k— 1 WZW 
theory (in the spin sector with SU(n) symmetry). 

Using the BA solution, Frahm and Schadschncider 
(1993) and Kawakami (1993) have calculated finite-size 
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corrections and the spectrum of the low-lying excitations 
in the ID multi-component degenerate Hubbard model. 
The asymptotic behaviour of various correlation func- 
tions has been determined from these low-lying spec- 
trum. In the same fashion, the asymptotics of corre- 
lation fu nctions for the Bose - Ferm i mixtures have been 
studied ( Frahm and Palaciosl 120051) . The general struc- 



ture of the critical exponents for the multiple nested BA 
solvable models has been well understood, see a review 
by Schlottmann (1997). 

The low-lying excitations above the groundstatc 



have been obtained dFrahm and Schadschneidei 


1993: Kawakami. 19931: 


Schlottmann. 1997 


Schlottmann and Zvvaeinl 2012all 





AE = T (l '^N)'^ 1 ^ 2 " 1 ^ 



+ '(AD)ZVZ'AD + ^v r (N+ + N~) , (126) 

r=0 / 

AP = — 'ANAD + ]T A r AD r + ^(A+ - N~) , 

\ r=0 r=0 ) 

where V = diag(u c ,?;x, • • • ,i> re _i) and iVjr are positive 
integers characterizing the excited states at the branches 
r = 0, 1, ...,«; — 1. Here r = stands for the charge 
degree of freedom. AN and AD are vectors, where AN r 
in vector AN denotes the changes of the total numbers 
in each branche. The values of AD r in vector AD are 
integer or half-odd integer depending on AN r , i.e., 



AA e + AN, 
AD C := (mod 1), 

AN r - 1 + AJV r+ i 
AD r :~ (mod 1) (127) 



with r = 1, ...,K- 1 and AA = AA C , AN K = 0. 

The conformal dimensions A^r of the primary field 
are given in terms of the dress ed charge matrix D 
( Frahm an d Scha dschneiderl Il993l ) 



2A± = ( (Z*AD) r ± ^(Z-^NJr ) + 2A±. 



The dressed charges 
obtained by solving 
tions obtained from 



Zij 
the 
the 



Zij(Qj) can be 



dressed charge equa- 
dressed energy equa- 
tions by definition (IBogoliubov and Korepinl Il990t 
iFrahm and Schadschneidei , Il993l ~ lzergin el all Il989t) . 
Consequently, different dressed charge equations arc 
accordingly derived for the multi-component Fermi gases 
in the repulsive and attractive regimes. The dressed 
charges for the degenerate SU(k) Hubbard model was 
calculated explicitly in the absen ce of magnetic field 
( Frahm and Schadschneiderl Il993l) . The single particle 
Green's function and charge density-density correlation 
were thus studied. In general, using the CFT result, the 



asymptotics of correlators for primary fields at T 
are given by 



(<p(x,t)<p(0,0)) = 



exp 



-2i 



n: 



(129) 

Here k r F = Trn r is the Fermi momentum of each 
branch. The critical exponents of the correlation func- 
tions of ID multi-component Fermi gases at zero temper- 
ature in external magnetic fields can be calculated in a 
straightforward way. Some response functions in multi- 
component Luttinger liqui ds were computed analytically 
( Orignac and Citrol . 120121 ). The quantum numbers of 



low-lying excitations, such as AN r ,AD r and N^, com- 
pletely determine the nature of the asymptotic behaviour 
of the correlation functions. There are many superfluid 
orders in multi-component attractive Fermi gases, for 
which the operators for the superfluidity (p(x) are writ- 
ten in terms of fermion anniliation operators of length- 
r. The instability of superfluidity and normal phase can 
be determined by the correlation function (|129p . In the 
context of ultracold gases with higher spin symmetries, 
Schlottmann and Zvyagin (2012a, 2012b) have studied 
various superfluidity ordering in spin-3/2 and spin-5/2 
attractive Fermi gases, in which the relevant superflu- 
idity operators of few-body bound states were labeled. 
However, the study of colour superfluidity, the FFLO-likc 
modulated ordering in ID interacting fermions with large 
spins still remains preliminary. It remains to further in- 
vestigate magnetism, superfluidity, charge and spin den- 
sity waves by using CFT and the BA solutions . 



D. Universal contact in ID 

The universal nature and phenomena of interacting 
fermions, such as Landau's Fermi liquid theory, Luttinger 
liquid theory, and quantum criticality, have always at- 
tracted great attention from theory and experiment. Tan 
(2008a, 2008b, 2008c) has shown that a few universal rela- 
tions for two-component interacting fermions involve an 
extensive quantity called the universal contact C. The 
first Tan relation is for the tails of the momentum distri- 
bution which exhibits a universal n a (k) ~ C/fc 4 decay as 
the momentum tends to infinity. Here we denote a ^. 
The constant C measures the probability of two fermions 
with opposite spin at the same position. Secondly, this 
contact also reflects the rate of change of the energy or 
free energy due to a small change in the inverse scattering 
length 1 1 a for fixed entropy s or temperature T 



dE 



Airm 



-C. 



dF 
da~ l 



Airm 



C. (130) 



Tan's additional relation states that the pressure and the 



¥ 



-c. 



energy density are related by P - 3 ^ , 12mna - 

The Tan relations were found as a consequence of 
operator identities following from a Wilson operator 
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product expansion of the one-particle density matrix 
(Br aaten and Platted . l2008h . These relations hold more 
generally as long as the interaction range ro is much 
smaller than any other characteristic length scale like 
the ave rage interparticle distance and the thermal wave- 
length ( Zhang and Leggettl . 120091 ). Prior to Tan's results, 
the derivative of the energy E with respect to the in- 
verse scattering length a has been experimentally exam- 
ined from measurements of the photoassociation r ate of 
a trapped gas of 6 Li atoms ( Partridge et al\ . [200511 . The 
static and dynamic structure factor for the trapped gas 
of 6 Li ato ms has been studied the by using Bragg spec- 
troscopy ( Kuhnle et all 120101) . In particular, Tan's uni- 
versal contact and thermodynamic relati ons have been 
confi rmed for a trapped gas of 40 K atoms ( Stewart et ali 
I2010T ) . Recent developments on the Tan relations for 
fermions with large scattering length have been reviewed 
by Braaten (2012). 

The generalization of Tan relations to ID interact- 
ing fermions has been studied by Barth and Zwerger 
(2011) using the Operator Product Expansion method 
developed by Kadanoff (1969) and Wilson (1969). The 
Tan adiabatic theorem has been generalized to the ID 
Gaudin-Yang model for which the ID analog of the 
Tan adiabatic theorem is given in the form -r^E = 

° aa 1D 

Tan's universal contact also exists in the asymp- 
totic behaviour of the tail momentum distribution of 
the ID Fermi gas. In terms of the fields tp^(R) and 
?/v(i?) the momentum distribution given by h a (_k) = 
JdR jdx e- ikx (ipi(R)ipa(R + x)}. Using the Operator 
Product Expansion method, Barth and Zwerger (2011) 
found that the momentum distribution of the ID two- 
component Fermi gas behaves as n a — > C/fc 4 as k — > oo. 
Furthermore, using scaling analysis, they also derived 
the pressure relation which connects pressure and energy 



C(ai) 
4m 



for 



density £ via the universal contact p = 2£ + a\o 
the ID Fermi gas. 

Tan's universal contact also arises in the additional re- 
lation for the pair distribution function and the related 
static structure factor. From the short distance expan- 
sion of the two particle density matrix Barth and Zwerger 
(2011) found the total pair distribution function at short 
distance 



(131) 



This short distance singularity gives rise to a 1 /q 2 power 
law tail in the associated static structure factor S(q) = 
1 + nJ dxe-' iqx (g i2) (x) - 1), i.e., 



S(q 



) = l-4n t n iig $(0)/q' (132) 



with the dimensionless interaction constant 7 = 
—2/{naxo). In fact, for the ID two-component Fermi 
gas, the universal contact is obtained by calculating the 
change of the interaction energy with respect to inter- 
action strength by the Hcllman-Fcynman theorem, i.e., 



C 



4 ( 2 )/n\ 



(2) 

The local pair correlation j is 



accessible via the exact BA solution through the relation 
-dE/dc. Here c = — 2/ am and E is the 



groundstatc energy per unit length. 

In a similar way, for the ID K-component Fermi gas, 
there exists a ID analog of the Tan adiabatic theorem 
where the universal contact is given by the local pair cor- 
relations for two fermions with different spin states. The 
two-body local pair correlation function is given by the 
expectation value of the four-operator term in the second 
quantized Hamiltonian, explicitly, g^ a , (0) = / w _*) n a ^~ 



with k > 1 (|Guan et q/J . l2012D . E is again the ground- 
state energy per unit length. From the asymptotic ex- 
pansion result for the groundstate energy of the bal- 
anced Fermi gas obtained in Section |y.C.2l we easily find 
(2) 

g aa ,(0) — > 1 as |c| — > 0. The local pair correlation can 
be obtained from the groundstate energy (|95p of the bal- 
anced gas with strong repulsion, 



<e«>) 



4k7T 



3(k- 1)7 2 



6Z 2 24 / . Z 3 n 2 
7 7 2 V 1 15 



(133) 



This local pa ir correlation reduces to tha t for the spin- 
less Bose gas ( Guan and Batchelorl 1201 lh as k — > 00. 
For strong attractive interaction, the local pair correla- 
tion fo^_two_f^mions with different spin states is given 
by dGuan et alJ . l2012ft 



,(2) 



(0) 



k{k + 1)\j\ 



4tt 2 



3k 5 (k — 1)7 2 



A K 



6AI 

K 2 M 



24 f , 3 B K 
k 4 7 2 V K 15 



+ 0(l/ 7 4 ). (134) 



Here A K = Y,r=i I/ 7 ' an( l B * = J2r=i V 7 " 3 - ^ is noted 
that for the whole interaction regime, the local pair cor- 
relations for the balanced K-component Fermi gas tend 
to the limiting value for the spinless bosons as K — > 00. 
Fig. [27] shows the local pair correlation for two fermions 
with different spin states in the multi-component Fermi 
gas. 



VII. EXPERIMENTAL PROGRESS 

In order to realize ID systems in experiments, one has 
to apply strong confinement in two transverse directions 
and allow free motion along the longitudinal direction. 
ID systems have been simulated in various settings from 
solid electronic materials to ultracold atomic gases. Al- 
though some ID features have been observed in solid 
electronic materials, it is still challenging to control ma- 
terial parameters and separate imperfect effects caused 
by defects and Coulomb interactions etc. In contrast, ul- 
tracold atomic gases trapped in external potentials have 
high controllability and clean environment. In this sec- 
tion, we briefly review the experimental developments in 
confining quantum systems of ultracold atoms in ID. 



40 



3.0 



2.0 



1.0 



0.5 



» , , , , . , 




\ \ 




- \ \ 










\ \ 








- \\ \ 




— Bosons 




\\ I 
















\\ \ 
- ^\ 1 




— K=10 




\\ 1 
















: \\ \ 




K=4 




















- N \ 








N \t\ 




K=2 
























\ i 




\ \ 
\ \ 




\ 





FIG. 27 The local pair correlation g^ a ,(0) vs 7 = —2/ (nam) 
for the balanced multi-component Fermi gas for k — 2,4, 10 
and for the spinless Bose gas. The solid lines are the numerical 
solutions obtained from the two sets of Fredholm equations. 
The dashed lines are the asymptotic limits obtained from the 
first term in (fT34)l . From Guan et al. (2012). 



A. Realization of ID quantum atomic gases 

Neutral atoms can be trapped by coupling their per- 
manent or induced dipole moments to electromagnetic 
field gradients. By using a laser field (or an inhomoge- 
nous magnetic field), it is possible to trap neutral atoms 
via the interaction between the induced electric dipole 
moment and the laser electric field (or via the interaction 
between the permanent magnetic dipole moment and the 
magnetic field). Correspondingly, there are two typical 
techniques for realizing ID quantum atomic gases: opti- 
cal lattices and atom chips. 



1. Optical lattices 

Optical lattices arc created by superimposing one or 
more pairs of counter- propagating laser beams ( Blochl . 
2005), see Fig. EHl The laser-atom interaction results in 
an ac-stark shift dependent on the laser intensity and 
the detuning from resonance 5 = u>l ~ w„, in which uil 
is the laser frequency and u) a is the atomic transition 
frequency. Therefore, due to the reason that optical lat- 
tices have periodic laser intensities, their ac-stark shifts 
provide periodic potentials for confining atoms. For a 
negative detuning (red detuning), atoms are attracted to 
the region of large laser intensities. On the contrary, for 
a positive detuning (blue detuning) , atoms are attracted 
to the region of small laser intensity. 

Similar to the electrons in a conductor, even if the lat- 
tice well depth exceeds the kinetic energy of the atoms, 
the atoms confined in optical lattices can move among 
neighbouring lattice sites via quantum tunnelling. In 
the case of Bose atoms, the quantum phase transi- 
tion between supcrfluid and Mott insulator phas es has 
been experimentally observed ( Greiner et aLl .[2~002j). The 




FIG. 28 Optical lattices of different spatial configurations, 
(a) For a 2D optical lattice formed by superimposing two 
orthogonal standing waves, the atoms are confined to an array 
of ID potential tubes, (b) For a 3D optical lattice formed by 
superimposing three orthogonal standing waves, the atoms 
are approximately confined by a 3D simple cubic array of 
harmonic oscillator potentials at each lattice site. From Bloch 
(2005). 



groundstate appears as a supcrfluid when the system 
is dominated by the quantum tunnelling among neigh- 
bouring sites. Whereas the groundstate appears as a 
Mott insulator if the system is dominated by the on- 
site interaction between atoms. For fermionic atoms, the 
metal-insulator transition and Neel anti-ferromagnetic 
states etc have been demonstrated in recent experi- 
ments ( Jordcns et al 1. 120081 : ISchneider et al 1. l2008h . 

Ultracold atomic systems confined in optical lattices 
have highly controllable parameters such as the inter-site 
hopping strength, on-site interaction strength and dimen- 
sionality. The inter-site hopping strength can be adjusted 
by tuning the laser intensity of the optical lattices. The 
on-site interaction strength can be tuned from positive in- 
finity to negative infinity by using Feshbach resonances. 
In particular, the dimensionality can be tuned from 3D to 
ID by controlling the spatial configuration of the optical 
lattices, see Fig. [28l 

The ID lattice can be created by a pair of counter- 
propagating laser beams to form of a standing-wave, 



Vldl{x) = V sin 2 (k L x) 



(135) 



The 2D square lattice can be formed by two orthogonal 
standing waves, 

V 2 dl {x,y) = V Q [sin 2 (k L x) + sin 2 (k L y)} . (136) 

The 3D simple cubic lattice can be produced by three 
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orthogonal standing waves, 

V 3DL {x, y, z) = V [sin 2 (k L x) + sin 2 (k L y) + sin 2 (k L z)] . 

(137) 

By imposing laser beams under different spatial configu- 
rations, it is also possible to make more complex lattices 
such as the triangle, honeycomb and kagomc lattices. 

It has been demonstrated that an array of ID systems 
can be created by superposing a harmonic potential onto 
a 2D optical lattice. The potential for such a system is 
of the form 



V(x,y,z) = V 2D L(x,y) + V har {x,y,z). 



(138) 



Under the tight-binding condition, each lattice well can 
be regarded as an independently harmonic potential. 
That is, the potential (x, y, z) around the lattice site 
[xi , yj ) can be expressed approximately as 



Vij(x,y,z) 



1 2 

2 mUJ xy 



(x-Xi) +{y-yj)' 



1 2 2 

—mui^z . 
2 



(139) 

If the lattice depth is sufficiently large, we have Lu xy lj z 
and so all atoms will always stay in the lowest transverse 
vibrational state along the x and y directions. The ab- 
sence of transverse excitations means that each lattice 
tube along the z-axis can be viewed as a quasi-lD sys- 
tem. Experimentally, the longitudinal frequency uj z is 
around 2n x 10 ~ 200 Hz and t he transverse frequency 

40 kHz (|Blochl . l2005l ). 



uj X y is around 2tt x 10 



2. Atom chips 

Atom chips arc a kind of nanofabricated atom-optical 
circuit used to trap and manipulate neutral atoms. The 
basic idea of atom chips is to build atomic traps through 
a bias magnetic field and a system of electric wires fab- 
ricated on chip surfaces. Thus the magnetic potential of 
a straight wire has a hole along the wire. To close the 
magnetic potential with cndcaps, there are two simple 
schemes: one scheme is achieved by combining a straight 
wire with an inhomogeneous bias field, the other scheme 
is achieved by combining a bent wire with a homoge- 
neous bias field, see Fig. [29] The basic property, prac- 
tical design and experimental proc edure of atom chips 
have been reviewed in the l i terature ( Folman et jZl. l2002l: 
iFortagh and Zimmerman 3 120071: iReichelL 12002V Here, 
we briefly review how to use atom chips to realize ID 
traps with neutral atoms. 

One can obtain a highly elongated trap if the longitudi- 
nal confinement is very weak compared to the transverse 
confinements. As illustrated in Fig. [29l the wire with 
current I z and the homogenous bias field -Bbias form a 
waveguide along the longitudinal direction. This waveg- 
uide is closed by two perpendicular wires with currents 
ii and I3. Near the trap center, the magnetic field is 
approximated by 



1 



B = B a + -(3(z-z a y + -[—-^)p 



ft 



(140) 




FIG. 29 ID magnetic potential generated by an atom chip. 
The transverse confinements are formed by the bias mag- 
netic field Bbias and the chip wire 7 Z . The longitudi- 
nal confinement is generated by the other two chip wires 
Ii and I3. From the PhD thesis of A. van Amerongcn 
(University of Amsterdam), which can be downloaded from 
http://staff.science.uva.nl/~walraven/walraven/Theses.htm 



with Bo denoting the value of the magnetic field in the 
trap center. The corresponding magnetic potential is 
given as 

V(p, z) = p m B + \^mul (z - z f + -mwi/) 2 (141) 



with the trapping frequencies 
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where p m is the atomic magnetic moment and m is the 
single atomic mass. At finite temperature, a real sys- 
tem enters its ID regime when both the residual ther- 
mal energy (of order ksT) and the chemical potential fi 
are far less than the transverse excitation energy hujj_, 
i.e., {ksT,fi} <C hujj_. Experimentally, the atom chips 
with a Z-shaped wire can successfully generate quasi-lD 
potential tra ps of high frequency ra t io oj \ /u> \\ up to a 
few hundred dBouchoule et all 20111: Hofferberth et al. 



20071: I Jo et all l2007al iReichel and Thvwissenl . 12004; 
Trebbia et all 120061: Ivan Amerongen et aZ.L l2008h . Typi- 
cally, the transverse frequency wj_/(27r) is in an order of 
kHz and the longitudinal frequency W||/(27r) is about 10 
Hz. It has also been demonstrated that a ID-box trap, 
which has nearly constant potential at the trap mini- 
mum combined with tight harmonic confinement in the 
transverse directions, can be f ormed by positionin g two 
wiggles in a long straight wire (|van Es et al I I2010D . 
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B. Tuning interaction via Feshbach resonance 



C. Data extraction 



The ID quantum gases of ^-function contact in- 
teraction U\d{z) = gmS(z) are characterized by 
the dimensionless parameter t dBergeman et all 120031 : 
iDuniko et all . l200ll: lOlshaniil . ll998HPetrov et aU^OOoT . 
which is defined as the ratio between the interaction en- 
ergy e int and the kinetic energy e kin , 7 = = 
. Here, gm denotes the coupling constant and ri\u is 
the ID number density. The coupling constant guj is 
determined by the 3D scattering length a and the trans- 
verse width of the wave function ( Bcrgeman et al ]. l2003t 
l01shaniiL[l99"8ri . 



There are two usual probes for exploring the physics of 
quantum atomic gases. One probe is the correlation func- 
tion, which measures correlations between the atomic 
fields at different positions and times. The correlation 
function has been extensively used for exploring many- 
body coherence and distinguishing order and disorder. 
The other probe is the dynamical structure factor, which 
describes the total probability to populate any excited 
state by transferring both momentum and energy from 
an initially equilibrium state. The dynamical structure 
factor is a powerful quantity for characterizing the low- 
energy excitations. 



9lD 
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where t he constan t A = 1.0326 and the transverse width 
l± = \J hj (mujj_). If \a\ <C l±, the atom-atom scattering 
acquires a 3D character, the interaction strength gm is 
given as g 1D = ^-jr- Theoretically speaking, by vary- 
ing the interaction strength guy, the gradual transitions 
from quasi BEC and the TG gas appear in ID Bose sys- 
tems, the BCS-BEC-like crossover takes place in ID two- 
component Fermi systems. Experimentally, to explore 
different regimes of ID quantum atomic gases, one may 
tune the 3D scattering length a via the well-developed 
techniques of Feshbach resonance. 

Feshbach resonance takes place when the energy of a 
bound state (a state belonging to a closed channel) for 
the inter-particle potential is close to the kinetic energy 
of the two colliding particles. The scattering length of the 
two colliding ultracold atoms undergo a Feshbach reso- 
nance if the energy of a molecular state is close to the 
kinetic energy of the the colliding pair of atoms. The 
interaction strength between a pair of ultracold atoms is 
proportional to the scattering length. A positive (nega- 
tive) scattering length gives repulsive (attractive) inter- 
action. The scattering length a near a magnetic Feshbach 
resonance point Bq is given by 



a(B) = dbg I 1 - 



A 



B-Bc 
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with the background scattering length a^g , the magnetic 
field B and the resonance width A. Obviously, the sign 
and strength of the scattering length can be tuned by ad- 
justing the magnetic field B. Feshbach resonances pro- 
vide an excellent tool for controlling the ultracold inter- 
action between atoms and have been widely used for ex- 
ploring different regimes of quantum atomic gases. The 
theory of Feshbach res onances in ultracold ato mic gases 
has been introduced in ( Duine and StooJ . 120041) . The ex- 
perimental developm ents in this field can be found in a 
recent review paper ( Chin et all l2010l ). 



1. Detecting correlation functions via optical imaging 

The first-order correlation function (also called the 
single-particle correlation function) is 



= (*t( ri )^-(r 2 )). 



(146) 



Here ^(r) and Vt'j(r) denote the field operators for atoms 
in states \i) and Obviously, the first-order corre- 

lation function becomes the density n.i(r) = (pi(r)) = 
(tyf (r)i£ i(r)) if i = j and r = rj. = r 2 . The density 
profile can be mapped out by the well-developed tech- 
niques of optical imaging, in which information about 
an atomic gas is encoded onto the absorption, phase 
shift or polarization information of the probe laser. One 
popular optical imaging technique is absorption imaging, 
which is implemented by comparing imagines taken with 
and without the atomic gas within the field of view and 
recording the fractional absorption of the probe laser. 
However, absorption imaging shows large uncertainty in 
the measured column density in probing high density 
atomic gases. Therefore, absorption imaging is usually 
taken after the atomic gas is released from the trap 
and expanded for a certain time. This imaging method 
is known as time-of- flight imaging ( Altman et all [20041: 
iKetterle etdl , Il999h . 

In time-of-flight imaging, the fractional absorption is 
proportional to the expectation value of the atomic den- 
sity operator m(r, t) = (pi{r,i)) = (*+(r, t)*,(r, t)). By 
assuming negligible atomic collision effects in the ballistic 
expansion during the time of flight, the density rii{v, t) is 
proportional to the initial momentum distribution rii(k) 
with the corresponding wave vector k = mv/t, i.e., the 
time-of-flight imaging gives the initial momentum distri- 
bution 



i(k) = (ttt(k)*i(k)> <xm(r,t). 



(147) 



However, if the time-of-flight is not long enough and 
the atomic collision effects cannot be ignored during the 
time of flight, the measured density distribution n,(r, t) 
will not be proportional to the in itial momentum dis- 
tribution n.j(k) with k = mr/t (jGerbier et al ., 2008: 
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iPedri et all l200l[) . Recently, different from free expan- 
sion via time-of- flight, a controlled ex pansion via short- 
cuts to adiabaticity has been p roposed ( del Cam"pcll201ll 
Idel Campo and Boshierl . 120121) . 

By analysing the correlation of different timc- 
of-flight images, one may reconstruct the density- 
density correlation function (?i i (r 1 )n :; (r2)) = 
( (r ! ) * t (r 2 ) * j (r 2 ) ^j(ri)). This t echnique is known 
as no i se spectroscop y (lAltman et all 12004 IChuu et all 
120051: iFolling et ail , 120051: iGreiner et all l2005al) . In 
different runs of experiments, the time-of-flight images 
have technical noises and quantum fluctuations at the 
same time. Each time-of-flight imaging is a measure- 
ment of the density operator. To explore the quantum 
correlation of the density operator, all technical noises 
should be reduced below the quantum fluctuations. 
Therefore, the measurement times M of imaging must 
be sufficiently large so that the statistical error 1 / \fM is 
small enough. The density-density correlation function 
is a special case of the second-order correlation function 

Gg(n,r 2 ,r3,r 4 ) = (^(n)*/^)*,-^)**^)), 

(148) 

which is also called the two-particle correlation function. 

To reconstruct the full correlation functions, in ad- 
dition to the measurement of density (diagonal corre- 
lations), one has to measure the non-diagonal correla- 
tions. It has been suggested that the non-diagonal corre- 
lations ( ^(k')^(k)) can be measured by atomic inter- 
ferometrv (iPolkqvnikov et al. 1. 120061: ISteneer et q*J . ll999t 
iTorii et aUl2000j) . The Fourier sampling of time-of-flight 
imagines may also be used to reconstruct the fu ll one- 
parti cle and two-particle correlation functions ( Duanl . 
l2006h . In this detection scheme, two consecutive Raman 
pulses at the beginning of the free expansion are used 
to induce a tunable momentum difference to the corre- 
lation terms so that both diagonal and nondiagonal cor- 
relations in the momentum space can be detected. The 
first Raman operation is implemented by two travelling- 
wave beams of different wave vectors ki and k 2 . The 
corresponding effective Raman Rabi frequency has a spa- 
tial dependent phase with fi(r) = floe l (* k-r+Vl ', where 
<5k = k 2 — ki and ipi is a constant phase. The sec- 
ond Raman operation is implemented by two counter- 
propagating laser beams of the resonant frequency for 
the transition between the two involved hypcrfine levels. 

In contrast to the first Raman operation, the ef- 
fective Raman Rabi frequency for the second Raman 
operation is a spatial independent constant Qoe tlfi2 . 
Therefore, through time-of-flight imaging, the real and 
imaginary parts of the nondiagonal correlation function 
(^ r +(k)^ , a (k — (5k)} can be measured by choosing the rel- 
ative phase dip = ip 2 ~<y2i = and 7r/2, respectively. Com- 
bined with the techniques of noise spectroscopy, through 
analyzing the correlations of imagines corresponding 
to (*+(k)* a (k-ak)) and (*+(k')* Q (k'-5kO), 
respectively, the two-particle correlation function 
(*+(k)* Q (k-(5k)*+(k')*a(k'-(5k')) can be re- 



constructed. If the atomic gases have two spin 
components ct\ and a 2 , the spin-spatial correlations 
(\E'+ 1 (ki)\E' a2 (k 2 )) can be reconstructed by a combina- 
tion of Fourier sampling with a pair of Raman pulses 
which mixes the two spin components. Therefore, in 
addition to the single-spin density-density correlation 
functions (*+(k)* Q (k - <5k)tf+(k')*a(k' - 5k')), the 
opposite-spin density-density correlation functions 
<*+ (ki)* Ql (ki -<Jki)*+ (k 2 )* Q2 (k 2 -<5k 2 )) can be 
reconstructed. 



2. Detecting the dynamical structure factor via Bragg 
scattering 

Based upon the backward-scattering of incident parti- 
cles (such as photons, electrons and atoms) from a tar- 
get sample, Bragg scattering is a powerful tool for deter- 
mining both the momentum and the energy absorbed by 
the target sample. In experiments with ultracold atomic 
gases, unlike photons arc diffracted by an atom grating, 
atoms are diffracted on a grating of coherent photons 
(lasers). Bragg scattering provides an effective tool of 
spectroscopy, called Bragg spectroscopy, which can be 
used to detect the dynamical structure factor and explore 
the intrinsic mechanism of the low-energy excitations. 

The dynamical structure factor S(q, uS) describes 
the total probability to populate an excited state by 
transferring a momentum hq and energy Hu>. At 
no nzero temperatures, the dyna mical structure factor 
is (jPitaevskii and Stringar il 120031) 



S(q,w) 



r ^e-^|(0 / |^ + (q + k)^(k)|^)| 2 

x5(E f (q + k) -Ei(k) - Hlj) , (149) 

where \4>i) and \4>f) are initial and final states of the 
many-body system, with energies Ei and Ef. The op- 
erator t/j + (k) creates a particle of momentum k. The 
function Z = £\ e~ l3E ' with /? = l/(k B T) is the parti- 
tion function. At zero temperature, the equilibrium sys- 
tem can only occupy its groundstate and the dynamical 
structure factor is expressed as 



E 



(0/l^ + (q- 
xd(E f (q + k) 



kty(k)|GS>|" 

- E GS Qs) - fiu) .(150) 

In a two-photon Bragg transition, the momentum shift 
q = k 2 — ki is controlled by wave vectors of the two Bragg 
lasers and the frequency cj = w 2 — uj\ is determined by 
the frequency difference. 

The techniques of Bragg spectroscopy have been 
successfully applied to measure the dyna mical struc- 
ture factor of Bo s e con densed atoms ( Ozeri et all 
120051 : IStenger et all 19991) and strongl y interacting 3D 
Bose ( Papp et all l2008h and Fermi ( Veeravalli et all 
I2008D atomic gases near Fcshbach resonance. The tech- 
niques of Bragg spectroscopy have been employed to de- 
tect the elementary excitations in an array of ID Bose 
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gases ( Clement et aZ.I . [2009T ) . By varying the detuning of 
the Bragg lasers, it is possible to probe both the spin and 
density dynami cal structure factors for multi-component 
Fermi systems ( Hoinka et all 120121 ) . 



TABLE I Key experiments in ID quantum atomic gases. 



D. Experiments with ID quantum atomic gases 

Over the last decade, there have been many break- 
through experiments with ID quantum atomic gases. 
These experiments are implemented by using an ensem- 
ble of ultracold Bose or Fermi alkaline atoms occupying 
one or multiple hyperfme levels. The key experiments in 
ID quantum atomic gases are listed in Table HJ 



1. Bose gases 

The ID Licb-Liniger Bose gas is a prototypical many- 
body system featuring rich many-body physics. This 
gas features three different regimes in quasi-lD traps: 
a true condensate regime, a qu asi-condensate regim e and 
the Tonks- Girardeau regime (jPetrov et al I 1200011 . In 
the limit of weak interaction, 7 = mgiu/(h 2 ni]j) <C 
1, i.e., in the quasi-condensate regime, where Bosc- 
Einstein condensation may take place and mean-field 
theory well describes the low-energy physics. In this 
regime, due to the intrinsic nonlincarity from s-wave 
scattering between atoms, matter- wave solitons have 
been observed in sever a l exp er iments (|Becker et al 



20081; iKhavkovich et all |2002| ; IStellmer et all I200S ; 



Strecker et al. , 20021 ). In the limit of strongly repul- 
sive interaction, 7 = mgiu/(h 2 niu) 3> 1, i.e., in the 
Tonks-Girardeau regime, fermionization of Bose atoms 
occurs. By loading ultracold Bose atoms into a 2D op- 
tical lattices, the effective mass of quasi-particles can 
be increased by applying an additi onal periodic poten- 
tial along the longitudinal direction (IParedes et al ]. l2004l) 
and therefore an array of Tonks-Girardeau gases of quasi- 
particles can be prepared by increasing the effective mass. 
The dimensionless parameter 7 can also be changed with- 
out modulating the longitudinal trapping enabling the 
realization of a set of parallel ID Bose atomi c gase s 
in the Tonks-Girardeau regime ( Kinoshita et all 12004) . 
In the quasi-condensate regime of intermediate values 
of 7, the experimental measurements suggest that both 
mean-field theory and fermioniz ation fail to describe this 
crossover (jTrebbia et all 120061 ). The in situ measure- 
ments of the linear density of a nearly ID trapped Bose 
gas on an atom chip indicates a good agreement with the- 
oretical pr ediction from the Yang- Yang thermodynamics 
equations ( van Amerongen et all I 
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Breathing modes and dipole modes are two usual col- 
lective excitations in trapped quantum atomic gases. 
The breathing modes can be excited by time-periodically 
modulating the longitudinal harmonic potential. The 
dipole modes can be excited by suddenly displacing 
the longitudinal harmonic potential. The ratio of 



Group (Leader) 


Research topics 


Amsterdam 


Yang- Yang thermodynamics 


(van Druten) 


non-equilibrium spin dynamics 


Cambridge (Kohl) 


quantum transport 


CNRS 


density fluctuations 


(Armijo, Bouchoule) 


phonon fluctuations 




1D-3D crossover 




three-body correlations 




mean-field breakdown 


ENS (Salomon) 


matter-wave solitons 


ETH 


confinement induced molecules 


(Esslinger, Kohl) 


p-wave Feshbach resonance 




1D-3D crossover 




Bragg spectroscopy 




collective oscillations 


Hamburg (Sengstock) matter-wave solitons 


Innsbruck (Nagerl) 


super- Tonks-Girardeau gases 




sine-Gordon quantum phase transition 




confinement induced resonance 




three-body correlations 


Kaiserslautern (Ott) 


spatiotemporal fermionization 


LENS (Inguscio) 


Bragg spectroscopy 




low-energy excitations 




dynamics of impurities 


Mainz/MPQ (Bloch) 


relaxation dynamics 




squeezed Luttinger liquids 




Tonks-Girardeau gases 


MIT (Ketterle) 


atomic interferometry 




phase fluctuations 




number squeezing 


NIST 


dipole oscillations 


(Phillips, Porto) 


three-body recombination 


Pennsylvania (Weiss) 


Tonks-Girardeau gases 




local pair correlations 




quantum Newton's cradle 


Rice Uni. (Hulet) 


spin-imbalanced Fermi gases 




matter-wave solitons 



Vienna quantum correlations 

(Schmiedmayer) twin-atom beams 

atomic interferometry 
quantum and thermal noises 
non-equilibrium coherence dynamics 



the frequencies of the lowest breathing modes and the 
dipole modes lub/^d has been measur ed in an ar- 
ray o f ID gases in 2D optical lattices ( Moritz et all 
2003). The experimental data show wg/uo ~ 3.1 
for a Lieb-Liniger gas and ljb/^d — 4 for a ther- 
mal gas, which are consistent with the theoretical 
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results (jPedri and Santosl . 120031) . The experimen- 
tal observation of strongly dam ped dipole oscillations 
due to quantum fluctuations (iGea-Banacloche et all . 
120061 : iPolkovnikov and Wand . 12004ft has been re- 
ported (|Fertig et aLl . l2005h . 

The 1D-3D crossover and quantum phase transitions 
from a ID superfluid to a Mott insulator have been ob- 
served in quasi- ID systems of an addi tional lattice poten- 
tial along th e longitudinal direction dFabbri et all . l2012t 
lHaller et all . l2010al IStoferle et all . l2004|). For weak in- 
teraction, the system is still a superfluid at finite lattice 
depth and the transition to the Mott insulator is induced 
by increasing the lattice depth. For strong interaction, an 
arbitrary perturbation by a lattice potential may induce 
a sine-Gordon quantum phase transition from a super- 
fluid Luttinger liquid to a Mott insulator. These quan- 
tum phases can be well distinguished by detecting their 
low-energ y excitations via th e technique of Bragg spec- 
troscopy ((Haller et all . l2010aft . 

The techniques of atomic interferometry have been 
widely used to explore the phase co herence and fiuctu 
ations between diff erent ID ga ses (|Hofferberth et al. 



20081: IJo et all. l2007al lbllcl: iKriiger et all I201C : 



Schumm et all . 120051 ). The second-order correlation 
functions, which relate to the density fluctuations and 
spatial correlations, have been pro bed by in sit u mea- 
surements of density fluctuations (lEsteve"eF Ml l2006t 
IJacqmin et aUl201lHPerrin et al. L l2012h and ex situ mea- 



surements of photo-association rates ([Kinoshita et all . 
12005ft . These experiments confirm that the local second- 
order correlation functions are about 2, 1 and for 
an ideal ID Bose gas, a quasi-condensed ID Bose gas 
and a Tonks-Girardeau gas, respectively. In addition 
to the density correlations, phase correl ations have 
been probed by matter wave interferometry ( Betz et ~all , 
1201 ll ). Moreover, the three-body correlation functions 
of ID Bose gases have been obta ined by measuring 
the three-body recombinati on rate ( Haller et all . 1201 it 
lLaburthe Tolra et al. , 2004) and the in situ measure- 
ments of third-order number fluctuations ( Armiio et all . 
120101 ). In addition to the measurements of local 
correlation functions ( Kinoshita et all . 120051 ), the 
temporal two-body correlation function has been 
measured by the techniques of scanning electron 



microscopy (jGuarrera et al l l2012t ). 

Beyond investigating equilibrium behaviour, there 
have several experimental studies on non-equilibrium 
behaviour in ID Bose gases. The coherence dynam- 
ics in both isolated and coupled ID Bose gases have 
been explored by using an atom chip ( Hofferberth et all . 
I2007D . The absence of thermalization in a ID Bose 
gas has been confirmed by the time evol u tion f rom an 
out-of-equilibrium state (jKinoshita et al I l2006h . The 
non-equilibrium dynamics of an impurity in ID Bose 
gases, such as large density fluctuations, multiple scat- 
tering events and interaction dependent quadruple os- 
cillatio ns, have been studied in som e rece nt experi- 
ments (jFabbri et aZ.l . l2012l ; IPalzer et all . [2009). Most re- 



cently, the experimental observations of fast relaxation 
towards equilibrium of quasi-local densities, currents and 
cohe rence in an isolated strongly correlated ID Bose 
gas ( Trotzkv et all . l2012h give a precise understanding 
of quantum dynamics and correlations, including obser- 
vation of_rtiB_sp2n_dyirainics in ID two-component Bose 
gases (|Widera et aZ.1 . 12008ft . 



As remarked, quasi- ID trapped systems are created 
by tight transverse confinement. A confinement-induccd- 
resonancc (CIR) takes place when the incident channel of 
two incoming atoms and a transversally excited molecu- 
lar bound state become degenerate. It has been demon- 
strated that a CIR takes place if the 3D scattering length 
approaches the length scale of the transverse trap. The 
CIR has been used to drive a crossover from the Tonks- 
Girardeau gas with a strongly repulsive interaction to a 
super Ton ks-Girardeau gas with str ongly attractive in- 
teraction ( Haller et al I l2009l l2010bh . In particular, the 
stable highly excited gas-like phase known as the super 
Tonks-Girardeau gas has been realized in the strongly 
attractive regime of bosonic Cesium atoms. 



2. Fermi gases 

In a 3D trapped Fermi gas, the bound diatomic 
molecules exist only when the scattering length be- 
tween the atoms is positi ve under s-wave in teraction, 
i.e., when l/(k F a s ) > (|Regal et all . l2003h . In the 
limit \/(kpa s ) <C 0, the system is a weakly attrac- 
tive Fermi gas. Thus the groundsta te is a BCS pair 
state for a negative s cattering length (IChin et al 1 12004 
iGreiner et oil l2005bf ). However, in a ID two-component 
trapped Fermi gas, the scattering properties of two col- 
liding atoms are altered by the tight transverse con- 
finement. The existence of a bound state does not 
rely on the sign of the scattering length. Using ratio- 
frequency spectroscopy in an array of ID Fermi atomic 
gases trap ped within a 2D op tical lattice, Moritz and 
coworkers ( Moritz et all . 120051 ) reported this that the 
two-body bound state exists irrespective of the sign of 
scattering length, see Fig. [3(11 In this experiment, they 
further demonstrated that the bound states for nega- 
tive scattering length ca n only be stabilized by the tight 
transverse confinement ( Moritz et al ]. l2005l ). 

The strong transverse confinement of a waveguide 
not only gives rise to an effective ID s-wave interac- 
tion, but also alters the p-wave interaction in spin- 
polarized fermions. Due to the absence of s-wave inter- 
action in the spin-polarized Fermi gas, the p-wave in- 
teraction b ecomes dominant under resonant scattering 
cond itions ( Granger and Blumd . l2004lllmambekov et all . 
12010ft . In a spin-polarized Fermi gas, the angular part of 
asymptotic collision wave functions is either the spheri- 
cal harmonic Y/ = o,m=o if the scattering state is parallel to 
the quantization axis or Yj = o jm =±i if the scattering state 
is perpendicular to the quantization axis. Therefore, in 
both the 2D and 3D cases, due to the coexistence of two 
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FIG. 30 Two-body bound states in ID and 3D. In the ID 
case, confinement induced molecules exist for arbitrary sign 
of the scattering length. Whereas in the 3D case, there are no 
bound states at magnetic fields above the Feshbach resonance 
(vertical dashed line). From Moritz et al. (2005). 



collision channels and the br eakdown of degenera cy be- 
tween two collision channels ( Ticknor et al\A2004) . dou- 



blet struc tures of p-waye Fes hbach resonance have been 
observed ( Giinter et al I l2005h . However, for the ID spin- 
polarized system, only one of two collision channels is in- 
volved. Thus there is a single peak structure of p-wave 
Feshbach resonance in ID spin-polarized fermions. The 
experimen tal observation has c onfirmed such a particular 
signature (jGiinter et al ., 2005). see Fig. [3TJ By loading 
the spin-polarized Fermi atoms into a deep 3D optical lat- 
tice, one can prepare a band insulator of localized atoms 
in potential wells, which is viewed as a zero-dimensional 
(0D) system. There is no resonance feature in such a 
0D system because the p-wave scattering is completely 
inhibited. When the geometry of the gas is tuned from 
3D to 2D, the shift of the p-wave Feshbach resonance de- 
pends on the depth of the optical lattice. In contrast, ID 
Fermi gases show a further shift of resonance and give 
rise to a broadening of the loss feature. 

The ID two-component Fermi gas defined by the 
Gaudin-Yang model is an ideal system for exploring 
novel pairing mechanisms. In particular, the formation 
of a FFLO-likc state with a nonzero centre-of-mass 
momentum gives rise to a precise understanding of the 
coexistence of BCS pairs and polarizations. In 3D, the 
FFLO state occupies a tiny portion of the phase diagram 
and it is very difficult to observe in ex periments. In con- 
trast, in t he ID Gaud i n-Yan g model (jGuan et "Zl. 120071: 
iHu et all 120071: lOrsol 120071) the FFLO state becomes 
much more robust due to the band fillings of two Fermi 
seas related with one another, and it occupies major 
parts of the phase diagram, as demonstrated in Fig. [7] 
and Fig. |H1 See also the experimentally measured phase 
diagram Fig. 1321 The system has spin population imbal- 
ance caused by a difference in the number of spin-up and 
spin-down atoms. The key features of these T — phase 
diagrams have been experimentally confirmed using 
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FIG. 31 Dimension dependence of the p-wave Feshbach res- 
onance in spin-polarised Fermi atomic gases. In both (a) 3D 
and (b) 2D Fermi gases, the coexistence of two collisional 
channels (m = and \m\ = 1) give rise to a doublet feature, 
(c) In a ID Fermi gas with the spin aligned orthogonal to the 
atomic motion, the existence of only the collisional channel of 
\m\ = 1 gives rise to a single peak feature, (d) In a ID Fermi 
gas with the spin aligned orthogonal to the atomic motion, 
the existence of only the collisional channel of m — gives 
rise to another single peak feature, (e) In a 0D Fermi gas 
within a deep 3D optical lattice, there is no resonant peak 
due to the absence of all collisional channels. From Giinter et 
al. (2005). 



finite temp erature density p rofiles of trapped fermionic 
6 Li atoms ( Liao et all . 120101) . Experimental observation 
reveals that the system has a partially polarized core 
surrounding by either fully paired or fully polarized 
wings, which is in agreem ent with theoretical predictions 
by diffe rent group s (IFeiguin and Heidri ch-Meisner . 



2007 



Guan et all 120071: IHu et all 12007 



Kakashvili and Bolechl l2009t lOrsol 120071: iParish et al 



20071 ). In a 3D trapped Fermi gas, the fu lly paired core is 
surrounded by a shell of exc ess fermions ( Partridge et all 
l2006t IZwierlein et adl2006t ). 

The key method to ma p out the phase diagram in 
Fig. [32] in the experiment ( Liao et all . |2010( ) is to mea- 
sure the in situ densities of the two spin species. The ID 
spatial density profiles nifi(z) can be expressed in terms 
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FIG. 32 Experimental phase diagram as a function of the 
central tube polarization. The red diamonds and blue circles 
denote the scaled radii of the axial density difference and the 
minority state axial density, respectively. The solid lines are 
given by the TBA equation (|47l) at temperature T = 175nK. 
Experimental observati on is in reasonable agreement w i th the 
theoretical p redict i on dFeiguin and Heidrich-Meisnerl. 12007 ; 



Guan et al. 
20091 ; lOrso 



iHu et al 



iParish ~ 



20071; 
[J, l2007f ) 



Kakashvili and Bolech 



From Liao et al. 



(2010). 



of chemical potential fx = fj,o — V(z) and effective exter- 
nal field h = ho- Here /in is the chemical potential at the 
trapping centre and ho is the effective magnetic field of 
the homogeneous system. They can be obtained from the 
relations (TT4")) or the equation of state (|4"5)) . Within the 
LDA (|39[) . theoretical density profiles can be used to fit 
the experimental data obtained by inverse Abel transfor- 
mation of the radial profiles, as per the Method section 
in Liao et al. (2010). By changing the polarization, the 
threshold values of different phases in the density pro- 
files read off the phase boundaries of the phase diagram 
of Fig. [32j From the density profiles given in Fig. [33l 
we see that at low polarization below the critical value, 
the system has a partially polarized core surrounded by 
fully paired edges. At the critical polarization, almost 
the whole region is partially polarized. At high polariza- 
tion above the critical value, the system has a partially 
polarized core surrounded by fully polarized edges. 

Systems with a small number of trapped atoms 
are also experimentally feasible. Most recently, the 
ID fermionization of two distinguishable fcrmions has 
been e xperimental l y stu died by using two fermionic 6 Li 

The energy of the two-particle 
t4-) is determined by tuning the 



2012D 



atoms (jZiirn et al 
system in the state 
trapping potential barrier through which the particles 
can tunnel out. The fermionization of two distinguish- 
able fermions was identified by measuring the tunnelling 
time constants for different values of the ID interacting 
strength. It is particularly interesting that for a mag- 
netic field below the confined induced resonance two in- 
teracting fermions form a Tonks-Girardeau state whereas 
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FIG. 33 Integrated axial density profiles of an array of ID 
spin-imbalanced two-component Fermi gases for different cen- 
tral polarization P. The black circles represent the majority, 
the blue diamonds represent the minority, and the red squares 
show the difference, (a) At low P(= 0.015), the partially po- 
larized core is enclosed by the fully paired edge, (b) For in- 
creasing P(= 0.055), the partially polarized core grows and 
the fully paired edge shrinks, (c) Near P C (P = 0.10), where 
almost the entire cloud is partially polarized, (d) For the po- 
larization exceeding the critical value P C (P — 0.33), the edge 
of the cloud becomes fully polarized. From Liao et al. (2010). 
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FIG. 34 Interaction energy of two distinguishable fermions 
trapped in ID across a confined induced resonance. The blue 
line is the theoretical result of the energy shift. The blue 
points show experimental data for the interacting energy of 
the two distinguishable fermions. The black vertical line is the 
value of the magnetic field at the confined induced resonance. 
From Zurn et al. (2012). 



a super Tonks-Girardeau gas is created when the mag- 
netic field is above the resonance value, see Fig. [34] In 
this case, the two-particle super Tonks-Girardeau state 
is stable against three-body collisional losses since there 
is no third particle present. 
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VIII. CONCLUSION AND OUTLOOK 



In previous sections we have seen how results for ex- 
actly solved models of ID Fermi gases provide valu- 
able insights into a wide range of many-body phe- 
nomena including Fermi polarons, Fulde-Ferrel-Larkin- 
Ovchinnikov-like pairing, the few-body physics of trions, 
Tomonaga-Luttinger liquids, spin-charge separation, uni- 
versal contact, quantum criticality and universal scaling. 
This included discussion of the exotic many-body physics 
of ID Fcrmi-Bosc mixtures and ID multi-component in- 
teracting fermions, in particular with SU(3), 50(5) and 
SU(N) symmetries. We reviewed experimental progress 
on the realisation of ID quantum atomic gases and the 
experimental confirmation of the phase diagram of the 
Gaudin-Yang model in a ID harmonic trap. In the 
present section, we discuss some of the promising devel- 
opments and provide an outlook for future research. The 
key points we have identified are as follows. 

(a) High spin symmetries. Very recent experimen- 
tal explora t ion o f highly symmetric Mott insulators 
(ITaie et Ml, 12012ft gives insight into the Pomcranchuk 
(1950) type of cooling due to the fact that the spin de- 
grees in the Mott insulating state can hold more entropy 
than the Fermi liquid does. In this Mott insulating state, 
the entropy per site increases as the spin degrees of free- 
dom increase. It in turn leads to a temperature reduction 
through transfer of entropy from particle motion to spin 
degrees of freedom. This exotic nature is different from 
large spins in solids where the quantum fluctuations of 
large spins are weak. Such higher symmetry opens up 
further study of magnetically ordered phases in ultracold 
fermionic atoms. In particular, loading large spin ultra- 
cold fcrmionic atoms into a ID quantum wire geometry 
provides exciting opportunities to test multi-component 
TLL theory and competing supcrfluid ordering in ID in- 
teracting fermions with higher symmetries. E.g., SU(4), 
50(5) and 50(4) symmetries in spin-3/2 fermions. It is 
natural to expect that ID interacting fermions with large 
spin resulting from higher mathematical symmetries will 
greatly expand our understanding of many-body physics. 
However, comprehensive understanding of these models 
still poses theoretical and experimental challenges due 
to the more complicated grand canonical ensembles in- 
volved. 

(b) p-wave BCS pairing and synthetic gauge fields. 
The traditional s-wave BCS pairing models have vari- 
ous applications t o prob l ems in condensed matter physics 
(iDukelskv et all 12004 iLinks et all 120031: lYing et all 
2008al Tbl). nuclear physics (|Pan and Draaverl 20021) and 
ultrasmall metallic grains (Ivon Delft and Ralphl . I200I 
etc. More work along this line has b een achieved in 
the s t udy of a p + ijo p a iring model (|Dunning et all 



- confined strong-pairing and deconfined weak-pairing, 
separated by a third-order confinement-deconfinement 
quantum phase transition. Moreover, a type of electron 
pairing model with spin-orbit interactions shows that the 
pairing order parameter can always have p+ip- wave sym- 
metry regardless of the strength of pairing interactions 
( Liu et all . l201ll ). More work is required along this line. 



In particular, recent developments in the study of ul- 
tracold atoms have opened up new opportunities to sim- 
ulate synthetic external abelian and non-abelian gauge 
fields coupled to neutral atoms b y controlling atom- 
light interactions ( Lin et all . 120091) and the ultracold- 



l2010t llbanez et all . 120091: iRombouts et all . l2010h . In 



contrast to the s-wave pairing model, the p + ip pair- 
ing model has an exotic zero temperature phase dia- 
gram in terms of density and attractive coupling strength 
( Rombouts et aZj . l2010h . It shows two supcrfluid phases 



atom analog of the mesoscopic conductor (jBrantut et all 
|2012| ). An important application of this synthetic gauge 
field is the realis ation of spin-orbit coupling in degener - 
ate Fermi gases (jCheuk et all 120121: IWang et all . l2012h . 
The Fano-Feshbach resonances can be used to modify 
the strength of atomic interactions wi thout changing 
the c haracteristic range of the potential ( Williams et all . 
2012). Raman laser beams coupled to Zeeman states of 
ultracold fermionic atoms can create synthetic magnetic, 
electric fields and spin-orbit c oupling, see a brief review 
dBloch et aZ. l. l2012tlZhall2lH2l) . These new developments 
inspire further study of spin-orbit coupling, p-wave and 
d-wave pairing by exact solutions of new mathematical 
models, with clear applications in physics. 

(c) Universal Wilson ratio. The low-energy physics 
of interacting Fermi systems exhibits universal phenom- 
ena, such as TLL in ID and Fermi liquids in higher di- 
mensions. It is highly desirable to find an intrinsic con- 
nection between these low-energy theories. The Wilson 
ratio is the ratio of susceptibility to specific heat. De- 
spite these two physical quantities having different low 
temperature behaviour, the Wilson ratio is a constant 
at the fixed po int of interacting fermionic systems in 3D 
(IWilsonl . ll975l) . For noninteracting electrons, the Wilson 
ratio is unity. The value of the ratio indicates interac- 
tion effects and quantifies spin fluctuations. In contrast 
to the phcnomcnological TLL parameters, this ratio gives 
a measurable physical quantity which manifests univer- 
sal TLL physics - the effective fixed point model of the 
TLL universality class. Experimental measurement of 
the Wilson ratio of the TLL in a spin- 1/2 Hei s enber g 
ladder has been recently reported ( Ninios et all 120121) . 
It is naturally to believe that this ratio can be used to 
quantify the magnetic phases (the FFLO-like phases) in 
ID attractive Fermi gases of ultracold atoms with differ- 
ent symmetries. We remark that at the critical point it is 
a constant solely dependent on the TLL parameter and 
onset charge velocities of different states. This ratio is 
experimentally measurable with ultracold Fermi atomic 
gases trapped in ID. 

(d) Diffraction vs non-diffraction. Non-diffraction in 
the many-particle scatt ering process is a un i que feature 
of in t egrable systems (iGu and Yang . 19891 iLamacraftl . 
120121 : iMcGuirel . 11964 ISutherlandl . 120041) . In principle, 
diffractive and non-diffractivc scattering can be tuned 
via controlling inter- and intra-species scattering lengths. 
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Experimental exploration of a violation of this charac- 
teristic would yield valuable insight into understanding 
the non-diffractive form of the Bcthe ansatz wave func- 
tion, i.e., how weak violations of non-diffractive scat- 
tering still give rise to the Bethe ansatz wave func- 
tion in the asymptotic region. ID many-body sys- 
tems with a finite range potential between atoms are 
likely to be an ideal simulator for identifying the con- 
sequences of diffractive vs non-diffractive scattering. In 
this regard, across a narrow resonance, the ID effec- 
tive potential is determined by not only the scattering 
length, but also the effective range which is introduced 
by a strong energ y- dep endent scattering amplitude of two 
collid i ng atoms (ICui). l2012at [G urarie and RadzihovskyL 
120071 : iQi and Guar! l2012h . Be yond the fermioni zation 
of two distinguishable fermions ( Ziirn et all I2012T) . it is 
an immediate goal to experimentally probe three distin- 
guishable fermions in a ID harmonic trap. A possible test 
of the fundamental nature of diffraction vs non-diffraction 
for systems with three particles may ultimately lead to 
experimental tests for Yang-Baxter intcgrability in quan- 
tum many-body systems. 
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